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Aerodynamic Interference Between 


Moving Blade Rows 


NELSON H. KEMP? ann W. R. SEARS: 
Cornell Unwersity 


ABSTRACT 


The blades of a typical turbomachine move through a nonuni 
form field of flow, disturbed by the induced effects of other sta 
tionary and moving blade rows. In this paper the concepts and 
results of the theory of single thin airfoils in nonuniform motion 
are applied to calculate the resulting nonsteady effects on such 
blades. Incompressible inviscid flow is considered. The usual 
representation of blade wheels as infinite cascades of two-dimen 
sional airfoils is adopted 

In this investigation, attention is confined to the basic case of 
a single stator row followed by a single moving motor row. As 
a first approximation, it is assumed that, to calculate the non- 
steady effects at any given blade, it is permissible to neglect the 
nonsteady parts of the circulation of all other blades. However, 
for the rotor blades, the effects of their passage through the vortex 
wakes shed by stator blades are also investigated. It is found 
that, with this approximation, the only results of thin-airfoil 
theory required can be obtained by a simple generalization of an 
earlier study of an airfoil flying through a sinusoidal gust pattern 
Expressions are worked out for the unsteady components of lift 
and moment on stator and rotor blades. 


Numerical calculations are carried out for some examples typi 
cal of compressor cascades. The results show fluctuating lift 
values with magnitudes as great as 18 per cent of the steady lift 
The effects of stator wakes on the unsteady lift of rotor blades 
are found to be comparable to the induced effects of the stator 
blades themselves 


The rate of energy transfer to the wake pattern of the rotor 
Stator stage has also been calculated within the approximations 
of this theory, but the details are not presented here. For the 
compressor stages of the present numerical examples, it is found 
that this rate of energy transfer is a very small fraction of the 


power required to turn the rotor 


Presented at the Aerodynamics Session, Twenty-First Annual 
Meeting, IAS, New York, January 26-29, 1953 
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Force, Air Research and Development Command, Office of 
Scientific Research. 
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1953. Now Aeronautical Engineer, Research Department, United 
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List OF SYMBOLS 


axial distance between centers of stator and rotor 
rows (Fig. 1) 

axial distance between trailing edge of stator row 
and leading edge of rotor row (Fig. 1) 

blade semichord (Fig. 1) 

pitch of blades in direction of rotation (Fig. 1) 

parameter giving relative orientation of rotor 
and stator blades (Figs. 3b and 4b) 

time 

horizontal and vertical perturbation velocity com 
ponents, respectively 

coordinate along blade from center, positive toward 
trailing edge 

coordinate perpendicular to x, positive up 

x+y 

coefficient of steady bound-vortex distribution, 
Eqs. (27) and (38) 

Theodorsen function, K\(iw)/|Ko(tw) + Ayi(iw)| 

energy coefficient, E/pV% 

steady lift coefficient, »L./pV*% 

average energy per unit length per unit width in 


wake 





STATOR 


Notation and configuration of elementary turbo 


machine stage 
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F = coeflicient of quasi-steady circulation, Eqs. (3), 
(48), and (53) 

G = coefficient of velocity induced at blade, Eqs. (37 
and (41) 

Il = steady bound-vortex distribution function, Eqs 
(31) and (42) 

a. = Bessel function of the first kind 

J = J,-—1); 

i ¢ = modified Bessel function of the second kind 

K,, Ku = lift and moment functions defined in Eqs. (16) and 
(17) 

E = lift, positive up 

M = moment about the center of the blade, positive 
clockwise 

R = real part 

S(w) = fiw[K (iw) + Ki(iw)]! 

U = rotational speed of rotor, positive down 

V = speed of flow relative to blade 

a = angle of stagger of blades (Tig. 1) 

B = angle between stator and rotor blades, a, + a 

Y = strength of bound vortex distribution, positive 
clockwise 

€ = strength of wake vortex distribution, positive 
clockwise 

¢ = complex coordinate 

t,n = Cartesian coordinates 

r = uc/V 

m = complex frequency of chordwise velocity distribu 
tion, Eq. (13) 

v = circular frequency of motion in time, v 2rl'/d,, 
vy = 2rU/ds 

p = air density 

o = solidity of row, 2c/d 

w = reduced frequency, vc/ | 

r = circulation, positive clockwise 

> = summation from » = 1 ton = 4 

Nn 

7 = summation from ” = ton fo, omitting 

ad n = (term 

Postsubscripts 

Ss = stator 

r = rotor 

k, m,n = summation indices 

n = nth blade 


0 = steady quantities 

q = quasi-steady quantities 

l = apparent-mass contributions 
2 = wake contributions 

t = time-dependent quantitics 


s = stator 

} = rotor 

rl = effects at rotor due to motion of steady stator cit 
culation 

r2 = effects at rotor due to stator wakes 


INTRODUCTION 


acne HAVE LONG RECOGNIZED that the airfoils 
arranged in rotating and stationary cascades of 
axial-flow machinery do not operate in steady flow. 
A typical airfoil of such a cascade acts in the disturb 
ance field of all the other airfoils of its own blade row 
and of near-by blade rows. Since these rows rotate 


relative to one another, each airfoil acts in a time-vary 


FIC AL 
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1953 
ing disturbance field; thus, its circulation varies with 
time, it experiences nonsteady forces, and it produces a 
varying disturbance field itself. The periodic changes 
of its circulation involve the production of a wake of 
vortices shed into the fluid at its trailing edge, as wil] 
be discussed below, and this phenomenon represents 
an undesirable energy exchange. The probable mag 
nitude of this energy exchange, in a typical turbojet 
engine, has been discussed briefly by Wu,! based largely 
on the work of Keller.” 

The intention here is to attack the problem of mutual 
interference of blades in cascades, due to relative mo 
tion of successive blade rows, on the basis of the theory 
of thin airfoils in unsteady motion. The objectives of 
this calculation will be to evaluate the magnitudes of 
the unsteady forces on the blades and the rate of energy 
transfer to wake vortices and to determine how these 
are affected by the major geometrical and aerodynamic 
parameters. 

To make any real progress in this complex problem, a 
number of simplifying approximations seems to be re 
quired from the start. Thus, the fluid will be assumed 
incompressible, and the usual representation of a blade 
wheel by an infinite cascade of two-dimensional air 
foils will be adopted. The airfoil theory to be used 
here is based, furthermore, upon the approximation of 
an inviscid fluid. This will not eliminate the possibil 
ity of accounting for some of the effects of viscous 
wakes in fact, this may be one of the important ap 
plications of the methods of this study. Nevertheless, 
in the present report all viscous effects are neglected 
to avoid needless complications in the development of 
the methods and in the belief that the extension to ac 
count for passage of the blades through the viscous 
wakes of blades upstream can be accomplished later 
without serious difficulty. 

Since there is no usable nonsteady-airfoil theory 
without the approximation of thin slightly cambered 
lightly loaded 


turning angles through the blade rows seems to be an 


(1.e., airfoils, the restriction to small 


inevitable one. Thus the study is limited to cascades 
of blades whose chord lines lie approximately in the di 
rection of the relative stream vector. It is therefore 
more directly applicable to typical compressor stages 
than to turbines. 

Finally, for the sake of simplicity in this initial attack 
on the problem, attention will be confined to the ele 
mentary case of a single stator cascade followed by a 
This figure 
thin 


single rotor cascade, as sketched in Fig. 1. 
represents two infinite cascades of nearly plane, 
airfoils of unequal chords, gaps, and incidence and 
The angle a, rep 
The 


rotor row moves (downward) relative to the stator with 


spaced an arbitrary distance apart. 
resents the swirl of the entering basic stream. 
velocity Ll’, and a, denotes the apparent swirl of the 
basic stream relative to the rotor, as indicated in Fig. 
1. It is to be understood that, in the usual fashion o! 


airfoil theory, all aerodynamic effects of the blade rows 
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RODYNAMIC INTERFERENCE 


AE 


be represented as induced velocities super 


The lightly loaded airfoil 


are to 
imposed on this basic flow. 
approximation means that such induced velocities are 
always small compared with |’, or V7. 


UNSTEADY AIRFOIL THEORY 


Before attacking this problem, it may be helpful to 
review briefly the major assumptions of the theory of 
thin airfoils in unsteady motion and some of its con 
clusions that are applicable here. For a more complete 
account, the reader is referred to references 3 and 4. 

In this theory the airfoil 1s visualized as being re 
placed by a plane vortex sheet (Fig. 2) whose strength 
y(x, t) is such, at every instant, as to cancel the relative 
velocity component normal to the airioil surface, thus 
satisfying the boundary condition there. The Kutta 
Joukowski condition of finite velocity at the trailing 
edge is imposed at every instant. According to the 
theorem that the total circulation about the whole sys 
tem is invariant, variations of the airfoil circulation 
l(t 


vortices 


must be accompanied by the shedding of free 


the so-called “‘starting vortices’ at the 


trailing edge. Since no forces act on these shed vortices, 


they are carried away by the flow 1.e., in first-order 


theory, they are carried rearward with the stream 


velocity V. 
Thus, the determination of flow about 
nonsteady motion or subjected to nonsteady flow con 


an airfoil in 


ditions is similar to the same problem in steady-airfoil 
theory, except that the transient velocity field of the 
trailing wake of free vortices must be accounted for in 
That is, the wake induces velocity 


calculating y(x, ¢). 
components normal to the airfoil, which are superim 
posed upon other relative normal components because 
of the airfoil’s camber, incidence, and motion and due 
the All of these com 

j 


involved in the determination of y(x, ¢), 


to nonuniformities of stream. 
ponents are 
as mentioned above. 

The only additional complication caused by the 
transient nature of the flow arises in the calculation o: 


pressures and forces alter y(x, /) is known. It will be 
recalled that the formula for pressure in terms of veloc 
ity components is different in steady and nonsteady 
Che familiar Joukowski theorem, which equates 
the hit is 


which 


flows. 
the lift to pl, is no longer valid. Instead, 


given by more complicated expressions, are 


worked out* * by separating the lift into three recog 


terms, 
At 


is the lift that would be calculated for the 


nizabl« 


»L(t) + Lil L(t 


where ,L(t 
instantaneous motion and flow conditions, neglecting 
entirely the presence of wake vortices (this is called 
,L(t) is the lit corresponding to 
variations of the “apparent mass’’ of the airfoil; and 
Lit 


The moment «bout the airfoil center is given simi 


the quast steady lift 
contains all wake effects. 


larly by three inalogous terms, 


SETWEEN MOVING BLADE ROWS DST 


| 
Yb) € (x, ¢) 
Raper 
| x 


Airfoil and wake in unsteady motion represented 
is vortex sheets 





Fic. 2 


W(t W(t) + MWe Mit (2 


LL, .M, and .M are 
given in references 3 and 4, and some of these will be 
reviewed the 
quasi-steady lift and moment, ,Z and ,.\/, are available 


Formulas for calculation of iL, 


below. Formulas for determination of 


aerodynamics textbooks and standard reference 


L and ».\/ require knowledge of the 


ill 


works. Only 
wake configuration and therefore depend on the history 
of the motion. 

Fortunately, considerable simplification of this whole 
matter is possible if the motion is sinusoidal in time. 
Then as can easily be found,’ all three parts of the lift 
and moment are sinusoidal, the wake vortex strength is 
sinusoidal in both x and ¢, and ».L and »./ can be cal 
culated directly in terms of ,L. 

Because the expressions for lift and motion for a cer 
tain case of sinusoidal motion will be used extensively 
The 


familiar complex notation will be used to represent 


in the present analysis, they will be reviewed now. 


sinusoidally varying terms; it must be remembered 
that, wherever physical quantities are given complex 
expressions, it is the real part that is to be taken unless 
otherwise specified. 

If the motion of the airfoil of Fig. 2 is such that its 


quasi-steady circulation is given by 


A(t Fe” (3) 

then reference 3 shows that 
L(t) + oL(t L(t) C(w) (4) 
where w denotes the “reduced frequency,” ve I, « 
being the semichord, and C(w) is the Theodorsen 


function,® which plays such an important part in wing 
flutter analysis. The Theodorsen function is a com 
plex function whose real and imaginary parts are tabu 
lated in various references. One analytic formula for 
it is 

[A (tw) + Ay(tw) | (5 


( ( Ww K 1\ Iw) 


where A, and A, are modified Bessel functions of the 
second kind.’ 

Since the expression for ,L(t) is not needed in the 
present paper, it will not be given. Instead, the results 
of the special case of sinusoidal motion of interest for 
the purposes of this paper will now be summarized. 
Reference 4 shows that, if an airfoil operates in a sinu 
soidal gust pattern such that the relative upwash at 


the airfoil is given by 
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v(x, 1) = nel”! — *//) (6) 


then the total nonsteady lift and moment about the 
center are given by 

L(t) we | 
M(t) = 


2m pc Viy.S(w)e 
»L(t)(1/2)e 


Il 


| 


7) 
{ (4 


where c is the semichord (not the chord, as in reference 
1) and 


S(w) = {tw[K,(tw) + Ki(tw)]}~! (S) 
|The identification of the function S(w) with that given 
in reference 4 was made in reference 11.] 

Inspection of the original form of ,L(¢) in Eq. (14) of 
reference 4, shows that, since the coefficient of C(w) 
in the expression for ,L(¢) must be ,L(t) [see Eq. (4)], 
the corresponding quasi-steady circulation is given by 


a(t) = 2mcv,J (we (9) 


where /(z) is a combination of Bessel functions defined 


as 


J(z) = J,(z) — tJ,(2) (10) 


From ,I'\(¢) the strength of the vortex distribution in the 
wake, e(x, /), can be obtained. Reference 3 shows that, 
for general sinusoidal motion in which ,I'(¢) is of the 


form of Eq. (3), the wake strength is 


a) Fe’” — x/V) F ; 
ex, t) = sme ee = — iwS(w)e? ~ *//) 
c[K,(iw) + Ky(iw) | Cc 

(11) 


Applying this to Eq. (9), the wake strength for the 
sinusoidal gust of Eq. (6) becomes 


e(x, t) = ee arr? (12) 


—2rv,J(w)iwS(w)e 
These results were derived in reference 4 only for the 
an airfoil flying through a stationary 
It will be shown in the present 


special case of 
sinusoidal gust pattern. 
report that, by means of a simple generalization, they 
provide just the formulas needed to calculate the non- 
steady aerodynamic quantities in the present rotor- 
stator problem. The generalization required has al- 
ready been pointed out in reference 11; it consists of 
replacing the relative upwash expression (6) by 


tut tux /V 


v(x, t) = v,e e€ (13) 


where v is real but » may be an arbitrary complex num- 
ber. As will be seen farther below, the induced velocity 
fields in which the rotor and stator blades operate can 
be expressed in this special form. 

The extension of reference 4 to this more general case 
is not difficult; it is only necessary to change the argu- 
ment of the Bessel functions J, in Eqs. (13), (14), and 
(15) of that reference to the complex parameter \ = 
uc/ V. 
Bessel functions, the expressions for lift and moment 


Upon application of the recurrence relation for 


are then immediately found to be 
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L(t) = 2apcVu,Kz(w, r\)e” (14 
M(t) = wpc?Vv,Ku(o, Ae” (15 

where 
Kz(w, A+) = J(A)C(w) + 4(w/A) Si (A) (16 


Ky(w, +) = J(A)[C(w) — 1] + (w/A)J,(A) 4 
(2 A)(CI — w/d)Ji(X) (17 
(In reference 11, A, was denoted by A.) 
Again inspection of Eqs. (14) and (16) shows that 
the corresponding quasi-steady circulation for this case 
must be 


al’ (t) = 2mcv,J (re (18 


so that the wake vortex strength follows from Eq. (11 
as 


vil x/V) 


e(x, 4) = —22v,J(A)iwS(w)e’ (19) 


The results expressed in Eqs. (7), (9), (12), (14), 
(15), (18), and (19) are those that will be used in the 
present attack on the two-cascade problem set up in 
the preceding section. 


APPLICATION TO STATOR-ROTOR PROBLEM 


In basic idea, the extension of thin airfoil theory to 
the stator-rotor interference problem of the present 
paper is fairly obvious: One must consider each airfoil 
to be acting in a velocity field induced by (a) its own 
wake, (b) the (variable) bound vortices of fellow mem- 
bers of its own blade row, (c) their wakes, (d) the 
(variable) bound vortices of members of the other blade 
row, and (e) their wakes. 

It is not surprising, however, that this approach leads 
to rather discouraging mathematical complexity—in 
fact, to a pair of intractable, simultaneous integral 
equations for the stator and rotor bound-vortex distri- 
butions. 
scheme, based upon an idea of successive approxima- 


It seems more profitable to adopt a different 


tions, which will now be outlined as follows: 

(1) Let it be assumed, subject to later verification, 
that the entire unsteady effect on the circulation of any 
blade is small compared with the steady circulation 
carried by the blade. Symbolically, this means that, 
for any airfoil 


reh= 1 + Th: ad iF 


where ,I° is the steady circulation of the blade, con- 
sidering its camber, incidence, rotational speed (for 
rotor blades), and the steady induction effects of its 
own cascade and the other cascade. In other words, 
ol’ is the circulation ordinarily estimated by the de- 
signer, neglecting all nonsteady effects. Thus, ,I'(/ 
contains all the effects due to the relative motion of the 
cascades. 

(2) As a first approximation, to calculate the non 
steady effects on a typical stator blade, such as ,I'*(¢), 
let it be assumed that the only significa 1t contribution 
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arises from ,l'’. This means that the effects of the 
ynsteady terms ,I’(f) of rotor circulation, as well as the 
ynsteady parts of the circulation of all other stator 
blades and also all vortex wakes except that shed by the 
blade itself, are neglected. (Presumably this approxi- 
mation breaks down for stator cascades of extremely 
large solidity, in which the adjacent blades are so close 
that their nonsteady circulation components induce 
large effects. ) 

(3) In the same approximation, to calculate the un- 
steady effects on a typical rotor blade, let it be assumed 
that the significant contribution arises from ,I°’. For 
the rotor blades, however, one additional effect will 
also be investigated—-viz., the effect of their passage 
through the vortex wakes shed by the stator blades. 
To be sure, this is, in the present system of approxima- 
tion, a sort of second-order contribution, since it is due 
to the unsteady part of the stator circulation. But 
the rotor blades pass directly through the stator wakes 
encounter appreciable induced 
For this reason it would not 


and therefore 
velocities due to them. 
be surprising if this ‘‘secondary”’ effect were found to be 
comparable in magnitude to the “primary” effects of 


may 


the steady circulation components. 

In principle, the approximate theory outlined here 
might be the first step of a series of successive approxi- 
mations—1.e., the second '*(é) 
might include the effects both of ,I’’ and of the first- 
approximation unsteady However, the 
mathematical complexities of the first approximation, 
as will be seen shortly, suggest that a second approxi- 
mation might be extremely difficult. At any rate, the 
first step will permit approximate evaluation of some 
useful engineering data, the absolute accuracy of which 
may subsequently be checked by experiment or by closer 


approximation to 


quantities. 


approximations. 


CALCULATION OF PERIODIC QUANTITIES 


Induced Velocity Due to ,V 


In this subsection the velocity at a rotor or stator 
blade induced by the steady circulation distribution of 
the blades of the other row, which is time-dependent 
because of the motion of these blades relative to the 
blade in question, will be found. 

The upwash velocity on the x-axis and normal to 
that axis, produced by equal vortices of strength I’ 
located at any .V complex points ¢,,, is given by* 


\ 


—(I'/2r) > 1/(x — ¢,) 


n l 


v= (20) 
If this equation is applied to the induction effect of the 
rotor on a stator blade and the rotor circulation is tem- 
porarily thought to be made up of single point vortices, 
as shown in Fig. 3a, then the induced velocity at the 


Stator may be written as 


* As mentioned above, the real part is implied 
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Fic. 3. Diagram used in computing velocity induced at stator 
blade 
I’ So 
eo 9) - M4 . ta (21) 
27 n=—o X, — (E + ty + tnd,e™ 


where x, is measured parallel to the stator blade from 
its center and £ and 7 are defined in Fig. 3a. 

The summation in Eq. (21) may be performed using 
the well-known relation’ 


+- = 


> 1/(2 — ine) = 


I'’e ~1Q@s ,; we las 1 . ™ 
coth (& in — X;) (22) 
2d, oo 


This form of the velocity induced by an infinite row of 
equal, equally spaced vortices is well known but is not 
It can, however, be 


coth z 


which gives 


so useful in the present case. 


transformed into a convenient form as follows: If 


R(s) > 0, 
coth z = | 
Moreover, when the line of vortices (Fig. 3a) does not 


intersect the stator blade, which is true in the present 
case, it is easily verified that 


R{le '(& + in — x,)] > 0 


and therefore Eq. (23) is applicable to Eq. (22). The 
resulting expression for v* is 
I’’e ~ tas j 
7 = 41+2> exp X 
2d, | m 
2om ) 
— eo "NE FH OY — Xp) ( (24) 


d, 


If, now, I’ is replaced by a distribution of vortices, 


oy’ (x,) dx, along a line representing a typical rotor 
blade, as shown in Fig. 3b, Eq. (24) need only be inte 


grated over x, from —c, to c,. The resulting velocity 


+t Since summation signs will many times in this 


analysis, it will be convenient at the outset to adopt the following 
: 1 


conventions: will denote > , while re ’ will denote > + 
1 


appear 


nN n n ’ 


© 
» . If any other type of summation is used, it will be ex- 
” 1 


plicitly indicated. 
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v’ is the induced velocity at the stator blade due to the 
whole row of rotor blades in a given position. In 
order to perform this integration, the complex coordi- 
nate — + 7m must be written in terms of the geometric 
parameters and of x,, which is measured along the rotor 


blade. 


27m 


ar’ las 
= od 1 oo 3 > exp | - (b + th, 


2d, m 


where ,l’ = i oF (Xr) Ey. 


ad, 
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With the aid of Fig. 3b, it is seen that 
E+ in = (64 the + x,e°% 7% 95 


where b is the spacing of the two blade rows and h, is de. 
fined in Fig. 3b, positive when measured upward. Use 
of Eq. (25) and replacement of I’ by ,y’(x,) dx, yields 


upon integration, | 


ex, dx, ( (6 


2mm 


~~ | - d, 


To eliminate the integral from Eq. (26), it is convenient to assume Oy’ in the Glauert form 


oY = 


where x, = c, cos 6. 


2V,[A,’(1 — cos 6)/sin 6 + 2 DB A,’ sin n@| yi 


n 


In this form, ,I’ and ,j/’, the moment about the center, are given by 


ak’ = Bac V4S + .A,’) Ig 


oM” = 


2c,7pV,7(A,’ — A2’) 29 


Furthermore, the integral appearing in Eq. (26), by virtue of some trigonometric identities, becomes 


we oy (x7) 2am 1 f{? , ee oe, 
/ 17 XP | - cx, dx, = H,! = i+ > . cos n0 exp [— mma,e'* cos 6| dé 
I E d, us 0 Nn Ao i A, 


~Cr o 


where a, is written for the rotor solidity 2c, d,. 


Since 


tr t{( m/e) ar} 


—71a,e'° es ima, 


the typical integral in the expression for //,,” is 


l ” = 
exp | —t1o,me il(=/2 
WT Jo 


The well-known relation J/,(—2z) 


“"l' cos n6é dé 


—1)"J,(z) com- 


Il 


bined with"? 


. = ” Samed 
ee? ©" cos nOdd = 
Tv 0 


J,(2) 


leads to 


l <r 
e 8% eos n6 dd = (—i)"J,(z) (30) 
wT Jo 
With the result of Eq. (30), /T,," becomes 
HH,” = Jo{ amor!) — awhy 4 
yy Sent — Ant’ ,, i[(x/2) — a ‘ 
> (-1)" ; _ Srirmae i ih Me 
- A, + A, 


In terms of this notation, Eq. (26) becomes 


ale tats 
vs = J, +2>0 H,,’ exp X 
od. { m 


2am : ; 
- (6 + th, — x,e7 '*) 
d, 


a | 
~— 


*(x5, ¢), It 
is only necessary to insert the time factor by assuming 


In order to complete the expression for v 


that the rotor row, hitherto taken as stationary, is in 
It is convenient to take the positions of the 
0), for which the 


motion. 
zero-th rotor and stator blades (n = 
calculations are done, to be as shown in Fig. 1 at ¢ = 0 
The trailing edge of the zero-th rotor blade just touches 


a line drawn through the zero-th stator blade. In this 
position 

h, = —btan a, — c, sin B/ cos a, 
where 8 = a, + a,;. Since the rotor moves downward 


at a speed Ll’, h, increases negatively with time as Ut. 
Thus, 
(33 


b+ th, = 6 —2(btan a, + c,sin B/cos a.) — it 


Reference to the velocity triangle of Fig. 1 gives the 
following relations: 


V,sinB = Ucosa,; V;sina, = U — V,sina,} 
V, cos a, = V, cos a, { 
34 
so that 
b+ th, = b{1 + i(tan a, — U/V, cos a,)] — 
1Uc,/V, —4Ut (35 


Inserting Eq. (35) in Eq. (32), the final expression for 
v* is obtained. 
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AERODYNAMIC INTERFERENCE 


r(x, ) = + —— ) G,,* exp X 


2d, 2 re s osm 
™ Sie. tyvgml 
Xe € 


2am 
e 
d, 


36) 


where 


od, 
Gt ae exp Bx mmo, X 
d, 


k ( ; 1U ) il 
1 +2ztan a, — — — — 
C, V, cos a, ] 


o,, the stator solidity, and », = 


| (37) 


with 2c, d, = 24 U/d,, 
the stator frequency, which is the circular frequency at 
which the rotor blades move past a stator blade. 

The advantage of the series form of v* over the hyper- 
bolic-cotangent expression is now clear. The vari- 
ables x,, Xs, and ¢ were separated so that the integration 
over x, was easily done; furthermore, in each term of 
the series, x, and ¢ appear in the same form as x and f 
do in Eq. (13). Since the problem deals with thin air- 
foils and is therefore linear, the results of Eqs. (14), 
(15), (18), and (19) for the aerodynamic quantities may 
be applied directly to each term of v* and then summed 
to give the total effect of all terms. This procedure 
will be followed in the next section. 


It should be noted that the time-independent term 
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(b) 


Diagram used in computing velocity induced at rotor 
blade 


(a) 


tities by the turbomachinery designer. If the rotor 
and stator rows are far apart (b ~ ©), only this term 
remains. In order that this steady-flow effect is not 
inadvertently included twice, the time-independent 
term in v* will be omitted hereafter. 

The calculation of the velocity v’' induced at a rotor 
blade by the motion of the steady-stator circulation 
exactly and can be carried out 
Besides the interchange of sub- 


ol’ parallels that of v° 
with the aid of Fig. 4. 
scripts on the geometric parameters, — a, must be re- 
placed by a, and a, by —a, to obtain the expressions 
for v’' from those for 7°. With these replacements and 
the introduction of x, = c, cos 8 

oyY* = 2V,[A,*(1 — cos@)/sin@ + 2 y ¥ A, sin n@] (38) 

mn 


I = 2ec,V({A,' + A;") 1 


in v* is precisely that which arises from the steady ue ee 2 V.AA.! 12°) (39) 
: a ‘ of = 29C,” s\Ag — Ao’) 
induction effect of the rotor row on a stator blade and , 
is usually taken into account in the steady-flow quan- the final expression for v’'(x,, ¢) is found, 
— Ie’ bg 2am . 
0 0 : 2 mn _ 
o"(x,, t) = 4+ > Gn’ exp { — fhe (40) 
2d. 2 rc, m d, 
Here 
, —1o/d, ; {—1o,d, ‘ 1U iu |) 
G, = e'"H1,* exp 1+itana,—- = -— (41) 
d, d, - V. cos a, Vv, js 
§ { i[(r/2 as} *\n Ant" as A, " { i[(3/2) as) ) 
IT, = J,\ro,me t+ » i (—1) : : J, {mame ; (42) 
n A, + A;' 
and the rotor frequency »v, is given by v, = 2rl’/d,. 


The velocity at the rotor just calculated is given the 
double superscript because, as explained in the preced- 
ing section, two types of nonsteady effects on the rotors 
will be investigated. The first kind, due to the motion 
of the steady stator circulation past the rotor, will be 
denoted by the superscript r1; the velocity of Eq. (40) 
The second kind—1.e., the effects of 
will be given the super- 


is one of these. 
the stator wakes at the rotor 
script 72. These latter are discussed farther below. 

Again, the time-independent term in Eq. (40) will 
be omitted from further consideration here, its effects 
already having been accounted for in the steady-flow 
quantities. 


Now that v* and v’! are available, the application to 
them of the formulas developed in connection with the 
relative velocity of Eq. (13) can be made to calculate 
the nonsteady lift and moment and the properties of 
the vortex wakes. This will be done for the stator in 
the next subsection and for the rotor in a later subsec- 


tion. 


Nonsteady Lift and Moment of Stator 


As has been mentioned above, the similarity of the 
form of each term of v*, Eq. (36), to that of v in Eq. 
(13) leads immediately to the use of Eqs. (14) and (15) 
to obtain the nonsteady lift and moment on a stator 








blade. Furthermore, the linearity of the problem 
permits one to apply these equations term by term to 
Eq. (36). 

The argument of the exponential involving x, in Eq. 
(36) can be written 


> 


(2amx,/d,)e '™ —i(2amV,/d,)(x,/V,)é'*’” al 
(43) 
Let 
Ms (24 V,/d,)e"/? ~— a!) 14 
(w/2) ats | \ } 
h, = (2xc,/d,)e f 


Then application of Eqs. (14) and (15) to Eq. (36) gives 


1°) ol” et ns se — - 
oe ae z GK (wm, A,m)e”’ (45) 
gh aa : m 
W(t) r’ | 
te 0 . - y 
— eee > — Gn'Ku(weam, Aym)e"" (46) 
ye Cs an ‘mm 


Here the steady-stator lift ,['pV, is denoted by ,L’, 
and w, is the stator reduced frequency »,¢,/ V’,. 

Furthermore, use of Eq. (18) gives the quasi-steady 
circulation as 


(47) 


. ygmt 
I” Z, F,, se" m 


m 


A(t) = 


where 
F,,' =G,,°J (Am) (48) 
This expression for ,I'* will be needed later. 

Eqs. (45) and (46) are to be used to evaluate, within 
the approximations of the present theory, the non 
steady oscillatory parts of the lift and moment on a 
stator blade. In addition to the various geometrical 
parameters characterizing the stage, this calculation 
requires knowledge of the coefficients A,’ of the rotor’s 
steady circulation distribution. In actual use, of 
course, one will approximate to this distribution by a 
i.e., one will terminate the 


small number of terms 


series over ” in Eq. (31). Moreover, one must decide 
how many harmonics are of interest in the lift and 
moment formulas and will retain only that number of 
terms in the summations over m in Eqs. (45) and (46). 
This process will be made clearer by means of some ex- 
amples in a later section of this paper. 


Nonsteady Lift and Moment of Rotor 


The lift and moment of the rotor blade due to the 
motion of the steady-stator circulation ,I'*, which are 
denoted by the superscript 71, can be found from v’! 
of Eq. (40) exactly as the analogous quantities for the 
stator were found above. 

The argument appearing in Eq. (40), —(24mx, d,) X 


[ > } 


e'” can be written —1(2mmV,/d,)(x,/V,Je *"’* i 
Let 
ur = (24V,/d,)e **’* - 9 
h, = (2ac,/d,)e~*!"*/? oe (49) 


Then, from Eqs. (14) and (15), 


- "— SS Gn'K (wm, X,m)je 5() 
E,  S 
Mf) ali. lo 1K X 
= tm Ay(w,m, A,;m)e”” [- 
2,L'c, I" 4 . ' 


Here the steady-rotor lift ,I'’pV, is denoted by ,L’, and 
w, is the rotor reduced frequency »,c,, V,. 


From Eq. (18), the quasi-steady circulation is 


T"(0) rs >» F,, reivrmt 32 
where 
fy G,,’J(\,m) 53 


The practical calculation of ,L7! and ,4/’' is analogous 


to that of ,L* and ,.\/°, described above —i.e., it requires 
knowledge of the steady stator circulation-distribution 
coefficients A,,*, in addition to geometrical quantities 
Again, a decision must be made concerning the number 
of harmonics required. 

As mentioned above, the additional effects of the vor 
tex makes shed by stator blades are also taken into ac 
count in calculating conditions at the rotor blades 
These give rise to the contributions identified here by 
the superscript 72. The details of this calculation ar 
not reproduced here because of the necessity of con 
serving space, the calculations being somewhat lengthy). 
The detailed calculation may be found in reference 13 

To determine these r2 quantities, it is necessary first 
to determine the strength of the stator vortex wakes, 
which is easily done by means of Eqs. (19), (36), and 
(48), and then to work out the phase relationships be- 
tween the successive stator wakes. This second item 
requires a study of the phase relationships existing 
within the stator row itself, the results of which are of 
general interest and will be of use later on. It is found 
that, if /,*(¢) is any function of time associated with 
the nth stator blade, such as its lift or circulation, it 
must be that 


fran'(t) = f,*lt + (Rd,/U)] o4 


Eq. (54) establishes not only the periodicity of effects 
in the stator row but also, because it applies to the 
stator circulation, the periodicity and phase relation 
ships between various stator wakes. 

[It may be mentioned here in passing that the analo- 


gous general relationship applying to any function 


¢ 
f,’(t) associated with the nth rotor blade is 


Tare) = Pa |: — (kd, U)] 59 
This result is needed in determining the energy transfer 
to the wake. ] 

The next step is to calculate the velocity components 
induced by the infinite row of sinusoidally varying vortex 
This 
is greatly simplified if the vortex sheets can be assumed 
+ © 
the stator trailing edges as they actually do 


sheets which trails back from the stator blades. 


instead of originating at 
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approximation is made; it is probably a good one, 
since the portions of stator wakes close to the rotor con- 
tribute the largest part of the velocity, while the con 
tributions of portions farther away, where the error in 
the approximation lies, are probably negligible. With 
this simplification the calculation of the velocity com 
ponent v”? normal to the rotor blades at the position of 
a blade is simple except for the geomettical complica 
tions. The distortion of the stator wakes due to the 
passage through them of the rotor blades will be ig 
nored, as is consistent with the approximations of thin 
airfoil theory. 

The velocity component v” for the 
zero-th rotor blade during the time interval 0 < ¢ < 


is calculated 


L(t) eis, SER ely ; bl’ 
- >> exp immo, ( : 
Py mk. Us COS Qs V.c, COS a, 
»M™*(t)/2,1'¢ gl.’ *(8) / 4,5," (58) 


From Eq. (9), the corresponding quasi-steady circu 


lation is 


a) SE", « kt 
MW’ ?= Yodo’ 2nc0mJ (wrk)e” 
” A 


and v,,, can be found by comparing Eqs. (56) and (57) 


with Eqs. (6) and (7). 


These formulas complete the calculation of unsteady 
lift and moment on stator and rotor blades. For prac- 
tical computation of ,L’? and ,\/", the same decision 
as mentioned previously, regarding the number of 
Here both stator and rotor 


More- 


harmonics, must be made. 
circulation coefficients A ,,* and A, are required. 
over, there is a summation over m in Eq. (57), and 
judgment will be required regarding the number of 
terms of this series to be carried in computing each 


harmonic. 


CALCULATION OF WAKE ENERGY 


The wake of trailing vortices behind any airfoil in 
unsteady flow involves motion of the fluid far behind 
the airfoil and corresponding kinetic energy that the 
passage of the air foil has produced in the flow. As 
part of the present investigation, the energy left in this 
flow behind the elementary turbomachine stage (stator 
and rotor) has been calculated. This calculation is 
given in reference 13 and will be the subject of a sub- 
sequent short paper. For completeness, some remarks 
about the numerical values obtained for wake energies 


are included in the closing section below. 


d,/U and, finally, is expressed in complex Fourier 
series form. The result is thus obtained in the form 


v(x, 1) = DO do'vn, exp [iv,k(t — x,/V,)] 
A 


m 


(56) 


where v,»,; 1S a certain complex geometrical coefficient 
and x, is the coordinate measured along this blade 
(Fig. 1). Expressions for the aerodynamic quantities 
associated with v” are now found immediately from 
Eq. (7), since each term of v”” is of exactly the same 
form as Eq. (6) and the linearity of the present problem 
|The Fourier expansion men 
in order to throw this 


permits superposition. 
tioned above was carried out 
velocity into the form where Eq. (7) would apply. | 
The lift and moment are, finally, 


Uv 
4 .) + 1w,k x 
F,, .S(@.m)S(w,k)(V.d,/ Vid.) (R/m) i 
1 + cot? B[1 — (V,d,/ V,d,)(k/m) sec Bl? 


D7) 


NUMERICAL CALCULATIONS 


Convergence of Series 


Since the theoretical results presented in this paper 
are given in the form of infinite series, the first question 
to be answered when numerical calculations are under 
taken is: How many terms of the series must be re 
tained? i.e., how rapid is the convergence? How 
ever, the successive terms in ,L*, ,\/*, ,L", and ,/"! 
and in the summations over & in ,L” and ,.\/"* have a 
special significance: they are the successive harmonics 
of the forces and moments calculated. 
the form in which these results have been obtained 1s 
that these periodic quantities are explicitly obtained 
Thus, each 


The virtue of 


as sums of these successive harmonics. 
term will be of interest in itself, particularly in studies 
of blade vibration. The question of rapidity of con- 
vergence of these series is therefore simply the question 
of how rapidly the higher harmonics diminish 1n size. 
On the other hand, the various harmonics of ,L” 
and ,M”™ involve infinite sums over that the 
practical problem of approximating to infinite series by 
finite sums does arise here (and also in all the summa- 


m, So 


tions in the wake-energy expressions). 
Investigation of the expressions for ,L*, ,.\/*, .L 

:M"' shows that, when allowance is made for the ex- 

pontential nature of /,,(s), the principal factor tending 


mb 


' and 


to decrease the size of successive harmonics 1s é ; 
where 5’ is the axial distance from stator trailing edge 
to rotor leading edge (Fig. 1). For present-day axial 
compressors, b'/d = 0.2 so that exp (—2rb' d) = 
0.285. This indicates that the successive harmonics 
ultimately decrease by a factor of at least one-quarter 
This same factor is involved in the summations over 
min ,L™ and ,\/". Therefore a few terms should give 
sufficient accuracy for engineering purposes 


5Y8 
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A similar analysis of the harmonics of ,L” and ,/” 
shows that they behave like Rk *, which is a rather slow 


convergence. 


To be sure, almost all the functions involved in these 
sums are oscillatory in nature, so that these approxi- 
mate factors of convergence obtained by the use of 
asymptotic expressions may not hold exactly in any 
particular case. It is believed, however, that they give 
lower limits to the speed of convergence, and the actual 
convergence is even faster than they indicate. Most 
of the numerical examples worked out so far indicate 
that such is the case. 


Evaluation of H-Functions 


Before the calculations are started one more sum 
must be approximated, and that is the one in the H- 
functions, which comes from the Glauert expansions 
for the steady bound vortex distributions, ,y(x), Eqs. 
(21) and (33). The first coefficient, A,, gives a dis- 
tribution corresponding to angle of attack on an iso- 
lated flat-plate airfoil, while the second, Aj, gives a 
symmetrical elliptic distribution corresponding to cir- 
cular-arc camber. Since they are the only coefficients 
that affect the total steady circulation ,I’ and since 
they represent the effects of both camber and incidence, 
it is felt that the most important features of the flow 
will be retained if the series for ,y(x) is terminated 


after the second term—i.e., A, = 0,” > 2. 


With this simplification and the use of the recurrence 
relation for the Bessel functions, it is easy to show that 


A. 
A o 2 a A, 
ye" =A 


rom A, + A, : 


+ 51 So) on ] . _ 
H,, | rmoet/? — ah (J, + id) ¥ 


J, (59) 


This result can easily be applied to /7/,,° or II,’ by 
proper choice of the sign and the parameters, as indi- 
cated by Eqs. (42) and (31), respectively. In Eq. 
(59), the effects of incidence and of camber are sepa- 
rated. If there is zero incidence in the steady part of 
the flow, then Ao = 0 and the steady lift distribution is 
elliptical over the chord; in this case the first term 
vanishes, leaving //,, expressed as one term, independ- 
ent of A;. If, on the other hand, there is only incidence, 
then A; = 0 and the steady chordwise lift distribution 
is the flat-plate one; again, //,, is given by a single 
term, (J, + 7J;), independent of A,. 


The linearity of the lift and moment in //,, permits 
any combination of flat-plate and elliptical distribu- 
tions to be imposed by merely adding pure flat-plate 
and pure elliptical results in the proper proportions. 
Also, since /7,, appears as a factor in lift and moment, 
it is easy to convert elliptical results to flat-plate re- 
sults, and vice versa, merely by multiplying by the 
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proper (complex) factor. This fact proves useful in 
the numerical work. (This addition procedure will not, 
of course, give the correct wake energy, since the energy 
is nonlinear. ) 


The Bessel functions J,(re’”’) and J,(re’®) are tabu- 
lated’? for0O< r< 10 and0< @< 90° at intervals of 


5°. Since 0 < a,, a, < 90°, the argument (7/2) -— 
a, is in the table, and noting that J,(re~) = J,,(re®) 
the argument —[(72/2) — a,] can also be found there. 


(The bar denotes the complex conjugate.) 


Cases Calculated in This Paper 


There are five lengths that enter into the geometric 
configuration—namely, c,, c,, ds, d,;, and b’ (or } 
Therefore, in dimensionless form, four ratios made 
from these lengths must be specified for each configura- 
tion. The ratios oa,, o;, d,/d,, and b’’c, have been 
chosen as the most convenient. It was felt that the 
effects of variations of the pitch ratio d,/d, and the ratio 
of axial separation to rotor half-chord, b’/c,, would be 
the most interesting to investigate, and these two param- 
eters were varied in the numerical calculations, while 
o, and o, were held fixed at 1.0. 

Since only velocity ratios are involved in the theory, 
the velocity triangle is determined when a, and a, 
are specified. The only remaining parameter that 
enters the problem is the ratio of steady circulation, 
xi. 
steady part of the flow is such as to cause no net turning 
of the flow through the stage——i.e., the flow leaving the 
rotor is parallel to that entering the stator. Since the 
turning angle of a row of equal, equally spaced vor- 
tices, jl’, is proportional to ,I'/d, d being the pitch, it 


In the present paper it is assumed that the 


is easy to see that the condition of no net turning 
through the stage requires that 


s/n = —d,/a, (60) 


This relates the circulation ratio to one of the previ- 
ously chosen geometric parameters, the pitch ratio. 


The result is that, for no net turning through the 
stage, it is necessary to specify four geometric ratios, 
o,, os, d,/ds, and b’/c,, and two angles, a, and a,, in 
order to carry out numerical calculations. This is, of 
course, in addition to the condition of flat-plate or ellip- 
tical distribution discussed above. 


It is interesting to note that it is not necessary to dif- 
ferentiate between a turbine and a compressor in the 
calculations, Eq. (60) being valid for either. How- 
ever, the restriction of the whole theory to thin airfoils 
and solidities close to unity makes it more realistic for 
typical compressors than for typical turbines. 


For the numerical examples investigated in the pres- 
ent paper, values of the various parameters that are 
believed to be typical of present-day practice were 
Both rotor and stator solidities were fixed at 
A 50 per cent reaction stage was 


chosen. 
unity, ¢, = o, = 1. 
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AERODYNAMIC INTERFERENCE 
chosen so that a, = a,, and the angle a, = a, = 45° 
was used. As mentioned above, the effect of variation 
of the parameters d,/d, and b’/c, was investigated by 
varying each, holding the other fixed according to the 
following scheme: 


= 0.8, 1.0*, 1.25 
0.2, 0.4, 0.6 


b’/¢, = 0.2; d,/d, 
éJ/d,= 1.07 Ole = 


The first (fundamental) and second harmonics of the 
rotor and stator lift were calculated for these cases. 
The moments have not yet been calculated. 

The effects of changing from elliptical to flat-plate 
circulation distribution on the rotor blades were also 
investigated by performing these calculations for both 
A,’ = 0 (elliptical) and A,’ = 0 (flat plate). The con- 


* Although the case of equal pitch (d,/d, = 1) is commonly 
avoided in practice, this is dictated by considerations of mechani- 
cal vibration rather than by purely aerodynamic matters. No 
singular behavior occurs in the aerodynamics of the rigid cas- 
cades treated here at this value of d,/d,, and this value is a con- 


venient one to use. 
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Fic. 5. Absolute values of unsteady lift of rotor and 
stator blades for d,/d, = 1, as functions of axial spacing 
ratio b’/c,. Key: S denotes |;L*/,L*| and R denotes 
ol? / 91? Dash curves: elliptical steady-state chord 


wise load distribution on both rotor and stator blades 
Solid curves: elliptical steady-state chordwise load distri 
bution on stator blades; flat-plate distribution on rotor 
blades. m = 1: first (fundamental) harmonic. m = 2: 
second harmonic. 2c,/d, = 2¢,/d, = 1. ar = a, = 45 
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Fic. 6. Absolute values of unsteady lift of rotor and 
stator blades for b’/c, = 0.2, as functions of pitch ratio 
d,/d,. Key: Same as Fig. 5 


dition A,* = O (elliptical distribution on the stator 
blades) was maintained throughout. 

It might be mentioned that, since nearly all the num- 
bers involved in the calculations are complex and most 
of the manipulations involve multiplication, it was 
found most convenient to express the complex numbers 
in polar form in carrying out the numerical work. 

Numerical values of the Theodorsen function C(w) 
are available from many sources—for instance, from 
Luke and Dengler.'* The function S(w) may be cal- 
culated from its definition, Eq. (8), and the relations" 


= —(1/2)m[J,(w) — 1Y,(w)] 
—(1/2)rlJi(w) —7Vi(w)] 


K,, (tw) 
K (tw) = 


RESULTS AND CONCLUSION 


The results of the calculations for ,L* and for ,L’, 
which is the sum of ,L"! and ,L”, are shown in Figs. 5 
The magnitudes of the first two harmonics, 
land m = 2, of the lift ratio |,1/,L| are plotted 
against b’/c, and d,/d, for elliptical steady circulation 
distribution on the stator and for both flat-plate and 
elliptical distributions on the rotor. 


and 6. 
m = 
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' 
0 TABLE | 
Flat-Plate 
Elliptical Steady-State 
Steady-State Distribution 
Distribution on Rotor: 
on Stator Elliptical 
and Rotor Blades on Stator 
(1) C 
Stator wake: x 0.00017 0.0069 
05 (oCx")" 
Rotor wake: ( = - 0. 0000042 0. 000046 
L 7 
| In order to give some idea of the relative sizes of 
m= | .L’' and ,L”, the magnitudes of the first two harmonics 
=2 mes ° , | 9 | . ~* — 
ve (2) -——_—_# of | .L7'/,L"| and | ,L’*/,L"| are plotted in Figs. 7 and § 
Oo 4 f i siiateeitinllis ; ba a. 
7 or the same cases that were presentec ‘igs. 5 < 
0 02 0.4 O06 b 08 te t le same Case t at ere J ented in Figs. 5 and 6. 
m rhe first two harmonics of ,L”! are the m = 1, 2 terms 
Fic. 7. Absolute values of rotor unsteady lift com- of Eq. (50). The first two harmonics of ,L"? come from 
ponents for d,/d, = 1, as functions of axial spacing ratio ae saad PNG = a se ek ae ee 
b’/c,. Key: Same as Fig. 5 except (1) denotes |;L7'/L’ the k = +1, +2 terms of the sum over & in Eq. (90), 


and (2) denotes |,L7?/,L’ 


It is seen that the nonsteady part of the lift increases 
sharply with decreasing b’/c,, as might be expected 
from the exponential nature of the factor in which this 
parameter appears. The change of lift with variation 
of d,/d, is also appreciable. A flat-plate distribution 
on the rotor blades produces much larger effects on the 
stator than does elliptical distribution. This is to be 
expected, since the centroid of the rotor steady bound- 
vortex distribution is a quarter chord closer to the stator 
in the flat-plate case than in the elliptical case. The 
change of rotor steady circulation distribution also 
affects ,L’ somewhat, because ,L’”? is altered, although 
,L’' remains unchanged. 

A general conclusion to be drawn from the results 
shown in Figs. 5 and 6 is that the nonsteady part of the 
lift may be as large as 18 per cent of the steady lift 
and therefore may be of practical importance, espe- 
cially with reference to blade vibration. 
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| 
a 
O (2)8—~ —_ __ 4° = 
0.8 1.0 125 d, 
ds 
Fic. 8. Absolute values of rotor unsteady lift com 
ponents for b’/c, = 0.2, as functions of pitch ratio d,/d, 
Key: Same as Fig. 7. 


and in the present calculations two terms in the sum 
over m in that equation were used. It must be re- 
membered that | el” is not the algebraic sum of 
|,L’') and | ,L™), because these are the real parts of 
complex quantities. 


It is seen that ,L”' and ,L” 
especially for flat-plate distribution on the rotor blades, 


are comparable in size, 


which increases ,L’ considerably while, of course, 
leaving ,L’' unchanged. The slower decrease in size 
of the successive harmonics of ,L”, as compared with 


:L", as noted above, is also apparent. 


Wake Energy 


As has already been mentioned, the equations for 
the rate of transfer of energy to the vortex wakes have 
been worked out. These have been applied to two of 
the numerical cases treated above—viz., d,/d, = 1, 
b’/c, = 0.2, elliptical distribution on stator blades, and 
both (a) elliptical and (b) flat-plate distribution on 
The results are summarized in Table 1, 
where Cy* and Cz’ are the respective coefficients defined 


rotor blades. 


by 
Ce = E/pV°*c 
FE being the (average) rate of energy transfer per unit 
length of wake per unit width. 
These results can easily be converted into the ratios 
of the energy produced per unit time in the wake to the 
steady power necessary to turn the rotor, which is 


oL’U cos a,. If the ratios are denoted by P* and P’, 
then 
Cr* fi ¥,. he, Cr’ oC’ V, 
ih Fe be ee ; Fras — 
(07°)? \4k’ | V> U 68: a, (,CxL’)? U cos a; 


For the cases calculated here, 


U cos a, 


—d,/d, = —1 


so that 
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. § 


pe = [Ce*/(.Cx*)*)C.*, PY = [Ce"/(Cx")*),Cx’ 


and the numbers given above are precisely the wake- 
energy power to rotor-power ratios, except for a factor 
.C.. The power in the vortex wakes is thus seen to 
be small compared with the power necessary to drive 
the rotor. 

A striking feature of the wake-energy results is the 
fact that the energy in the stator wake is about 100 
times as large as that in the rotor wake. This implies 
that the nonsteady circulation on the stator blades is 
about ten times as large as that on the rotor blades, 
since the wake energy is proportional to | ,I°}*. Al- 
though this seems surprising at first, it has been verified 
that, because of the difference in chordwise velocity 
distribution on the stator (upstream) and rotor (down- 
stream) blades, as indicated by Eqs. (36) and (40), it is 
indeed true that the nonsteady circulation on the 
stator is considerably larger than that on the rotor, 
although the nonsteady lift is approximately the same 
on both. 
ie., that the circulation induced on a given row by the 


This can be expected to be a general result 


relative motion of an adjacent downstream row will be 
considerably larger than that induced by the adjacent 


upstream row. 


Verification of Basic Assumption 


The assumption that | I} < |,[f|, on which the 
present approximate treatment is based, can also be 
checked by numerical calculation in order to verify 
that the theory given in this paper is a consistent one. 
The expression for ,I’ for each blade can be found, as 
mentioned above, from Eqs. (47) and (52), and the 


general formula for any sinusoidal motion!'! 


T@) = THe “Sle) 
These give 
r/o? = (0/0) > Fae "S(wm) 
m 
py = Ca 3") Re Fg " Slaome) 


m 
and a similar expression for ,I’’? can be found. An upper 
is obtained by taking the absolute 
For the 


limit for | I(T 
value of each term in the sum on the right. 
case d, d, = 1, b' c, = 0.2, these upper limits, based 
on two terms of the series, are shown in Table 2. It 
appears that for the other cases these numbers would 
be much the same, so that the basic assumption is well 


verified by the results of the theory. 


CONCLUSIONS 


An approximate calculation has been carried out, 
based on the assumptions of incompressible inviscid 
two-dimensional flow, of the unsteady aerodynamic 
interference between two cascades of thin airfoils in 

The basic approximation, which is 
the results obtained for 


relative motion. 


confirmed by numerical 
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TABLE 2 
Elliptical Flat-Plate 


Steady-State Steady-State 


Distribution Distribution 
on Stator on Rotor; 
and Rotor Blades Elliptical on Stator 
0.02 0.13 
0.0027 0.030 


typical present-day compressor cascades, is that the 
unsteady part of the circulation about any blade is 
small compared with the steady part. Formulas for 
the unsteady components of lift and pitching moment 
of rotor and stator blades have been obtained in the 
form of successive harmonics. 

The limited numerical results presented here indicate 
that: 

(1) The fundamental harmonic component of vari- 
able lift varies in magnitude from | to 18 per cent of 
the steady value, depending especially on the axial 
distance between rows and on the chordwise distribu- 
tion of steady circulation of the blades of the other row. 

(2) The second harmonic component of variable lift 
js considerably smaller than the fundamental. 

(3) The contribution to of the 
wakes shed by stator blades, through which the rotor 
blades pass, is comparable in size to the effect induced 


rotor lift vortex 


by the stator blades themselves. 

(4) The fluctuating part of the circulation, expressed 
as a ratio to the steady part, is several times as large 
for stator (upstream) blades as for rotor (downstream) 
blades, even though the corresponding ratio of lifts 
is about equal. 

(5) The rate of transfer of energy into the vortex- 
wake pattern which results from the variable circula- 
tion is small compared with the power required to turn 
It is about 100 times as large for the stator 
the (downstream) 


the rotor. 


(upstream) blades as_ for rotor 


blades. 
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Stall-Flutter in Cascades’ 


F. SISTOt 
Wright Aeronautical Division, Curtiss-Wright Corporation 


ABSTRACT 


Assuming the complex of oscillatory aerodynamic reactions to 
be small compared with blade inertia and elastic forces, some re- 
sults from nonlinear mechanics are utilized to develop a con- 
venient method of studying stall-flutter oscillations. This repre 
sents a new approach to the problem. 

Data are presented from an exploratory experimental investi 
gation of a ‘‘two-dimensional”’ cascade in which one blade and 
later three blades are flexibly mounted. The variation of flutter 
amplitude with stream velocity and mean incidence is obtained 
for angles of attack near the static stall. Experimental deter 
minations are given for three different gap/chord ratios. The 
motion is predominantly torsional. 

The results confirm certain deductions from the nonlinear 
theory concerning transient behavior and equilibrium flutter 
amplitudes or “‘limit-cycles.’’ Other qualitative results are re- 
ported which confirm the necessity of adopting this general ap- 
proach in order to explain the important features of stall-flutter 
phenomena. 


INTRODUCTION 


_— FATIGUE AND EVENTUAL FAILURE of the aero- 
dynamic blading in axial-flow compressors is often 
attributable to stall-flutter oscillations. At the large 
mean angles of attack where this type of vibration is 
encountered, the flow about the blade is separated 
during the whole, or at least a part, of the cycle of os- 
cillation. Because of the difficulty involved in the aero- 
dynamics of unsteady stalled flow most of the previous 
work has been wholly experimental! or has consisted of 
empirical and semiempirical* * modifications of classi- 
cal flutter theory. 


In the sequel, the nonlinearities in the oscillatory 
aerodynamic reactions are considered directly by stand- 
ard methods of nonlinear mechanics. 


Finally, the results ef an exploratory experimental 
investigation are presented in corroboration of the 
gross features of compressor blade stall-flutter phenom- 
ena. Eventually, this general method, in conjunc- 
tion with suitable experimental work, should lead to 
reliable quantitative information. 


PRELIMINARY CONSIDERATIONS 


Certain simplifications that appear to be allowable in 
the turboblade problem consist in assuming that (1) 


Presented at the Flight Propulsion Session, Twenty-First 
Annual Meeting, IAS, New York, January 26-29, 1953. 

* Portions of this paper are from an Sc.D. Thesis investigation 
conducted at the Gas Turbine Laboratory, Massachusetts Insti- 
tute of Technology. Permission to include these portions is 
gratefully acknowledged. 

t Assistant Project Engineer. 


the ratio of air density to blade density is small and 
consequently, (2) the unsteady aerodynamic forces 
may be taken to be small compared with the elastic 
and inertia forces of the blade, (3) the frequency of 
flutter will be close to the natural frequency of the blade 
in the predominant normal mode, and (4) the equilib- 
rium motion of a sustained flutter condition will be 
“nearly sinusoidal” in the principal mode. 

For simplicity we consider only a single degree of 
freedom system acted upon by two-dimensional aero- 
dynamic forces. 

Under the previous assumptions, the differential 
equation of motion (either pure torsional or pure trans- 
lational degree of freedom) is 


¥ + win 'X + wx = Fix, X) (1) 


where w is the natural frequency of the system, 6 is the 
logarithmic decrement im vacuo, and F is the ratio of 
the small aerodynamic forces to the system inertia 
(mass or moment of inertia). 


Presupposing a ‘‘nearly harmonic’ solution of the 


form 

xX = Xp COS wi (2) 
where x is a slowly varying quantity, it is possible to 
arrive at an approximate expression for the rate of 
change of amplitude 


2a 

f;> = —(2 — 1 [ F(x» cos —- 34 

Xo 2TX yw) (Xp COS T, Xow Sin T) + 
0 


w*bxor' sin 7] sin rdr_ (3) 


by an energy consideration® or by one of the standard 
methods of nonlinear mechanics.® ? 

The equilibrium flutter amplitudes consequently are 
given by the condition 


Qn 
f [F(x9 cos r, —Xow sin r) + w*6xpr~! sin 7] sin 7 
0 


dr =0 (4) 


where only the real positive values of xo satisfying this 
equation are significant. Whether these values of 2 
represent stable{ or unstable motions is determined by 
the criterion 


OXo(Xo) 


0 stable 
Oxo _ 
(2) 
OXo( Xo) 
0 unstable 
Oxo 


t Stability according to Liapounoff.* 
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STALL-FLUTTER 
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Fic. 1 Mean dynamic normal force coefficient for the illustra 
tive example in pure bending stall-flutter 
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Fic. 2. Flutter amplitude in pure bending as a function of angle 


of attack and reduced velocity 


Hence, the stable equilibrium flutter amplitudes may 
be obtained to a first approximation within the accu 
racy of the aerodynamic reaction F(x, x). 
Experimentally,” § F is known to be a more or less 
complicated function that exhibits an aerodynamic lag 
or “hysteresis” effect. For the time being it is con 
venient to ignore this hysteresis and to suppose that a 
single-valued mean dynamic characteristic exists. 
This characteristic is either aerodynamic moment or 
normal force as a function of instantaneous angle of 
attack. The angle of attack, however, is a function of 
blade velocity transverse to the air stream, so that a 
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certain amount of “‘hysteresis’’ reappears in the aero- 


dynamic reactions. 


PURE BENDING STALL-FLUTTER 


The previous concepts are best illustrated by some 
simple examples. A pure bending stall-flutter will be 
considered first. In this case put 


“=f; Xo = Re, etc. (6) 


and 


Fih, h) = F(h) = 2gbf(a)/m ( 


where a is the instantaneous angle of attack and is 


given by 


a= aj; t (h V (S) 


In these expressions } is the blade semichord, |” is the 
mean stream velocity at upstream infinity, g is the 
dynamic pressure based on this velocity, a; is the mean 
value of a, and m is the blade mass per unit span. 
Assume, for the purposes of this example, that / may 
be expressed as a polynomial 
Ila) > a,a’ (9) 
where the a,, could be adjusted to fit the mean dynamic 
characteristic if this characteristic were known. Re 
taining the first three terms and applying the previous 
development, the expression for the stable equilibrium 
flutter amplitudes (limit-cycles) in pure bending is 


S{f’’’(a,)]-! X 


, ( m \(“") Pal 
—f'(a;) + ~ 1(6) (10) 
mpb*/\ 1 wh 


Only positive values of the 


(hy /b)* 


where p is the air density. 
right-hand side of Eq. (10) are significant. 

Taking a; = —5, a2 = 16, a; — 16, for instance, 
defines the mean dynamic normal force coefficient illus 
trated in Fig. 1. The stable equilibrium flutter ampli 
tudes are shown in Fig. 2 in nondimensional form for 
various values of the reduced velocity |’ (wb), and with 

(m_/ mpb*)(6) 15 

Further, if V (wh) is cross-plotted against a, for con 
stant values of /y)/b, the conventional ‘‘critical”’ flutter 
velocity diagram illustrated in Fig. 5 is obtained. 
The term ‘‘critical’’ is somewhat misleading, since the 
flutter condition may or may not develop explosively 
as a governing parameter is slowly varied. In the 
present instance of so-called soft flutter, the amplitude 
grows smoothly from zero as the flutter boundary 1s 
traversed by means of slowly varying the air-stream 
velocity or angle of attack. In hard flutter, on the 
other hand, the amplitude may suddenly jump to a 


finite value. 








PURE TORSIONAL STALL-FLUTTER 


Pure torsional stall-flutter, by a similar analysis, 
proves to be somewhat more complicated, since the in- 
stantaneous angle of attack due to rotational velocity 
is a strong linear function of position along the blade 
chord. By way of illustrating certain other features, 
however, put (with x = 6, x» = 4, etc.) 


a= a;+0+ (nb6/V) (19) 


where 7 is so chosen as to give the instantaneous angle 
of attack at the three-quarter chord point. (This con- 


(V wb)? i g 


? = 12 : ~ 
n? + (V/wb)? | g 
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Fic. 3. Conventional flutter boundary of “‘critical’’ reduced 
velocity versus angle of attack for the pure bending example 
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vention is taken over directly from classical flutter 
theory’ as a crude approximat:on in stall-flutter. ) 
Then, taking 
F(0, 6) = 4qb?g(a)/I (19) 


and assuming a fifth degree polynomial for g 
5 
g(a) = > ba” (13 
n= 1 


there results, after a similar reduction, the expression 
for the limit-cycles in pure torsional flutter 
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where / is blade moment of inertia about the elastic 





axis per unit span. 
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The results of calculations for an illustrative case in 


which 
n = 0.85 
b = h = 0 
b, = —3/5, b:=11, b = —30 
(1 wrpb*)(6) = 3 


are shown in Figs. 4, 5, and 6. Data similar to Fig. 5 
were obtained in the earliest stall-flutter experiments 
Fig. 6 is included to indicate the possibility 
The dotted lines are unstable limit- 
If, for a given 


by Studer. ! 
of hard flutter. 
cycles according to the inequalities (5). 
angle of attack and air-stream velocity, the operational 
point lies in this region of the figure, the blade will re 
main at rest unless impulsively started with an ampli- 
tude greater than the dotted portion of the curve. 
This situation is illustrated in Fig. 7, which indicates 
that it is possible to traverse a “‘hysteresis’ loop of a 
different kind by cyclically changing the velocity. 


EFFECT OF CASCADING 


Embedding a single blade in an infinite cascade of 
similar blades gives rise to two effects: 

(1) The shape of the mean dynamic force characteris- 
tic is changed because of the proximity of other blades. 

(2) The movement of these adjacent blades further 
modifies the “‘instantaneous’’ characteristic, in addition 


‘“Hysteresis’’ loop in amplitude versus velocity obtainable in hard flutter (a@ 


= ().175 


to giving rise to possible transient effects that travel 
along the cascade (i.e., progress from blade to blade). 

In the first case F becomes a function of blade spac 
ing and stagger in addition to x and x. In the second 
case the displacement and velocity of all the additional 


blades must enter into F. 


EXPERIMENTAL RESULTS 


Considering only the first type of cascading effect, 
experiments’ were conducted in a low-speed cascade 
tunnel in which all blades but the central one were 
rigidly attached to the tunnel walls. The central steel 
blade was flexibly mounted on springs to allow inde 
pendent motion in torsion and bending, the blade itself 


TABLE | 
Experimental Blade Constants 


a = —0.35 

b = 0.112 ft 

Xx, = 0.192 

Wp 170 rad. per sec 

w, = 48.6 rad. per sec 

I, = 1.74 (10~‘) slug-ft.?/ft. span 
M, = 0.065 slug/ft. span 

M, = 0.050 slug/ft. span 


Profile NACA 6 1/2 409 
Blade span 5.5 inches 
log dec 6, < 0.05 


log dec 6, < 0.01 
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being completely rigid by comparison. The constants 
of the system, in the Theodorsen notation,’ are listed 
Structural and aerodynamic conditions 


Data were ob- 


in Table 1. 
were substantially ‘‘two-dimensional.”’ 
tained by means of electric resistance strain gage pick- 
ups attached to the springs. The signals from the 
gage bridges were amplified and recorded on conven- 
tional oscillographic equipment. 

The results of this investigation are illustrated in 
Figs. 8, 9, and 10 for three different gap chord ratios 
at zero stagger. While exact correspondence with the 
illustrative example is not to be expected, general 
curve shapes and broad trends are substantially as 
predicted. Two influences that cause deviations from 
predicted variation are: 
motion in some regions of the curves (datum points 
with ‘‘tails’’) and (2) the degree of approximation in- 
(Time- 


(1) the coexistence of bending 


volved in using a fifth-degree characteristic. 
average values of moment in pure pitch determined 
experimentally in reference 2 indicate that the use of a 
seventh- or ninth-degree polynomial would greatly 
improve the representation of a mean dynamic moment 
coefficient.) The effects of bending motion on tor 
sional amplitude may be explained satisfactorily if it is 
assumed that the aerodynamic reactions due to this 
bending motion do positive work on the blade in the 
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regions where the torsional amplitude is larger than 
the theoretical type of variation would predict. 

The self-starting (e.g., point S in Fig. 7) and quench. 
ing (e.g., point Q in Fig. 7) limits could not be deter. 
mined accurately. A certain amount of residual coy- 
lomb damping in the torsional spring supports made it 
impossible to obtain a completely soft pure torsional 
flutter. 
the starting transients indicated that in some cases a 


Furthermore, a close examination of some of 


soft bending flutter was sufficient to start the torsional 
mode into vibration. The resulting steady amplitude 
torsional flutter completely obscured the important 
“triggering” effect of the bending instability. 

This characteristic of nonlinear systems is referred to 
as asynchronous excitation and quenching,’ the term 
“asynchronous” denoting that there is no relation be- 
tween the frequencies of the two oscillations (bending 
and torsion). 

The second type of cascading effect is illustrated by 
two series of experimental determinations conducted 
with three spring-mounted blades. The results of one 
such series is compared in Fig. 11 with the result for a 
single spring-mounted blade. In addition to showing 
a consistent reduction in amplitude at this particular 
mean angle of attack, the records for the three fluttering 


blades indicated either an amplitude or a frequency 
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amplitudes for gap/chord = 1.5 
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ing blade 


The dotted line represents mean values 
The frequen- 


modulation. 
of the amplitude modulated wave form. 
cies in the upper portion of the diagram refer to the 
configuration with three flexibly mounted blades, «, 
being the frequency of the central or reference blade 
and w, the frequency of the adjacent blade on the convex 
side of the reference. The point where the two fre- 
quencies separate corresponds to the onset of the ampli- 
tude modulation mentioned before. Below this point 
the frequencies recorded are long-time averages (i.e., 
The 
phenomenon exhibited by the two separate frequencies 


over several cycles of frequency modulation). 


w, and ws is known as synchronization or entrainment of 
frequency and is predicted by nonlinear theory.’ 

Fig. 11 also affords a particularly accurate determina- 
tion of amplitude variation with velocity. The similar- 
ity with predicted curve shape is apparent. 


EXTENSION OF THE METHOD 


The experimental investigation indicates a serious 
handicap in considering a single degree of freedom sys 
aerodynamic and mass coupling 
The analysis 


tem theoretically; 
phenomena are automatically excluded. 
for a two degree of freedom system, including first- 
order coupling effects, may be performed by applying 
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the Method of Linearization.” When 
sufficient experimental data become available to war. 
rant their use, the general method may be extended jn 


Equivalent 


this direction. 

Another approach to aerodynamic coupling in cas. 
cades presupposes each blade to have only one degree 
of freedom. In this case each blade is an “‘identica]” 
element, and the analogy may be drawn to infinite lat. 
tices of identical coup!ed circuits in network theory. 
The problem may then be attacked by the successive 
application of: the Method of Equivalent Lineariza- 
tion, the Method of Complex Amplitudes, and the 
theory o¢ linear difference equations. 

Finally, when the stage is reached where it is desired 
to determine the mean dynamic aerodynamic character 
istics (i.e., f and g or their counterparts) from experi- 
representations for these functions. On the contrary, 
it is desirable to expand each in Taylor series about 











a;, terminate the series at some convenient point, and 
pursue the analysis as before. The experimental data 
may then be used to solve for the first derivative, and 
the form of f (or g) may be obtained within an additive 
constant by graphical integration. 
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An Exact Solution for the Bucklin 


Flat Sandwich Panel 


» Load of 


Ww 


s with Loaded Edges 


Clamped™ 


S. V. NARDO? 
Polytechnic Institute of Brooklyn 


SUMMARY 


This paper presents an exact solution for the buckling load of 
flat rectangular plates whose loaded edges are clamped and un 
loaded edges simply supported. The theory is also shown to be 
applicable for the case in which unloaded edges are clamped and 


loaded edges simply supported. 


The solution is based on small deflection theory, and it is as 
sumed that the stresses are always below the proportional limit 
The material for both the faces and the core is considered iso- 
tropic. Three equations of equilibrium, derived by means of the 
principle of virtual displacements, form the starting point of the 
theory. The equilibrium equations are partial linear differential 
equations with constant coefficients 
and the third is of the fourth order. 
equilibrium and 16 boundary conditions are satisfied identically 


two are of the second order, 
Both the equations of 


Numerical calculations have been made for the case of loaded 
edges clamped. A series of design charts is presented which 
gives the buckling load as a function of the plate aspect ratio and 
the ratio of core to plate thickness. Each curve is valid for a 
particular value of a parameter representing the ratio of shear 
rigidity to bending rigidity. Values of this latter parameter, 
which are four in number, are chosen to span the range encoun- 


tered in most practical applications. 


A check with experiments run at the Forest Products Labora- 
tory indicates that the theory is in reasonable agreement with 


the experimental results 


SYMBOLS 


2a,2b = edge lengths of sandwich panel, in 

( = core thickness, in 

C = buckling coefficient, o,,/ Fo-,, 

E = Young’s modulus, lbs. per sq.in 

F = 1+ 3[(c/t) + 1]? 

G = shear modulus, lbs. per sq.in 

h = total thickness of sandwich plate, in. 


Presented at the Structures Session, Twenty-First Annual 
Meeting, IAS, New York, January 26-29, 1953. 

* This paper is based upon a thesis submitted to the Poly- 
technic Institute of Brooklyn in partial fulfillment of the require- 
ments for the degree of Doctor of Philosophy (in Applied Me- 
chanics). 
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a contract awarded to the Polytechnic Institute of Brooklyn by 
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technic Institute of Brooklyn. The author is grateful to these 
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Salerno, and to Burton Erickson, Ivan Rattinger, and 


edges his indebtedness to his thesis advisers, Profs. N 
and V. L 
Dominick Torrillo for their assistance in carrying out part of the 
numerical calculations 
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n = number of half waves in sandwich plate in x or 


direction 


' = critical value of edge load, lbs. per in 

R = relative rigidity parameter, 12(1 v*)Ga?/ Et 

t = face thickness, in 

u,v,’ = sandwich plate displacements, in 

x, ¥, 2 = coordinate system for sandwich plate 

~,n, ¢ = nondimensional coordinate system for sandwich 
plate 

r = critical load coefficient 

v = Poisson’s ratio 

a = buckling stress of sandwich plate, Ibs. per sq.in 

o = buckling stress of two independent faces, lbs. per 


sq.in., 7? Et?/3(1 v?) (2a)? 


INTRODUCTION 


i be TERM “SANDWICH CONSTRUCTION’ generally 
refers to a flat or curved load-carrying surface, 
composed of two thin high-strength faces on either side 
of a comparatively thick low-weight low-strength core. 
In aircraft practice, aluminum alloy has been widely 
used for face material, but also suitable in certain ap- 
plications are Papreg or Fiberglas. Balsawood, 
oriented in such a manner that the fibers are perpen- 
dicular to the faces, has been popular as a core material. 
Cores have also been made of expanded synthetic ma- 
terials and built-up square or honeycomb grids of the 
same material as the faces. Considerable research 
has been conducted to develop a more satisfactory core 
material—i.e., a core material with relatively high 
strength compared to its weight. 


The modulus of elasticity of the core in the plane of 
the faces may be as low as one-thousandth of that of the 
faces. Asa consequence, the normal stresses developed 
in the core can be neglected in comparison with those 
developed in the face material, even though the ratio 
of core to face thickness may lie between 10 and 100. 
The core, however, is counted upon to transmit shearing 
forces. Because of the core’s low shear modulus, con 
siderable shearing deformations take place which can- 
not be disregarded as they are in the case of thin plate 
theory. 

Toward the end of, and immediately following, World 
War II, considerable work was done in the field of 
sandwich construction. Many airplanes have been 
designed utilizing sandwich-type construction for the 


major structural components. For this reason the 
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strength of sandwich construction was also investi- 
gated, and various theories were proposed for the 
analysis of sandwich plates. 

The first among these was one developed by Leggett 
and Hopkins! in 1942, which was followed by the theory 
of Van der Neut? in 1943. 
derived differential equations, while March and Smith* 
used the strain-energy approach in 1945. A simplified 
procedure was suggested by Bijlaard‘ in 1946, and in 
1948 Reissner,’ as well as Libove and Batdorf® pre- 
These were 


These two investigators 


sented new sets of differential equations. 
followed by Hoff’s’ theory, which is the basis of the 
present investigation. 

Most of the work cited dealt with sandwich plates 
having simple supports along all four edges, because 
this condition leads to particularly simple solutions. 
The condition of simple supports along the unloaded 
edges and rigid clamping along the loaded edges has 
been given much less attention. 

This boundary condition, as well as the one in which 
the unloaded edges are clamped and the loaded edges 
simply supported, was first solved for ordinary (not 
sandwich-type) thin plates by Schleicher in 1931.* ° 
Schleicher assumed that the deflected shape of the plate 
could be expressed as the product of an unknown func- 
tion in the loaded direction and a trigonometric func- 
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tion in the unloaded direction. Substitution of this 
function into the partial differential equation of equilib. 
rium reduced it to an ordinary differential equation jn 
the unknown function. 
conditions led to a set of two linear homogeneous equa- 


Enforcement of the boundary 


tions whose denominator determinant was required to 
Since the 
elements of the determinant contained the edge load, 


vanish if the solution was to be nontrivial. 


this gave a condition yielding the critical load. 

An exact solution to the buckling of flat plywood 
plates was published in 1944 by R.C. T. Smith.'° 
it was assumed that in the expression for the de- 
The solution for 


Again 


flected shape the variables separate. 
this plywood plate-buckling problem was given for a 
number of boundary conditions exactly. In the case 
of all four edges clamped, however, a Rayleigh-Ritz 
procedure was utilized. 

In 1945, March and Smith* gave the strain-energy 
solution to the buckling of flat sandwich plates, loaded 
edgewise in one direction, under various edge condi- 
The plates were considered orthotropic. The 
straightforward energy 


tions. 


procedure involved was a 
method, the work done by the edge load being equated 
to the strain energy of bending of the panel. Curves 
were presented for all cases and for various degrees of 
orthotropy. 

PROBLEM AND Basic ASSUMPTIONS 


STATEMENT OF 


The problem investigated in the present paper is the 
exact solution to the buckling load of a flat sandwich 
panel, loaded edgewise in one direction, with loaded 
edges clamped and unloaded edges simply supported. 
The method used in the solution is also satisfactory 
when the boundary conditions are reversed— i.e., the 
loaded edges are simply supported and the unloaded 
edges clamped. 

Hoff’s differential equations were used in the analysis 
because they represent the behavior of sandwich plates 
more accurately than the theories in which the sand- 
wich plate is replaced by a homogeneous thin plate 
weak in shear. At the same time, they are not so dif- 
ficult to solve as the theories in which the sandwich 
plate is treated rigorously as the assembly of three 
interconnected plates. These equations have already 
been solved for the cases where all the four edges are 
simply supported’ and when the loaded edges are 
simply supported and the unloaded edges rigidly 
clamped."! 

This theory assumes that the deflections are small 
and the stresses are less than the proportional limit. 
The edge load is applied only to the faces and not to 
the core. 

In the derivation of the equations of equilibrium, 
the following strain-energy quantities were considered: 
(1) strain energy due to extension and shear of the 
faces in their own plane; (2) strain energy of bending 


of the faces; and (3) shear strain energy of the core. 
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Note that any strain energy due to transverse shearing sandwich plate coincide with and » = +1. A sketch 
of the faces and to extensions in the core parallel or per- of the loading, coordinate system, and displacements 1s 
pendicular to the faces was considered negligible. shown in Fig. 1. Note that u represents a displacement 


The displacements are described by means of three _ in the positive £ direction in the upper face and an equal 
variables u, v, and w as functions of the independent but opposite displacement in the lower face. The v 


variables, £ and n, the coordinates in the plane of the displacements are similarly defined, except that they 


plate. These latter coordinates are nondimensional obtain in the direction. The plane of the plate be- 
and are chosen in such a manner that the edges of the comes curved only during w displacements. 


THEORETICAL SOLUTION 


The equations of equilibrium as derived in reference 7, after setting the transverse load and the edge load acting 


in the x(£) direction equal to zero, may be written as follows: 
— By[2(02u Ox?) + (1 — v) (O2u Oy?) + (1 + v) (0°v/Oxdy)] + 2 Geu + Ge(e + t) (Ow/Ox) = 0 (1) 
— By[2(0°v/ Oy?) + (1 — v) (O°v Ox?) + (1 + v) (O°u/Oxdy)] + 2 Gev + Gelc + t) (Ow/dy) = 0 (2) 


B,[(O'w Ox!) + 2(0'w/ Ox*dy") + (O'w Oy')] — [2Ge/(c + t)] [(Ou/Ox) + (Ov/Oy)] — Ge[(O?w/Ox*) + 
(0°w/Oy?)] + P.-(0?w/Oy?) = 0 (38) 


where By = Et(c + t)?/2(1 — v®) and B, = Ef® 6(1 — v?). Making a transformation of coordinates — = x/a and 


n= y/6, Eqs. (1), (2), and (3) become 


— By} (2/a®) (0°u/dE?) + [(1 — v), b?] (07% /On?) + [C1 + v)/ab] (O°v 0£0n) {| + 
2Gcu + Ge(c + t) (Ow/ OE) = O (4) 


—By{ (2b?) (0°%v/On?2) + [C1 — v) a?] (0° dE?) + [C1 + v)/ab] (O%u/dEOn) fF + 
2Gcv + Ge(c + t) (Ow/On) = O (5) 


B,[(1/a*) (O'w/OE*) + (2/a2b?) (Otw OE7On*) + (1/54) (O*w/On*)] — Ge{2/(e + t)] [(1/a) (Ou/ OE) + 
(1/b) (Ov On) | — Ge[(1 a?) (O°w/OE*) + (1/5?) (O?w/On*?)] + (P.-/b?) (0°w/On*?) = 0 (6) 


A solution for the displacements is assumed in the form: 


u = a,e“"" sin (nw/2)E (7) 
vy = b,e*"" cos (nmr/2) & (8) 
w= c,e"" cos (nw/2) & (9) 


If the foregoing solution is substituted into Eqs. (4), (5) and (6), the equations of equilibrium are transformed into 


three simultaneous algebraic equations, 


—A,[—2(n7,2)*a, + (Ra,)? (1 — va, — (Ra,) (1 + v) (ue 2)b,] + a, — Az (nw/2)c, = O (10) 
—A,|2(ka,)°b, — (nmw/2)? (1 — vb, + (ka,) (1 + v) (nm/2)a,] + b, + Ae (ka,)c, = 0 (11) 
A;|(nm 2)? — (ka,)?]2c, — (1 As) [(nm/2)a, + (Ra) b,] + [(nw/2)? — (Ra,)*] On + MRa,)? 6, = O (12) 
where 
A, = 3[(c/t) + 1]? Ri(c/t), A = P.,/Ge 
As = (t/2a)/[(c/t) + 1], R = 12(1 — v*?)Ga?/ Et? 
Ag = 2/Ric/t), k = a/b 


The above equations represent three linear homogeneous equations whose denominator determinant, formed 
by the coefficients of the a,, 6,, and c, terms, must vanish. Enforcement of this condition leads to the following 
quartic algebraic equation: 


i 


[f, + 1/A,(1 — v)] | —2AjAsf,® — [24101 — A) + As] fi? — A[2Ai(mw/2)? — 1] f, — Maw 27} =O (13 


where 


{, = (nw/2)? — (ka,)? (Isa) 
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One of the roots of the above equation is immediately obtained, 
fin = —1/Ai(1 — ») 


Note that this root is dependent only on the physical and geometrical properties of the sandwich plate. The re- 
maining three roots, /2,, f,, and fs,, are obtained by solving the cubic expression in the brackets in Eq. (13). Once 
the four roots of Eq. (13) are determined, the eight corresponding values of a follow from the substitution pre- 
viously made [Eq. (13a) ]. 

The eight roots in @ are in the form +a;,, +a»,, +a3,, and +ay,. If we substitute these values into the dis- 
placement functions and stipulate that the displacements in the loaded direction be symmetrical, the expressions 
for u,v, and w can be written as follows: 


u = [Cyd Cosh ayn + CrpA2_ COSH a2 77 + C343, CoSh a3,n + Cy,04, Cosh as,n] sin (nw/2)é (14) 
v = [C,,)1, sinh ain + Cobo, sinh ae,7 + C3,b3, sinh a3,n + Cy,b4, sinh as,n] cos (nw/2)E (15) 
w= [Cin cosh Qin = i Crp cosh Q2nN sa Can cosh Q3n7 = Ce. cosh a4, | cos (m7, 2)é (16) 


In writing down the foregoing displacement functions, Eqs. (7), (8), and (9) were redefined by arbitrarily as- 
signing a value of unity toc,,. 

Thus, the a,,’s and 6,'sin Eqs. (14) and (15) correspond to c, = 1 in Eqs. (7), (8), and (9). By substituting the 
above displacement functions into the equations of equilibrium, it can be shown that 


a, = Ax(nmw/2)/2Aif, + 1 (17) 

b, = —A2x(ka,)/2Aify + 1 (18) 

BOUNDARY CONDITIONS maining four boundary conditions yield the following 
equations: 


For the present problem, the sandwich plate is 
loaded edgewise in one direction and has clamped _ (C,, cosh a, + C»,, cosh a», + C3, cosh a3, + 
loaded edges and unloaded edges that are simply sup- C,, cosh a4, = 0 (27) 
ported. The following boundary conditions are re- 


2 es Ci, 0%, sinh a,, + C.,a2, sinh as, + C3,a3, sinh a 
quired to be satisfied: Inn In 2n 2p 2n 3nQ3n an T 


C;,04, sinh ay, = 0 (28) 

w= 0 (até = +1) (19) 
w=0 #£(atn = +1) = (20) 
ow /On = 0 te= &1) (9) 
O°w/dt? = 0 (at & = +1) (22) Cibin Sinh ay, + Conde, sinh oo, + Cagbs, sinh as, + 


cosh A2n + C3,43y cosh a3, aa 
C4,d4, cosh ay, = 0 (29) 


CinQ1, cosh aj, + C2,d2, 


U 


u=0 (atyn = +1) (23) C4,,b4, sinh as, = 0 (30) 
v=0 (atn = +1) = (24) 


a The boundary conditions therefore yield the fore- 
v=0 (at & = +1) (25) . : 


Ou/OE + v(Ov/On) = O (até= +1) (26) 


going four linear simultaneous homogeneous equations 
in four unknown arbitrary constants. The only possi- 
Boundary conditions (19), (22), (25), and (26) at the | bility for a nontrivial solution occurs if the denominator 
unloaded edges ( = +1) are satisfied identically by determinant, formed by the coefficient of the C’s, is 
the trigonometric functions. Enforcement of the re- equal to zero. 


cosh a), cosh ay, cosh az, cosh ay, 
a, sinh ay, ay, sinh ay, a3, sinh a3, a4, sinh ay, -— (31) 
a), cosh aj, dz, cosh a», a3, Cosh a3, 4, cosh ay, 
b,, sinh ay, b., sinh a, 3, sinh a3, 04, sinh ay, 


By a manipulation of the elements of the above determinant, it may be reduced from one of the fourth order to one 
of the third order. 


a, tanh az, — a, tanh a, d2, — a), ob, tanh a2, — 5), tanh a, 
Q, a3, tanh a3, — a, tanh a, 43, — a, 63, tanh a3, — b,, tanh a, = O (32) 
ay, tanh ay, — a, tanh a, 4, — @), 04, tanh ay, — 0, tanh ay, 
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where 


Q, = cosh a, cosh ay 


-n 
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cosh a3, cosh ay, (32a) 


It will be recalled that the quantities forming the elements of the above determinant come from the solution of 


the quartic equation, Eq. (13), and that three of the four roots are dependent upon the critical load. 


The pro- 


cedure for obtaining the numerical value of the critical load involves assuming a value for the critical load, and, if 


the value happens to be correct, determinants (31) or (32) will be satisfied. 
value is not the correct one—i.e., the determinant yields a value other than zero. 


In general, the assumed critical 
If three or four values of the 


critical load are assumed and the corresponding values of the determinant are plotted against the critical load, then 


the value of the critical load may easily be obtained by graphical interpolation. 
In the foregoing analysis, it is assumed that the plate buckles symmetrically in the loaded (y) direction. 


There 


will be certain values of the plate aspect ratio (ratio of height to width), for which the sandwich plate buckles un- 


symmetrically. 


In these cases, the expressions for “, v, and w will take the following form. 


u = [C,,a), sinh a, + Co,d2, sinh as, + C3,d3, sinh a3,n + C4,d4, sinh as,n] sin (mw/2)E (33) 
v = [Ci,)1, cosh ay, + C2,b2, cosh ae,n + C3,b3, cosh a3,n + C1,b4, cosh ays,n] cos (nmw/2)é (34) 
w= [C,, sinh a, + C2, sinh az,n + C3, sinh a3,n + C;, sinh a4,n] cos (m/2)é (35) 


The boundary conditions remain unchanged, and their enforcement leads to four simultaneous homogeneous 


equations as it does in the symmetrical case. 


Following the same line of reasoning as before, the determinant 


corresponding to Eq. (32) giving the critical load is obtained. 


a, /tanh ae, a,,/tanh a,, G2, — Qi, 62,/tanh a», — b,,/tanh a, 
S, | a3,/tanh a3, — a,/tanh a, d3, — di, 53,/tanh as, — b,/tanh a, = 0 (36) 
ay,/tanh ay, a@,/tanh a, 4, — a, 54,/tanh as, — 6,/tanh a, 


where S, = sinh aj, sinh ae, sinh a3, sinh ay,,. 


For the case in which the loaded edges are simply supported and the unloaded edges clamped, an exact solution 


may also be obtained. 
will be discussed. 
respect to the coordinate system. 
supported edges. 


The procedure is similar to that described above and only the more important differences 
It can be assumed that the rectangular sandwich plate maintains its original orientation with 
The edge load, however, is made to act in the x or £ direction 


i.e., on the simply 


Under these conditions, the three equations of equilibriums, Eqs. (4), (5), and (6), remain unaltered except that 


the last term in Eq. (6) is replaced by (P.,/b) (0°w/0€?). 


fixed case [Eqs. (7), (8), and (9)] are again assumed as solutions. 


The displacement functions assumed in the loaded edges 


By following the procedure previously described, 


the quartic equation corresponding to that given in Eq. (13) is obtained. 


(fn + ] A,(1 _ v) | | —2A,B,f,3 —_ (2A, + Bij," f [2A \(n3 2)*r) fn + A(a 7, 2)?} = () 


From this point, the analysis is identical to that pre- 
sented for the loaded edges fixed case. Even the 
boundary conditions need no alteration because the 
sandwich plate was allowed to maintain its original 
relationship with the coordinate system. Hence, the 
determinants required to vanish are given by Eqs. 
(32) and (36). In reference 11, a solution is given for 
this problem. Lower bounds are obtained by Leggett’s 
method; upperbounds, by Galerkin’s method. 


NUMERICAL RESULTS 


Values of the buckling load, calculated according to 
the present exact theory, are presented in Figs. 2 
through 5. In Figs. 2 and 3, a buckling coefficient, C, 
is plotted as a function of the sandwich plate height to 
width ratio (aspect ratio). Each curve is valid for a 
particular value of the parameter R (relative rigidity 
in shear to that in bending) and a particular value of 
core to face thickness ratio. 


(37) 


The buckling coefficient, C, is defined by 
C = ocr! Foer, 


Physically, C is the ratio of the actual critical stress of 
the panel to the critical stress of a sandwich plate 
whose edges are simply supported, its core being infi- 
nitely strong in shear and without bending resistance. 
For values of the aspect ratio approaching infinity, C 
approaches 1, provided R is relatively high. 

Fig. 2 shows the relationship between C and the as- 
pect ratio for values of R equal to 50,000 and 25,000. 
The values chosen for c/t, the core to face thickness 
ratio, are 15 and 25. Differences in C due to the change 
in c/t were so small that a single curve is shown for each 
value of the relative rigidity R. 

The results of March and Smith (strain energy solu- 
tion) are also presented for comparison. In converting 
to C the values of the critical load coefficient, k,,, used 
by March and Smith in reference 3, the relationship 
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C = (3/s")he 


In deriving the above expression, D, 


was employed. 
the plate bending rigidity, was assumed equal to 
k(h® — c®)/12(1 — v*), which means that the contribu- 
tion of the core to the plate bending rigidity was neg- 


lected. 

In Fig. 3, curves of the same type are presented for 
values of R equal to 15,000 and 1,000. For each value 
of R, curves are presented corresponding to c/t ratios 
of band 25. As R decreases, the effect of the shearing 
deformations becomes noticeable. Note the changes 
that take place in the buckling coefficient at values of 
R between 15,000 and 1,000. At R equal to 1,000, the 
results are almost independent of the sandwich plate 
aspect ratio, and smooth curves are drawn rather than 
the ‘‘scalloped”’ characteristic curves corresponding to 
larger values of R. 

To extend the usefulness of the curves, the numerical 
results of Figs. 2 and 3 were cross-plotted to yield curves 
of the buckling coefficient C versus the sandwich plate 
relative rigidity R for selected aspect ratios and a fixed 
core to thickness ratio. These results are presented in 
Figs. 4 and 5. Also shown in these figures are curves 
corresponding to infinite aspect ratio (all edges simply 
supported) obtained by cross-plotting the results pre- 
sented in reference 7. At values of the aspect ratio of 
3.6, it was found that the values of the buckling coef 
ficient were within 6 to 8 per cent of the values given for 
the simply supported case of reference 7. 


COMPARISON WITH EXPERIMENT 


The theory presented in this paper was checked 
against the results of experiments conducted by the 
Forest Products Laboratory.'* The critical stress was 
computed by the present theory for 42 sandwich plates 
whose faces were composed of aluminum-alloy sheet 
and whose cores were of either balsa or cellular cellulose 
Only those panels that suffered a glue or 
Fig. 6 


acetate. 
compression failure were disregarded. shows 
the comparison of the theoretical results with the ob- 
served values. 

On the whole, it is considered that the theory is in 
reasonable agreement with the experimental results. 
Although three panels indicated a difference between 
theory and experiment of over 20 per cent, three- 
quarters of the panels agreed within 15 per cent, and 
one-half of the panels agreed within 10 per cent. Fig. 
§ shows a tendency for the theory to be conservative 
i.e., theoretical critical stresses lower than those ob- 
served experimentally. 

A discrepancy of 10 to 15 per cent between theory 
and experiment is judged to be reasonable on the basis 
of the many factors that affect the observed critical 
load. One of these factors is the unavoidable variations 
in the mechanical properties of the face and core ma- 
terial. In computing the theoretical critical stresses, 
the modulus value used for the faces was taken from 
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data published by the Aluminum Company of America 
for 24S-T Alclad aluminum alloy. For sandwich 
panels buckling at stresses exceeding the elastic limit 
of the aluminum alloy, the Young’s modulus in the 
The 
values recommended in reference 12 by the Forest 


equations was replaced by the tangent modulus. 


Products Laboratory were taken as the shear moduli of 
the core materials. 

There are factors that would tend to make the ex- 
perimental critical loads fall below theory. Probably 
the most important factor is the inability to obtain 
perfect rigid end fixation in practice. On the other 
hand, experiment is apt to be higher than theory if the 
unloaded edges are not perfectly hinged and if the core 
absorbs part of the edge load (in the theory, only the 
faces are assumed to absorb the edge load). In ad- 
dition, there was a '/4-in. overhang on each side of the 
simple supports, which also contributed toward in- 
creasing the observed critical values over those ob- 
tained by the theory. 

The scatter of 10 to 15 per cent mentioned should be 
compared with the scatter in the experimental results. 
For example, there were three pairs of panels for which 
all the geometric and physical properties within each 
pair were identical with the exception of the core thick- 
ness. With thickness, the critical 
load would be expected to decrease. In one of the 
three cases mentioned, the experimental critical load 


decreasing core 


remained constant, and in the other two cases the ex 
perimental critical load rose approximately 10 per cent. 
In the case of one pair of identical panels, there was a 5 
per cent difference in the observed critical loads; in a 
second identical pair, the difference was 30 per cent. 
It is not the intention of the author to cast aspersions 
on the experimental results. On the contrary, it is 
appreciated that reliable plate buckling tests, especially 
on sandwich plates, are difficult to perform, and it 1s 
considered that the experiments of reference 12 were 
well conducted. The comments in the foregoing para- 
graphs were solely intended to justify the conclusion 
that the theory was in reasonable agreement with ex- 


periment. 


CONCLUDING REMARKS 


An exact theory has been presented for the buckling 
loads of flat rectangular sandwich panels with loaded 
edges clamped, unloaded edges simply supported, and 
loaded edgewise in one direction. Briefly indicated is 
the method of solution for the case of a sandwich plate 
whose loaded edges are simply supported and unloaded 
edges clamped. 

A series of design charts are presented which give 
the buckling loads of sandwich plates for various rela- 
tive bending rigidities, aspect ratios, and core to face 
thicknesses of 15 and 25. It is appreciated that many 
plates may not fall within the range of parameters 


covered by the curves. The parameters were chosen 
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with a view toward encompassing most plates en- 
countered in practice without unduly lengthening the 
time required to perform the numerical computations. 
A check of the theory with published test data indi- 
cates reasonable agreement between theory and ex- 


periment. 
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Thermal Stresses in Conical Shells 


J. H. HOUTH* 


The RAND Corporation 


Shiite 


Graw SUMMARY 


This paper considers the stresses in a conical shell of revolution 


tes i) 3 é : 
es resulting from aerodynamic pressure and thermal gradients. In 
A-12, particular, the results are applicable to the conical tip of a missile 

flying at zero angle of attack; a numerical example is given for 
ng of a temperature distribution typical of a laminar boundary layer. 
SSUOR It is concluded that from the standpoint of eliminating high 


sl . 
n ped, stresses, large temperature gradients both along a ray and across 


the shell thickness are most undesirable. 


NOMENCLATURE 





11-4 = constants of integration 
D = Eh3/12(1 py?) 


modulus of elasticity 
= Poisson’s ratio 
shell thickness 
\ . \ (Eh 
= cone half-angle 


D) ctn? a 


= aerodynamic pressure 








\,, Me = bending moments per unit length of shell section 
V,, Ne membrane forces per unit length of shell section 
5 in p = density of shell material 
ical w = angular velocity of shell 
= distance along cone element measured from vertex 
(2 Q shear per unit length of shell section 
sity VY, WV = components of displacement in the direction of 
increasing / and the inward normal to /, re 
jer spectively 
. t = outward radial displacement 
I, = axial displacement 
a a = thermal coefficient of expansion 
cal 1 = temperature 
I-( = temperature levels 
Al = mean shell temperature (at any point /), above a 
reference temperature 77 
67) = temperature differential across the shell at any 
point (assumed linear) 
»¥ = dummy variables 
- i = V-1 
r ( d? d ') : 
L an operator defined by: {/ T x= L(y) 
di? dl l 
n AV il 
m = 1” 
z = 2Vi1 
I, Y2 = Bessel functions of the first and second kind and 
of the second order 
Z\(x = ber(x); Z(x) = bei(x); Z; = —(2/r)kei(x); 
Z(x) = —(2/7)ker(x), where the ber, bei, ker, 
| and kei functions are Thompson’s functions. 
Z'(x) refers to the derivative of the function 
with respect to x 
Received November 24, 1952. 
* Missiles Division. 








INTRODUCTION 


—- HEATING along the forward portion 
of a high-speed missile often results in large 
temperature gradients. In many cases such tempera- 
ture variations may have a pronounced effect on the 


structural design of a missile. 








Fig. |— Forces and bending moments acting on an 
element of a conical shell of revolution 





e 
shell surfoc 
e on outer 

Temperotur 











\. 10 in. 
——_ —__—— - —_ - ——__ 
7é in — 





Fig. 2 — Temperature variations along the inner 
and outer shell surfaces 
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The case of thermal stresses in cylindrical shells ANALYSIS 


has been treated by Timoshenko and others," ? while 4 , 
In reference 3, Meriam has restated the general 


equations from the bending theory of shells specialized 
to rotating bodies of revolution. With the inclusion 
conical nose of a supersonic missile at zero angle of of appropriate temperature and aerodynamic pressure 
attack. <A further application can be made to aircraft terms they appear as follows (refer to the Nomen. 


the present paper considers conical shells. In par- 


ticular, the analysis will consider conditions in the 


gas-turbine tailpipes. clature and Fig. 1): 
(d/dl)(Nil) — No + pw'/? sin? a = 0 (] 
(d/dl)(Qil) + Nactn a — pw*/* sin acos a+ pl = 0 (2 
(d/dl)(M jl) — Q, — M, = 0 (3 
V = Fe "(V — Wet (1 + v)a,AT 
l — ( saa 5 — a; { 
! a — y?) dl / ctn @) Via 
We a Eh | dV 4 I 7 W et a+ AT : 
y = ( —_— . — y P (+) 
4 (1 — vp’) v dl / ctn @) via 
uv D Ke in vy dW 14 a (dT) , 
‘ , = -_ ( (0 
di?" 1 dl glee 
7 p| ew 4 | div (14 a (67) - 
4 = —_ ( \/ 
; "ae 1 di "4 


Here it is assumed that temperature and pressure depend only on /. This would be the case for a missile at 
zero angle of attack. 

These equations may be handled in the standard manner developed by Meissner‘ and illustrated by Meriam.' 
In particular they lead to the following fourth-order differential equation for (Q,/): 


‘(0,1 1*(Q,1) K(1 a, d(6T kha, 1?AT PAT 
r di(Q1) + yf! (Cr +MOD = — (1 + via, d(é6T) ‘ ho EB ign ne | be 
dls dl’ tan” a dl tan a dl? dl 
3pw"(3 + v) sin acos a = [I(/ (S 


Designating the particular integral of Eq. (8S) by (Q,/),, the complete solutions in terins of five constants of integra 
tion may be expressed as follows: 


») ») ») ») 
(Ol) = Ai(Z, 4 * Zi!) + As (z. - 4) " As ( Zs + = Zi!) + Ai(Z, ; = fi) + (Ql), ( 
Xx \ x x ; 


x 
(NJ) = A; tan a — (Q,/) tan a — (pl*?/2) tan a (16 


. ene pr d _ 
No = pw’/* sin? a — , tana — tana (Ql), — tan a} A, (2, 
al 


5) | 2 | 2 f 7 ox? 
A: (Z' - Z. + -, Zi’) As (2 Zs ai) + A (2 a 2) | Yu 
— 2 oF X i —" ss x 
WV Du r a (67°) 4 —' it es 3 4 — | 4 D , As tan- a 
ee - p) - (2+ p)/ t: ¢ (33 y) sin- « - ( — pv) 
v h 1 in a Eh 1 a 1 Eh 
d(AT) d?(AT) tan a 4d ' y Le 3D 
D(1l + via, tan a dl + Da,/ tan a Al? - D Kh dl [L(Qi) p] + / L(Qil) » ¢ OKI (1 + v)p tan a 
2 V4 g 2(1 — ¥) 9  £(f = ” 21’ | 14 | Zz 4 2(1 — "> 1(1 —- a | 
y y 1 <2 xX 227 re “1 412 “1 x aT te “2d 


24 ay 211 — vy). 1H1—~yp),, = 2(1 — pv). 4(1 — ») _ th 
y ‘ie Zé — _ Zs 4. Y? 23 + A, = 23 + < Zs + x? Zs ( (12 
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M, = (d/dl)(M I) — (Q,) (13) 

dW/dl = [(l/DQi — v*)\@vM, — My) — ( h)al6T) (14) 
W = (1/Eh) tan a(vN,; — No) + V tan a — aJAT tan a (15) 
¢ = Vsina — Weosa (16) 

§é= Vcosa+ W sina (17) 


¢ = a/AT sin a + (/ sin a/Eh)(Ng — vN,) (18) 


In the general case of a truncated cone the five constants of integration could be evaluated by specifying either 
M,and Q, at both ends and NV, at one end, or appropriate combinations of stresses and displacements. However, 
in the following, the discussion will be limited to a complete cone. In this case A;, Ay, and A; are all zero. 

If AJ and 67 can be sufficiently well represented by a polynomial of only a few terims then clearly a series 
In the case of transition to a turbulent boundary layer, 


expansion for (Q,/), can be easily obtained from Eq. (8). 
However, as pointed out by 


the temperature variations may be too rapid to permit such a direct approach. 
Meriam, Eq. (8) may be factored into the two Bessel equations, 


, (at dy ; n* 
n° ~}+n + (nm? — 4)¥ =— A) = MN) (19) 
dn? dn ir? 
d*(Q,1) d(Q,1) —m?* 
n” e +m ‘ + (m*? — 4)(Q/) = —— ¥,(m) (20) 
dm* dm id? 
The particular integrals of these second-order equations are given by Kamke® and are of the form, 
r_ = T ye ; 
y, = ; V.(n) Jo(s)ITy(s) ds — : Jo(n) V2(s)F1,(s) ds (21) 
F Je S§ Z ( Ss 


Consequently, the given temperature distribution may be approximated, for example, by parabolic segments, 


and the resulting particular integrals may be computed with the aid of Eq. (21). 











































































































Dimensions of shell and properties ‘| F 
| of shell material 
s° in | 
[re th Ft 
90 — 7 4 04 —— 
———— Gi | a 
Z ——y | _ 
fs | iin | 
s v= + (Poisson's ratio) ' — 
- | c 
5 80 E= 30x 10° Ib/sq in. = Oo ) 
re a 
8 Outer edge of the shell is “ 
id assumed clamped so that 2 
= , av. 7 | o 
2 Crag 0 at £8, | $-00 - 
» 70 — _—_——— = 
7 se 
° | a — 
c< 7 
= ° 
¢ | £ -.08 -8 
* 60 —_—_— ~—_—- 2 
= al = 
3 | Pa s 
) tc 4 
. } > =-12 —— * 
a | - 
: a i 8 
¢ 50f——__+____1_____ : 3 
3 Temperature distribution 2 = 
2 ra 
é — fT o 716 -16 9 
2 ISOO'F / ,-T Inner 5 =~ 
5 aol 7 a 7 be = 2 
s ae : 
: . E -20 - 20 
= ° a * = 
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Shell assumed stress free at 
reference temperature of |OO F | -24 -24 
L L | 
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Lyin. 
Fig. 3a— Bending moment, M,, as a function Fig. 3b— Membrane force, N, ,as a function of 2 


of 2 in a heated conical shell in a heated conical shell 
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Fig. 3c—— Membrane hoop force, N, ,as a function of 2 
in a heated conical shell 


EXAMPLE 


Consider a temperature distribution as shown in 
Fig. 2. Such a variation might occur in the case of a 
laminar boundary layer and is given analytically by 


Eqs. (22). 


6T = C— [b+ Kil — 1,)?] 
AT = 1/2[C+6+ Kil — 4)]? (22) 
K = (a — D)/I,’ \ 


For this simple case (and for a nonrotating shell in 
the absence of aerodynamic pressure), the particular 
integral is given by Eq. (23), 
0), = ke :+ ree) paix 
h 
[2A ees tana — 2KD(1 + eet (23) 
h 
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The shell was considered as clamped so that the 
boundary conditions are as follows: 


l=h, ¢=0, dW/dl =0 (24) 


The principal results for a typical choice of a-c are 
shown in Fig. 3. 


CONCLUSIONS 


The abrupt changes in behavior for / > 110 apparent 
in all figures is, of course, the effect of rigid clamping. 
In an actual missile configuration, rigid clamping 
certainly will not be achieved, so that the case under 
consideration represents an extreme. Moreover, the 
effects of clamping will always be extremely localized. 
The important thing is the appearance of a high 
bending moment at the tip of the nose. This bending 
moment is the result of the thermal gradients and will 
be affected little by any reasonable variation in clamp- 
ing conditions at the outer edge. It is apparent that, 
if high stresses are to be avoided, then large temperature 
gradients must be eliminated. 

The equations provided are sufficiently general to 
include aerodynamic pressure and rotational effects 
should this be found desirable in any given application 
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Diffusers for Supersonic Wind Tunnels 


J. LUKASIEWICZ* 


National Aeronautical Establishment, Canada 


SUMMARY 


A systematic review of published experimental data relating to 
performance and design of wind-tunnel diffusers for entry Mach 
Numbers up to 10 is presented. It is found that: 

(1) Compressibility effects are negligible in diffusers for M - 
0.9. 

(2) Shock compression is more efficient in constant-area than 
in divergent ducts and usually requires appreciable duct length 
for completion. 

(3) Pressure recovery of diffusers without contraction is im 
proved by use of a long throat at entry but is always smaller 
than the theoretical normal shock recovery 

(4) At Mach Numbers below 4, maximum or starting con 
traction of constant geometry diffusers is closely predicted by 
simple theory, and better than normal shock pressure recoveries 
are achieved. 

(5) With variable contraction diffusers, pressure recoveries of as 
much as twice the normal shock recovery are obtainable at 
hypersonic Mach Numbers. In general, optimum diffuser con- 
traction angles increase with entry Mach Number increasing. 


SYMBOLS 


D = diameter or side of square 

D = hydraulic diameter, 4(area/perimeter) = 
D (for circular or square sections ) 

L = length of shock compression region 

M = Mach Number 

M; = Mach Number based on static pressure and Fanno 
line of state 

P = absolute pressure, static 

r = absolute pressure, stagnation (at nozzle or diffuser 
inlet) 

P,’ = stagnation pressure at diffuser outlet 

Rey, = Reynolds Number based on free-stream velocity, 
density, viscosity, and length D 

0 = total angle of divergence or contraction (all angles 
quoted in the text are total angles) 

y = diffuser contraction ratio = 
(diffuser entry area) /(diffuser throat area) 

Ymar, = Maximum theoretical diffuser contraction ratio (for 


flow starting ) 


INTRODUCTION 


i es SUBJECT OF THIS PAPER was first reviewed by 
Crocco! in 1935, At that time little experimental 
data existed on supersonic diffusers, and most results 
quoted by Crocco related to shock compression in di- 
vergent ducts, such as Laval nozzles or subsonic dif- 
fusers. A second review? appeared in 1946; it included 
newer experimental data and dealt with the theory of 
supersonic diffusers. Since then scores of supersonic 
tunnels of various designs were built, and enough ex- 
perimental data were accumulated to warrant the 
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present attempt at a systematic treatment of super- 
sonic diffuser design and performance. While in no 
way complete, it indicates the extent of available data 
and gaps to be filled in. 

The subject has been divided up into a number of dis- 
tinct problems, as follows: 

(1) Performance of subsonic diffusers at high sub- 
sonic Mach Numbers (i.e., effects of compressibility). 

(2) Shock compression in a simple divergent and 
constant-area duct as a phenomenon common to all 
supersonic diffusers. 

(3) Performance of diffusers without contraction and 
influence of a long, constant-area section at entry. 

(4) Maximum contraction of constant geometry dif- 
fusers and their performance. 

(5) Contraction and performance of variable geom- 
etry diffusers. 

Only diffusers used with fully closed wind-tunnel 
working sections are considered, with no models (un- 
less otherwise stated) mounted in the working section. 
Insofar as possible, only reliable experimental data are 
quoted; in fact some well-known results, obtained 
with wet air or with an ill-defined diffuser setup, were 
omitted from this review. 

In all cases the Reynolds Number of tests is given in 
terms of free-stream conditions at diffuser entry Mach 
Number and hydraulic diameter D of diffuser entry. 
SUBSONIC DIFFUSERS 


COMPRESSIBILITY EFFECTS IN 


Although outside the scope of supersonic diffuser 
problem, the operation of subsonic diffusers at high 
subsonic velocities is briefly considered as a particular 
aspect of compressibility effects in diffusers in general. 
It is found that the performance of diffusers of small 
divergence angle (< 7° total cone angle) is practically 
independent of entry Mach Number (provided M < 
0.9) when expressed as a fraction of theoretical pressure 
rise. This is shown by experimental results plotted 
in Fig. 1; about 85 to 90 per cent of theoretical pressure 
rise is obtained in practice. The former corresponds 
to a stagnation pressure loss of 3 per cent at M = 0.6 
and 6 per cent at 7 = 0.9. Since these results were 
obtained with no model in the tunnel working section 
and with a particularly good velocity distribution at 
the diffuser entry, they must be regarded as upper 
limits of efficiencies attainable in actual subsonic wind- 
tunnel diffusers. 

In some high-speed subsonic tunnels, an adjustable 
sonic choke or throat is fitted between the working sec- 
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tion and the diffuser entry. Its purpose is to provide 
control over working section Mach Number and/or 
to prevent disturbances originating in the diffuser or 
compressor from reaching the model. Following the 
choke, an expansion to slightly supersonic velocity and 
shock compression take place in the subsonic diffuser 
entry. The diffuser data just reviewed are therefore 
not applicable to subsonic tunnels equipped with a 
choke which, insofar as diffuser is concerned, are equiva- 
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lent to supersonic tunnels operating at low Mach Num- 
bers. 


SHOCK COMPRESSION IN Ducts 


Experiment shows that, in general, a continuous de- 
celeration of supersonic flow to subsonic velocity is not 
possible, the transition always occurring through shock 
compression. The corresponding simple theoretical 
model of flow postulates a normal shock but neglects 
viscous effects. It is not surprising, therefore, that 
shock compression as observed in ducts, in presence of 
wall boundary layer, usually differs from the theoretical 
model. Observations of shock compression in diver- 
gent ducts (such as Laval nozzles, downstream of 
throat) show a gradual rather than sudden increase of 
pressure at the wall, the static pressure never attaining 
the values prescribed by simple theory and stagnation 
pressure losses appreciably exceeding those appropriate 
to anormal shock. Typical of the latter are results of 
Castagna’s® measurements of shock compression in a 6° 
conical pipe, shown in Fig. 9 in terms of Mach Number 
based on the lowest observed wall pressure (1.e., at the 
start of compression) and assumption of isentropic flow. 
Comparison with the normal shock curve indicates the 
inadequacy of the normal shock theory. However, 
the discrepancy is reduced, particularly at higher Mach 
Numbers, by making the angle of duct divergence 
smaller, as shown in Fig. 2. 

The above overall characteristics of shock compres- 
sion in divergent ducts are supplemented by schlieren 
observations, which indicate that, as a result of shock 
and boundary-layer interaction, the region of shock 
compression may involve several curved or oblique 
shocks, forming regular or Mach-type reflections, and, 
at the walls, extremely thick boundary layers and flow 
separations. The latter increase both with shock 
pressure ratio (i.e., flow Mach Number) and with duct 
divergence, approaching in extreme cases an essentially 
free-jet configuration, the flow being completely de- 





Detached flow in supersonic nozzle for M = 4.4. 
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DIFFUSERS FOR 


tached from the walls following shock compression. An 
example of this is shown in Fig. 3, which is a schlieren 
photograph of flow in a +.4 Mach Number tunnel nozzle 
operating at a pressure ratio insufficient to induce a 
fully supersonic flow. Both schlieren and pressure ob- 
servations suggest that higher efficiencies of shock com- 
pression are obtained in small divergence or constant- 
area ducts. In such ducts there exists in the subsonic 
flow, downstream of the shock compression region, a 
positive velocity gradient that tends to stabilize the 
boundary layer and reduce the separation and resulting 
losses. 


A number of observations of supersonic compression 
in constant-area ducts, in fact, indicate that this is 
the case, although the compression does not necessarily 
resemble a normal shock. Typical schlieren photo 
graphs taken in a l-in. square duct are shown in Fig. 4. 
The compression consists of several shocks closely fol- 
lowing one another, their shape presumably depending 
on boundary-layer thickness relative to duct diameter 
(at given Mach and Reynolds Numbers). 


A similar flow configuration may be expected to oc- 
cur in shock-tube flow and is believed to cause the 
hitherto unexplained appearancet of wave-speed cam- 
era photographs, examples of which are shown in Fig. 
6.t The wave-speed photographs are taken with a 
schlieren optical system and a drum camera, the flow 
inside the tube being viewed through narrow and long 
windows parallel to the tube axis. The photographs 
can be interpreted with the aid of Fig. 5, in which shock- 
tube flow is represented diagrammatically in x, ¢ co 
The region of tube flow corresponding to 
the tube 


ordinates. 


the photographs is marked in Fig. 5, and 
t First evident in original photographs obtained by Payman 
and Shepherd, reference 8, Figs. 3, 4, and 8 
t The author is indebted to Dr. I. I. Glass and the Institute of 
Aerophysics, University of Toronto, for permission to make use 


of wave-speed photographs originally published in reference 9 


Schlieren photographs of shock compression in 1-in 
square duct (Rep ~ 0.2 X 10°) in presence of thick boundary 


Fic. 4. 
c — M ~1.5, shock compression at increasing 
d- Mw~2.5 


layer. a, 5, 
distances from nozzle 




































SUPERSONIC WIND TUNNELS 619 
INTERACTION OF 
REFLECTED SHOCK 
WITH TUBE 
BOUNDARY LAYER 
RAREFACTION WAVE 
\ 7 ry y 3 
\ A fy Be / ' 
4 ~ 
om \ ie, 
\y \ \ i ~~ 
“\ \ \ a 
/ \ \ : 
K y 4 m | 
\ \ \ ee 
\ \ "| 
! 
\ H 
\ ce 
\ // Bye 
WZ 
Wy. 
WU 
x 
DIAPHRAGM REGION OF PHOTOGRAPHS 
=: SSS=—}« 
4H | @H Low 
'nwiriat AER PRESSURE 


Fic. 5. x-t diagram of shock-tube flow 























Wave-speed camera schlieren photographs of shock 


compression in a shock tube.’ 


length (x) and time (¢) coordinates are oriented in the 
same manner in Figs. 5 and 6. Following the bursting 
of the diaphragm, a shock wave is propagated in the 
low-pressure chamber of the shock tube and eventually 
is reflected from the closed tube end, as shown in Fig. 
5. On its return travel, the shock enters tube sections 
in which boundary layer has already developed and 
presumably causes interaction similar to the stationary 
shock-boundary-layer system of Fig. 4. The multiple 
shocks, sketched in Fig. 5, appear as series of oblique 
ripples in the wave-speed photographs, Fig. 6. It will 
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TABLE 1 Aes 
Shock Compression in Tube in Presence of Thick Boundary Layer (Computed from Experimental Data of Reference 12) iF 
1. Station Nozzle exit or 1 2 3 4 5 6 7 8 i 
tube entry 6k 
2. L/D 0 2 3.75 8 14 21.5 26.4 30.5 40 
3. (P;/Po) X 103 observed 2.2 2.0 2.9 3.5 4.4 5.5 6.6 7.0 10.4 
(shockless flow ) r 
4. Myanno (based on nominal 5.0 4.28 3.94 3.53 3.08 2.69 2.39 2.19 1.95 
conditions at tube entry, 
1 = 5, Pi/Po = 0.00189) 
5. M (based on measured P;/Po 4.87 4.76 4.64 4.49 4.31 4.14 4.0 3.9 3.67 cr 
and isentropic flow ) 
6. Ps/P, measured fe 19.6 16.6 13.4 10.1 7.6 6.1 5.1 
4; e/a teens We 21.2 17.9 4 10.9 8.3 6.5 5.4 
8. (P2/Pi)m i 26.3 24.9 23.3 21.5 19.8 38.8 38.4 0°} 
9. L/D (length of shock com- a 12 2.5 12 11.4 10.8 10 9 et 
pression ) 

6 b 
be noticed that the ripples are approximately parallel, the low Reynolds Number of tests (about 0.2 million ‘rT 
indicating that the whole compression region moves based on tube diameter), appreciable compression oc- 
with a substantially constant velocity. The phe- curred in the tube even for shockless flow, the static 
nomenon appears to gain strength with the initial tube- pressure P, (third line, Table 1) increasing about four _r 
pressure ratio (given in Fig. 6) or shock intensity and, and one-half times between the nozzle exit (or tube en- 
in fact, is hardly noticeable at extremely low tube- trance) and station 8, 40 diameters downstream. Since °° 

- ” a 
pressure ratios.” the flow was, to a good approximation, adiabatic, its 5 
-_ ‘ ee , effective mean parameters could be obtained from pres- 
So far, only few quantitative results are available on ene ? : : ; : FG 
eect it sure distribution and Fanno line of state (derived from n 
shock compression in constant-area ducts. Most ex- : : : ; 
: ‘ Perey ‘ equations of conservation of mass, energy, and area). 
tensive tests were carried out at M.1.T.!° with relatively Hs 
; : ies In order to assess the boundary-layer thickness, the 
thick boundary layers and therefore under conditions  ,._. k's ad ‘ ee eee shoc 
Rethas tes : ; Fanno line’? Mach Number 7/7; is compared in Fig. 7 
not representative of conventional wind-tunnel prac- aay ‘ee : is re 
z and Table 1, lines 4 and 5, with Mach Number J1/ de- 
tice. : : ; a i at 
rived from static pressure ratio P;/P 9 on the assump- 
as = ‘ : ae a inert 
rhe M.1.T. flow setup consisted of a circular super- tion of isentropic flow. The corresponding boundary- \ 
sonic nozzle, followed by a 50-in. long 1l-in. diameter layer thickness can be roughly evaluated from the fol- ne 
; ‘ ‘ ve : a ‘ ‘ y* the 
tube in which wall-pressure distribution was measured lowing consideration: Assuming that boundary layer 
s : ae : Rae aye 
for different back pressures. Typical results obtained does not fill completely the tube cross section, the isen- h 
. mn —— : : : . : : the | 
in M.I.T. tests are analyzed in Table 1. Because of | tropic Mach Numbers 1/ and Mach Number-area ratio 
: - ; : lap: 
relationship can be used to estimate the boundary-layer | 
: . : Nie tap é nor! 
displacement thickness at any station. Taking a one- 
, aad . : resp 
seventh power-law velocity distribution in the boundary 
sk /s , . Tah areé 
P, layer and 6*/6 ~ O.4.at \/ = 4, we find that at station 7 frict 
40 Tal - °. ‘ wr ric 
é ee : ary layer 1S < 0. ; . 
COMPUTED FROM (Table | _the boundary layer is bout M4 in. thick, ‘sth 
30 NEUMANN & LUSTWERK = almost filling completely 0.5-in. radius tube. 7 t 
MIT. (1949) care _ 
With increased tube back pressure, a shock compres- whi 
20 — sion develops; pressure ratios P2/P, recorded for com- tail: 
15 = pressions that started at stations | to 7 are given in mot 
Table 1, line 6. Here P; is the pressure immediately ing 
10 preceding the compression, and P2 is the highest pres- flux 
8 sure attained downstream. As shown in the last line loss 
of Table 1, the shock compression region (defined as the 
\ - tube length from P,; to P:) is some 10 diameters long, F 
NORMAL SHOCK/“  \y, / ISENTROPIC and appears to increase in length with Mach Number. in | 
(THEORY) ™ FANNO FLOW ' a pre 
Observed pressure ratios are compared in Table 1, t 
te cod stat 
3 lines 7 and 8, with the ones corresponding to a normal thi 
lil 
2 + Incidentally, from the observed pressure rise in the tube and latt 
Mr it is possible to determine! the average friction coefficient lish 
. For entry Mr = 5, the pressure rises (at 30 tube diameters) diff 
to three and four-tenths times the value at entry. This corre- tha 
' sponds to \L/D = 0.07, whence \ ~ 0.0023. This value agrees ‘ 
te) 1 2 3 4 5 wm With d determined at Reynolds Numbers from 0.2 X 10° to to é 
Fic. 7. Shock compression pressure ratio correlated to Fanno 5 X 10° at L/D ~ 30 and at inlet Mach Numbers between 2 by 
line and isentropic Mach Numbers. and 4.11 SiO} 
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sion in constant-area ducts, at WM ~ 2. 
shock at Mach Numbers /, and \/. The agreement 
isremarkably good with respect to normal shock at J/, 
(at all tube stations), the discrepancy with (P2/P1)1 
increasing downstream. This is also shown in Fig. 7. 

Although, as seen in schlieren observations of flow, 
the shock compression in presence of thick boundary 
layers differs vastly from normal shock configuration, 
the above agreement in the overall pressure ratio is per- 
haps not surprising. The theoretical assumptions of 
normal shock transformation are a priort met with 
respect to conservation of energy (adiabatic flow) and 
area (constant tube diameter) so that, provided the 
frictional forces are small in the compression region (as 
is the case with separated flow), the normal shock solu- 
tion is strictly applicable at tube sections between 
which compression takes place, irrespective of the de- 
tails of shock transformation mechanism. Further- 
more, the ‘‘Fanno line’? Mach Number /,;, represent- 
ing a mean value that corresponds to the actual energy 
flux and accounts for the upstream kinetic energy 
losses in the boundary layer, is the correct one to use in 
the calculations. 

From the above consideration one would expect, 
in general, a good agreement between measured and 
predicted (on basis of ./,) shock compression in con- 
stant-area ducts, both with thick (Mr < M) and 
thin (M, ~ M) boundary layers. However, in the 
latter case experimental data are scarce. Most pub- 
lished data, here reviewed in the next section, refer to 
diffusers that include a divergent subsonic portion, so 
that pressure ratios measured across them are subject 
to additional subsonic diffuser losses, in turn influenced 
by the velocity distribution following shock compres- 
sion in constant-area entry. 


FOR SUPERSONIC 
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On the other hand, more data, if less conclusive, 
are available on length of shock compression region. 
It was obtained directly from the above tube experi- 
ments (Table 1, bottom line), as well as indirectly from 
tunnel diffuser tests, in which the length of constant- 
area diffuser entry was determined to give maximum 


pressure recovery. 


On the basis of tube experiments, a specific variation 
of shock compression length with Mach Number (\/,) 
was suggested,'? whereas compilation of results ob- 
tained at Mach Numbers close to 2 was to determine 
the variation of compression length with Reynolds 
Number.'* This is seen in Fig. 8, which indicates, as 
expected, smaller compression length L/D at higher 
Reynolds Numbers. small L/D results 
need further corroboration since they were obtained 
with M.I.T. tunnels, which show higher pressure re 
coveries than other tunnels of similar geometry (cf. 
Fig. 10). 


However, 


Whereas pressure ratio of shock compression may be 
adequately correlated to Mach Number J/p, the length 
and structure of the compression region may be ex 
pected to depend on 1/7, on Reynolds Number based on 
duct diameter, and on Reynolds Number based on 
boundary-layer thickness.f The last parameter was 
not included in the attempted correlations we just 


mentioned. f 


The information available on shock compression in 
ducts, although inadequate in detail, clearly indicates 
that (1) the efficiency of shock compression is higher 
in constant-area than in divergent ducts and (2) several 
duct diameters are required for compression to run its 
full course. The advantages resulting from inclusion 
of a constant-area section at entry or at throat of super 


sonic diffusers are examined below. 


DIFFUSERS WITHOUT CONTRACTION 


Under this heading come subsonic geometry diffusers 
operating with supersonic entry velocities. Their 
performance is given in Figs. 9 and 10 in terms of dif- 
fuser entry Mach Number and stagnation pressure 
ratio P,/P,’' across the diffuser. 

Least efficient are diffusers without a long constant- 
area entry section and with large divergence. Cas- 
tagna’s results (6° conical) and Stewart’s tests (about 
12°, two-dimensional expansion), Fig. 9, are in good 
agreement, but the M.I.T. result, Fig. 10, obtained at 
M = 1.75 with a square, 6° divergent (in both planes) 
diffuser, falls a long way below and lies in line with 


+ Provided a fully turbulent pipe flow has not developed. 

t The phenomenon of shock compression in ducts is analogous 
to flow at the base of bodies moving at supersonic speeds in that 
it involves interaction of shock waves with boundary layers and 
separated regions. The base flow phenomena were correlated 
in terms of ratio of boundary-layer thickness to base diameter by 


Chapman." 
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Fic. 9. Diffusers without contraction: pressure ratio 


other M.1I.T. tests made with a small divergence (2.66° ) 
conical diffuser. t 

The expected favorable effect of a constant-area 
section at diffuser entry was demonstrated by NACA" 
tests of a 10- by 3.84-in. wind tunnel, fitted with a J/ = 
2 nozzle, which was followed by 46 in. of constant-area 
section joined to a 5.5° (average, in both planes) rec 
tangular diffuser. In Fig. 11, the pressure ratios re 
corded for various distances of the start of the pressure 
rise (or shock compression) in the constant-area duct 
from divergent diffuser entry are given. The required 
pressure ratio is reduced from 1.86 to 1.65, witb length 
of compression in constant-area duct increased from 
17 to 46 in. 

In Fig. 10, pressure ratios for diffusers with a long 
constant-area entry are indicated by filled-in symbols. 
M.I.T. 
closely approach the ideal normal shock values and 


tests show smallest pressure ratios, which 
indicate the favorable effect of a constant-area entry. 
Comparable JPL,'? NACA," and German'® results lie 
much higher, even exceeding already mentioned M.I.T. 
data for diffusers without a long constant-area entry. 

Diffusers of the geometry here considered are not 
subject to large hysteresis effects. The pressure ratio 
required to obtain and maintain a given Mach Number 
in the working section is practically the same whether 

7 In general (cf. below), the M.I.T. tests show a smaller pres- 
sure ratio than obtained in other comparable experiments. 
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Fic. 11. 





Effect of length of constant-area entry on diffuser 
pressure ratio. 


reached from larger or smaller values.{ This was de- 

t In intermittent tunnels run by means of large volume reser 
voirs, the pressure ratio across a tunnel diffuser decreases during 
the run. With diffusers of the above type, the minimum start- 
ing and the actual flow breakdown pressure ratios are practically 
equal. 
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TABLE 2 
Maximum Theoretical Diffuser Starting Contraction Ratio and 


Pressure Recovery (y = 1.4) 

M W mex (Po/ Po’ )W. maz 
1.0 1.000 1.00 
1.2 1.028 1.00 
1.4 1.068 1.01 
1.6 1.120 1.04 
1.8 1.170 Li 
2.0 1.216 1.20 
2.5 1.316 1.59 
3.0 1.391 2.27 
3.5 1.445 3.33 
4.0 1.488 5.00 
1.5 1.520 7.14 
5.0 1.548 10.75 
6 1.577 21.79 
7 1.598 $1.15 
s 1.614 73.53 
8) 1.624 124.69 
10 1.632 202.43 
x 1.667 x 


termined for M.I.T. circular diffusers'* and also in 


Stewart’s'* tests; of the two values of pressure ratios 
shown in Fig. 9 at JJ = 5.5 and 6.5, the higher (by 
less than 10 per cent) ones were required to start the 


flow. 


DIFFUSERS WITH CONSTANT CONTRACTION 


Diffusers grouped in this category are characterized 
by a convergent entry section in which appreciable 
compression is achieved at supersonic speeds. They 
are therefore usually known as supersonic diffusers, 
The 


final compression to subsonic velocity occurs in super- 


in distinction to purely divergent, subsonic ones. 


sonic diffusers by means of shock compression of the 
type already considered near the minimum or throat 
diffuser cross section, on its downstream side. 

One of the fundamental design parameters is the dif- 
A,/A*’, where A; = dif- 
fuser entry or nozzle exit area and A*’ = diffuser 


fuser contraction ratio y = 


minimum or throat area. 

The maximum possible value of y at any entry dif- 
fuser Mach Number was derived theoretically? from 
considerations of starting the flow in a supersonic 
nozzle-supersonic diffuser system and corresponds to a 
normal shock at the diffuser entry Mach Number and 
sonic flow in the diffuser throat. The theoretical re- 


sults are given in Table 2, where War, denotes the 
maximum theoretically permissible contraction and 
(Po/P'\)Wmar. is the corresponding stagnation pressure 
tatio across the diffuser on the assumption of isentropic 
flow except for shock located at throat. 

Typical variation of pressure ratio with contraction 
ratio for constant geometry diffusers is shown in Fig. 
12. The top curve marked ‘‘starting’’ shows that 
supersonic flow at diffuser entry cannot be obtained 
for y > 1.4, irrespective of the pressure ratio applied. 
This limiting contraction closely agrees with the theo- 


t See reference 2. This criterion was first established and 


experimentally checked in Germany.” '® 
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retical War. The optimum contraction, at which 


smallest pressure ratio (~ 3) is required to initiate the 
flow, is only slightly smaller than the limiting contrac- 
tion. 

Also in Fig. 12, a curve is included for a variable dif- 
fuser, whose contraction can be increased after super- 
sonic flow at entry has started. Comparison of this 
curve with the one for starting or, in effect, constant 
contraction diffuser, shows that at the limiting contrac- 
tion of 1.4, a pressure ratio of only 2.4 is required to 
maintain the flow compared with the starting pressure 
ratio of 3. 

In Figs. 13 and 14, experimental contraction and 
diffusers are 
At M < 3.5, 
maximum contraction ratios for starting,{ Fig. 13, 
a result 


pressure ratios of constant geometry 


given in terms of entry Mach Number. 


agree closely with the theoretical Ya, values 
somewhat unexpected in view of the already considered 
and actual 


large discrepancies between theoretical 


‘normal shock’’ compression in ducts. The agreement 
is presumably due to the favorable influence on the 
boundary layer of the accelerated subsonic flow (during 
starting) downstream of the shock compression, in the 
diffuser contraction. In general, a smaller contraction 
is possible with diffusers having a long constant-area 
throat, because of additional losses in the throat sec 


tion. 


t Practically equal to optimum contraction (for minimum pres 
sure ratio), as pointed out above 
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At Mach Numbers between 7 and 10 results from 
only one set of experiments” are available and give 
contraction ratios much in excess of mar. No expla- 
nation of this discrepancy with results at lower Mach 
Numbers can be offered at present. 

The diffuser pressure ratios are given in Fig. 14 and, 
in general, closely follow the theoretical normal shock 
values. Hysteresis effects, already considered, are 
indicated in Fig. 14 by arrows, which connect pressure 
ratios required for starting and for maintaining the 
flow. Although long-throat diffusers require a smaller 


operating pressure ratio, their starting performance is 
similar to diffusers without a long throat. 


1953 


In square or rectangular tunnels the supersonic dif- 
fuser contraction is usually two-dimensional and the 
question of optimum contraction angle arises. Since, 
for a given contraction ratio, shock losses increase and 
frictional wall losses decrease with increasing contrac- 
tion angle, an optimum contraction angle might be ex- 
pected to exist. The limited data available? indicate 
that at 17 < 4 small contraction angles, of the order of 
3° to 4°, are favorable, higher pressure ratios being 
required with large contraction angles. f 
agreement with further evidence obtained from tests of 


This is in 
variable diffusers. 


VARIABLE CONTRACTION DIFFUSERS 


Following the establishment of supersonic flow, a 
further contraction is possible with variable diffusers, 
the pressure ratio at which the flow breaks down being 
much lower than the starting pressure ratio. In order 
to take advantage of variable diffusers, a tunnel drive 
capable of providing momentarily (for starting) a 
higher than operating pressure ratio is usually required. 
Intermittent tunnels run by means of pressurized 
and/or evacuated reservoirs lend themselves very well 
to this type of operation. 


Typical characteristics of variable contraction dif- 
fusers are shown in Fig. 12 at JJ = 2.83. Comparison 
of ‘‘starting’’ and ‘‘variable contraction” curves indi- 
cates that in the latter case much larger contractions 
and smaller pressure ratios are possible. The optimum 
contraction, at which a minimum pressure ratio is ob- 
tained, is significantly smaller than the maximum pos- 


sible contraction. 


Experimental results for variable diffusers are col- 
lected in Figs. 15 and 16. From Fig. 15 the maximum 
and optimum contraction and the effect of model 
mounted in the working section on optimum contrac- 
tion can be estimated. Diffuser pressure ratios are 
The minimum values indicated by 
smaller than the 


given in Fig. 16. 
the broken line are 
theoretical ones for diffusers having Wmar. contrac- 


somewhat 


tion. 


With regard to the optimum contraction angle of 
variable diffusers, it appears that it increases with the 
Mach Number. In the NOL tests! (diffuser without 
long constant-area throat), an optimum throat loca- 
tion was determined, while the total (including con- 
vergence and divergence) diffuser length was kept con- 
stant. The optimum throat position was found to be 
at 5.35D from the diffuser entry for JJ = 2.48 and 2.83 
and at 3.9D for \J = 4.92, where D is the height of 
square tunnel working section. The corresponding 


t Insufficient to separate Reynolds and Mach Number effects. 

t This statement is based on comparison of NOL data with 
GALCIT results, Fig. 14. The M.I.T. results show lower pres- 
sure ratios with large contraction angles but also indicate the de- 
sirability of a small angle. '* 
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DIFFUSERS FOR 


entry angles were 4.5° and 5.6° at optimum contrac- 
tions for the two lower Mach Numbers. 

The above results were confirmed by R.A.E. tests 
at MV = 2.48 with incomplete diffuser,** Fig. 17, which 
show that contraction ratio (and, hence, pressure re- 
covery) increases with contraction angle decreasing. 

Galcit tests,?! Figs. 15 and 16, indicate relatively 
small optimum contraction ratios and pressure recov- 
eries, Which can be attributed to an excessively large 
diffuser entry angle (~28°). 

On the other hand, apparently extremely low pres- 
sure ratios were obtained in M.I.T. tests at similar 
Mach Number, Fig. 16, in spite of a 20° contraction 
angle. It is not certain whether this was due to a long 
constant-area throat used in these tests. 

Both NACA and NOL tests indicate that at hyper- 
sonic speeds (say J > 5) much larger diffuser entry 
angles (between 13° and 18°) should be used. The 
following results were obtained at NOL” with two dif- 
ferent entry angles, at 1/ = 7.2, Rep ~ 2 X 10°: 


Diffuser entry angle a (g* 12.9° 
Optimum ¥ = 5.5 5.5 
Length to throat from diffuser entry = 6.6D 3.6D 
Subsonic diffuser angle = 3° 3° 
Min. Po/Po’ = 58 36 


In the same series of experiments it was found that 
the effect of subsonic diffuser divergence (up to 5.5~) 
was negligible and that a constant-area throatt (3D 
length) had no effect on optimum contraction ratio and 
minimum pressure ratio; at contraction ratios smaller 


+t As pointed out in reference 20, in hypersonic tunnels oper- 
ating at high stagnation temperature, it is advantageous to have 
a short diffuser throat in order to reduce regions of high heat- 


transfer rates and thus ease cooling problems 
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WIND TUNNELS 625 
than optimum, parallel-throat diffuser gave slightly 
(about 10 per cent) smaller pressure ratios. 

In two investigations the effect of Reynolds Number 
on diffuser performance was found to be small at hyper- 
sonic speeds. Bertram?* found no change in the mini- 
mum diffuser pressure ratio with the variation of 
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Fic. 17. Effect of angle of contraction on contraction ratio 








Rep from 0.5 to 3 X 108, at M ~ 6.86 and y = 2:8. 


The NOL results, summarized below, indicate increase 


with Rep in the maximum possible contraction but 
small change in the pressure ratio: 


M=7.2 To = 15°C. 


Po at. 5 10 20 30 
Rey X 10-8 ] 2 } 6 
y 6.6 7.5 9.2 10.4 
Po/Po’ 36 33.5 


CONCLUSIONS 
(1) About 8S per cent of theoretical pressure rise is 
cone) diffusers at 


> 


obtained in subsonic small angle (7' 
M < 0.9. 

(2) Shock compression in ducts may involve com- 
plete flow detachment in divergent ducts or flow sepa- 
ration in constant-area ducts with multiple shock sys- 
tems extending over several diameters of duct length. 
The latter type of shock compression can be observed 
in Wwave-speed camera records of shock-tube flow. 
Shock compression losses decrease with duct divergence 
decreasing and static pressure ratios observed in con- 
stant-area tubes correspond (at least in one case) 
to theoretical normal shock values taken at mean or 
Fanno line Mach Numbers. In general, the phenom 
enon of shock compression in constant-area ducts can 
be expected to depend on Mach Number and Reynolds 
Numbers based on tube diameter and boundary-layer 
thickness (for a given boundary-layer flow). 

(3) Pressure recovery of diffusers without contraction 
is improved by use of a long constant-area entry but is 
always smaller than the theoretical normal shock (i.e., 
Pitot tube) recovery. 

(4) At Mach Numbers below 4, maximum contrac- 
tion of constant geometry diffusers (or starting con- 
traction of variable diffusers) is closely predicted by the 
simple theory. At optimum contractions (which are 
practically equal to maximum contractions) normal 
shock pressure recoveries are achieved, although ap- 
preciably larger starting pressure ratios are usually re- 
quired at MJ > 2. Small angles of two-dimensional 
contractions, of the order of 3° to 4°, are favorable at 
AJ < 4 and increase to about 9° at 1/ = 7 to 10. 

(5) With variable diffusers, optimum contraction in- 
creases with Mach Number and is appreciably larger 
than starting or Wy»ar. contraction but smaller than 
maximum possible contraction. Contraction angles of 
the order of 5° are best at J < 3 and increase to 15‘ 
and more at J = 7 to 10. Minimum pressure ratios, 
appreciably smaller than the normal shock values, are 
obtained at ./ > 2 and amount to only half the normal 
shock stagnation pressure ratio at JJ ~ 7. Although 
performance of variable diffusers deteriorates when 
models are mounted in the working section, it is ad- 
vantageous to use this diffuser type in hypersonic tun- 
nels. 
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Transonic Flow Past Simple Bodies’ 
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ABSTRACT 


Transonic flow past conical tips, a lifting wedge, and a lifting 
flat plate is investigated theoretically and experimentally 
Results on the drag and lift forces are presented. The existence 
of a closed smooth subsonic zone embedded in a supersonic 
flow field is demonstrated experimentally for transonic flow 
past conical tips. The cases discussed, together with the flow 
past wedge sections previously investigated, represent simple 
but fundamental cases for which the flow pattern from subsonic 
through sonic to supersonic flow can be obtained with relative 
ease. An understanding of these cases is extremely useful for 
handling the general transonic-flow problem. 


(1) INTRODUCTION 


, THIS PAPER THE RESULTS of some experiments on 
transonic flow past two simple shapes—namely, a 
cone-cylinder and a double-wedge airfoil—are presented. 
The main aim of these experiments and the small 
amount of theory given here is to bring out some of the 
general features of transonic flow, especially over the 
whole transonic range of Mach Numbers. The authors 
wish to thank Prof. H. W. Liepmann for many dis- 
cussions. 

The special problems studied are the finite cone at 
zero angle of incidence and the symmetric double-wedge 
airfoil at extremely small angles of attack. The experi- 
ments were carried out in the 4- by 10-in. transonic 
Tunnel at GALCIT, which is described in reference 1. 
Details of the experimental technique are given in 
references 2 and 3. The main problem studied here 
for the cone is the breakdown of conical flow, or Taylor- 
Maccoll theory,' and the details of shock detachment. 
The behavior of the head drag curve in passing from 
subsonic to supersonic speeds is also discussed. For 
the double-wedge airfoil, the main problem studied is 
the behavior of the lift-curve slope and center of lift as 
the Mach Number varies. All of the results are pre 
sented in Figs. 1-17, which will now be discussed. 


2) Flow Past CONE-CYLINDERS 
Consider first the cone at zero angle. Figs. 1 and 2 
are schlieren photographs showing attached shock 
waves on 25° semiangle cones at M/.. = 1.415 and 
1.328, respectively. In Fig. 1 the shock is straight 

Presented at the Aerodynamics Session, Twenty-First Annual 
Meeting, IAS, New York, January 26-29, 1953 

* Theoretical research supported by U.S. Air Force contract 
AF 18(600)383. Experimental research partially sponsored by 
the NACA, contract NAw-6154 

t Assistant Professor of Aeronautics and Applied Mechanics 


{ Research Assistant 


near the tip but curved at some distance away. In 
Fig. 2, at a lower Mach Number, the shock is curved 
right from the nose. Now, a curved shock indicates 
that the flow field is not conical and that the corner 
influences the nose. This implies the existence of a 
subsonic region extending from the corner to the nose. 
The location of the sonic line is easily found experimen- 
tally by reflecting weak waves from it as shown in Fig 
3. The sonic lines found that way are shown in Figs. 
t and 5. Thus, for a given cone, as A/., decreases 
conical flow breaks down as soon as the Mach Number 
on the surface reaches 1, as would be expected. It can 
also be noticed on Figs. 4 and 5 that the sonic lines run 
from the shoulder to the nose, that they are straight 
near the nose, and that their location agrees well with 
the conical theory. From this the possibility that the 
flow is locally conical near the tip can be inferred. 

For a more detailed investigation the entire density 
field was found by means of an interferometer. Also 





25° SEMI-ANGLE CONE-CYLINDER 
Me? 1.415 


FIG. | 
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pressure measurements were obtained on the surface 
of various cones by means of models with pressure taps. 
A typical axially symmetric interferogram (30° cone, 
J. = 1.280) is given in Fig. 6. Interferometric re- 
sults presented in Figs. 7, 8, and 9 show the Mach 
Number fields for a 25° semiangle cone for a series of 
three decreasing Mach Numbers 1.401, 1.273, and 
1.229. The theoretical Mach Number at which the 
flow on the surface becomes sonic for this cone is 1.500. 
Already at 1/.. = 1.401 there is a considerable sub- 
sonic region. The sonic line runs from the corner to the 
tip, and both the sonic line and the shock are straight 
near the nose. The pressures measured near the nose 
(Fig. 12) check the Taylor-Maccoll theory, thus veri- 
fying the idea of locally conical flow near the tip. As 
the Mach Number decreases to 1.273; Fig. 8 shows how 
the sonic region is enlarged and how the sonic line now 
runs from the corner to the shock. The shock is 
curved from the tip. Fig. 9 shows a further decrease 
in ./., and the pattern when the wave is detached. 
Some explanation of the various régimes observed 
is given in Fig. 10. In Fig. 10 a transonic form of the 
shock polar is shown, with abscissa u = component of 
velocity perturbation from sonic and v = radial com- 
ponent of velocity. The points on the shock polar 
represent all possible downstream states that can be 
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reached through a shock with the initial state u = u., 
v= 0. 
given wave angle. Thus the points immediately 


Corresponding to each point on the polar is a 


behind the shock wave S on the cone appear as a point 


or points on the shock polar. The curves I to IV are 


integral curves of the equations of motion and are 


1953 


SEPTEMBER, 


found as follows (see reference 4 or 5): It is assumed 
that the flow field is conical so that the physical proper- 
ties p, u, v, etc., are constant on rays, like R in Fig. 10, 
The partial differential equations of motion can be 
expressed as a second-order ordinary (nonlinear) differ- 
ential equation for v as a function of u. The solutions 
of this differential equation are the integral curves. 
Two boundary conditions for the integral curve are 
given on the shock polar—namely, v and dv/du. The 
slope of the integral curve at a point, dv du, determines 
the slope of the corresponding ray in the physical 
plane. An end point C for the integral curve is deter- 
mined by the requirement that the slope of the ray, 
depending on dv/du, should equal v, a* so that, approxi- 
mately, a solid cone surface can be inserted in the flow. 
The locus of all the possible end points of integral 
curves starting on the shock polar is called the ‘‘apple”’ 
curve. A typical integral curve corresponding to a 
purely supersonic field is represented by I in Fig. 10. 
The point S corresponds to the state immediately be- 
hind the shock. Each point R on the integral curve 
corresponds to the entire ray R in the physical plane. 
The end point C on the ‘apple’ curve represents the 
ray that lies along the cone surface. If now we think 
of increasing the cone angle, with ./., fixed, the integral 
curve moves to the left in Fig. 10 until a curve like II 
At II the end point is sonic so that sonic 
If the cone 


is obtained. 
velocity is reached on the cone surface. 
angle is increased further, subsonic flow occurs next to 
the surface, the corner influences the nose, and conical 
flow in the large is no longer possible. However, the 
flow may still be locally conical, as indicated by the 
integral curve III. As the cone angle is increased from 
II, the velocity at the corner remains sonic. The 
velocity near the tip is given by the end point of III, 
and the cone surface is represented py the horizontal 
line from the end point of III to the sonic line. The 
entire integral curve maps to the nose point in the phys- 
ical plane; each point on the integral curve corresponds 
to a ray in the neighborhood of the nose. The sonic 
line intersects II] or runs from the corner to the nose. 
Now as the cone angle is increased further, integral 
curves like IV are obtained which are the same as II], 
except that the sonic line runs from the nose to the 
shock. As the cone angle is increased still further, the 
cone surface cannot be reached as a point on the ‘‘apple”’ 
curve, or by a locally conical flow, so that the shock 
must be detached. Thus the detachment is given by 
the largest possible deflection that can be obtained. by 
conical flow or is correctly given by Taylor-Maccoll 
theory. It should be remarked that an experimental 
determination of the detachment Mach Number is 
difficult because of the slow rate of detachment for a 
cone. Thus, we have a description in the hodograph 
plane of all the régimes observed experimentally. 

The effect of these various régimes as .\/.. decreases, 
for a fixed cone, is shown in the head-drag curves, 


based on p., in Fig. 11. For high values of 7. there 
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is good agreement with conical-flow theory. As A. 
decreases, conical-flow theory breaks down when the 
surface Mach Number is |. As ./,, is decreased fur- 
ther, the pressure is increased only near the tip and 
Thus, the drag curve 
Subsonic points are 


remains sonic near the shoulder. 
falls below the conical theory. 
also given on the figure and show a gradual rise of drag 
as J, ~ 1. The slope of the head-drag curve at 
M.. = 1 can be given approximately by applying the 
Principle of Stationarity of Local Mach Number. This 
Principle was enunciated and verified for two-dimen- 
sional flows in reference 7. The Principle states that 
(OM OM.)u..-1 = 0; it is based on the idea that, as 
M.. — 1, the head shock is infinitely distant from the 
For \/..2 = 1 + ¢, the Mach Number behind 
Then the flow 


body. 
the normal shock is \/,? = 1 — e. 
field at large distances is effectively subsonic with a 
Mach Number \/,.2 = 1 — «. Thus the local Mach 
Number remains the same as the free-stream Mach 
Number passes through |. By applying the Principle 
and using the pressure-Mach Number relations, it can 


Co Mao=!l 
2 


Roughly estimat- 


be deduced that 


(26) _ 4 ( 
ol, Mo =1 er 


where C, = head-drag coefficient. 
1 gives a good estimate for the slope as 
The stationarity of the local Mach 
The re- 
Mach 


ing Cpy, 
shown in Fig. 11. 
Number was also checked experimentally. 
12, where the local 


sults are shown in Fig. 


Number is plotted as a function of the free-stream Mach 
Number for various chordwise stations. The flattening 
of the local Mach Number curves as J/.. — | is evident. 

The Principle of Stationarity can be extended to the 
curvature of local Mach Number (07.)//041/.*)y7.-1 = 
0. This result follows from the plausible assumption 
that the local Mach Number is a nondecreasing func- 
tion of the free-stream Mach Number near J/,. = 1 or 
O0M/OM,, 2 0. This assumption is verified experi- 
mentally in reference 7 and in Fig. 12.7 

It should also be remarked that Yoshihara has com- 
puted, by relaxation calculations, the pressure on a 
slender cone at JJ, = 1.00 and has experimentally 
verified the results.* 

The above experimental results cannot be used to 
test the transonic axially symmetric similarity rules as 
derived by von Karman®* and Oswatitsch and Berndt." 
The derivation of the similarity rules employs some 
linearized surface tangency condition that is probably 
valid for small cone angles. The nose angles of the 
experimental models were too large for the lineariza- 
tion of the boundary conditions to be valid. 

It can be shown, however, that the axially symmetric 
transonic equations are a good approximation to the 
exact equations. A second-order solution for the flow 
past a semi-infinite cone has been found* using the tran- 
sonic equations, and the results for local Mach Number 
agree well with the exact conical solution (see Fig. 14). 

+ The same results can be deduced analytically (for example, 


see reference 13). 
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Fig. 14, a schlieren picture of the typical flow pattern 
of acone at \/.. < 1, is the last figure pertaining to the 
cone. It is noteworthy that there is no shock of appre- 


ciable strength emanating from the corner. 


(3) DoUBLE-WEDGE AIRFOIL 


Surface pressure measurements were made for the 
completely symmetric 5 per cent thick double-wedge 
airfoil at small angles of attack by sliding a static pres- 
sure probe along a line touching the wedge shoulder 
and parallel to the wedge chord line. The details of 
the technique are given in reference 2. The lift coeffi- 
cients are obtained by integration of the pressure dis- 
tributions, and the results are plotted in transonic 
similarity form in Fig. 15. The lift-curve slope C,/a 
for small angles is given, and it should be remarked that 
over the range of angles tested, a = 0° to 2°, the lift 
was a linear function of a forall 1... The salient fea- 
tures of Fig. 15 are as follows: For low .\/., the lift- 
curve slope follows the Prandtl-Glauert correction, but 
the incompressible slope for this airfoil is somewhat less 


than 27. There is a striking drop in lift coefticient 


as ./., increases near 
t. = (M.2 — 1)/[(y + 1) (/o]” = - 1 


As M.. — 1, the lift-curve slope goes toward the theo- 
retical value given by Guderley and Yoshihara in refer- 
ence 11. An explanation of the drop in lift is the follow- 
ing: As .V/., increases, a supersonic zone forms behind 
Near —. = 


grows rapidly in size, and its terminating shock moves 


the shoulder. —1l, this supersonic zone 
from the corner to the trailing edge with a small in- 
crease in &, The rear half of the airfoil operating 
supersonically is inefficient in the production of lift. 
It should be emphasized that the phenomenon is not 
one of separation or stalling, since the boundary layer 
remained attached during the experiments. 

The arrival of the terminating shock at the trailing 
edge is interesting from the point of view of theory. 
It means that the conditions for the flow leaving the 
The problem 
This 


trailing edge can be satisfied locally. 
of satisfying the Kutta condition disappears. 
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fact was used by Guderley and Yoshihara in their solu- 
tion at /.,, = 1." 

Pressure distributions showing the three typical 
cases—practically no shock, shock partway back, and 
shock off the trailing edge—are shown in Fig. 16. It 
might be remarked that it is the effect of the boundary 
layer which prevents any supersonic zone from forming 
at M.. = 0.309. 

Another interesting quantity is the location of the 
center of lift. If the local pressure difference across 
the airfoil is linear in the angle of attack, so are the 
liftand moment. The center of lift is thus independent 
of the angle of attack. 
moves rearward to a maximum value and then forward 
At low subsonic speeds, the center of lift 


However, the center of lift 


as M, > 1. 
is located at one-quarter chord; for 1/.. = 1, the center 
of lift is located at 30 per cent chord, according to the 
solution of Guderley and Yoshihara.'' Fig. 17 shows 
the experimental result illustrating the rearward trend 
of center of lift as .\7.. increases. The maximum rear- 
ward position of center of lift (40 per cent chord) corre- 
sponds to .\/.. such that the terminating shock waves of 
the supersonic zone leave the airfoil. The center of 
lift travel can be explained as follows: As .\/.. increases 
from zero, supersonic zones terminated by shocks ap- 
pear behind the corner. Because of the small angle of 
attack, the terminating shocks occur at different chord- 
wise stations (Fig. 16) resulting in an appreciable lift 
over the rear half. This shifts the center of lift toward 
the rear. The effect increases with increasing free- 
stream Mach Number until the shocks leave the trailing 
edge. Any further increase in Mach Number produces 
relatively little change in lift on the rear half, and the 
center of lift moves forward. Theoretically, at J/.. = 
1, only 25 per cent of the lift comes from the rear half."' 
The theoretical center of lift for 1/.. > 1 is a continua- 
tion of the subsonic curve. The center of lift moves 
forward as ./.. increases until a Mach Number close 
to the attachment Mach Number is reached. The 
center of lift then moves rearward to the one-half chord 
point, its supersonic value. The theory for ./.. > 1 
is presented by Vincenti and Wagoner.'” 
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A nonlinear effect also appears for the higher Mach 
Numbers. When the terminating shock waves are no 
longer on the airfoil, the center of lift travels forward 
as the angle of attack increases. This means that the 
moment is not linear in a. The forward shift comes 
from the strong expansion evident on the upper surface 
0.921). It is 


interesting that there are deviations from linear behay- 


near the leading edge (Fig. 16, |7. = 


ior at angles of attack which are less than the wedge 
nose semiangle (3°). 
Two problems are interesting for future experi- 


ments—namely (a) evaluating the range of validity 
of the assumption of linearity in a@ (theoretically a flat 
plate has C, ~ a‘ at MJ. = 1); (b) behavior of lift 


near the wing tips. 
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Free-Flight Determinations of the Drag 


Coefficients 


of Spheres 


ALBERT MAY* ano W. R. WITT, JR.? 
U.S. Naval Ordnance Laboratory 


SUMMARY 


Drag coefficients have been determined in the NOL Pressur 
ized Ballistics Range for steel and aluminum spheres '/, to 4/4 
in.in diameter. The Mach Number ranged from 0.8 to 4.7, and, 
by using air pressures from 1/120 to 1 atmosphere, the Reynolds 
Number was varied separately over a range of 1.14 K 10° to8.4 X 
108 The drag-coefficient data obtained, together with the 
results of other investigators, permit a plotting of the variation of 
drag coefficient over a large region of Mach Numbers and Reyn 
The variation of Cp was found to be only about 


103 < Re < 


olds Numbers. 
10 per cent in the range 1.6 < M < 4.7 and 4.0 X 


1.0 X 108 


INTRODUCTION 


COEFFICIENTS of 


oe TIONS OF THE DRAG 
spheres at various Mach and Reynolds Numbers 
have been made in the NOL Pressurized Ballistics 


Range for the Engineering Research Institute at the 
University of Michigan in connection with upper- 
atmosphere research for the Meteorological Branch of 
the Signal Corps. Because of the low air densities in 
the upper atmosphere, the Reynolds Number is much 
lower for a given Mach Number and sphere size than 
it would be at sea level. The purpose of the present 
program was to provide data on the drag coefficient of 
spheres in a supersonic region (1.4 < AJ < 3.3, 10° < 
Re < 10°) in which reliable data were almost entirely 
lacking. 

The Pressurized Ballistics Range, shown in Fig. 1, 
is an enclosed tube 318 ft. long permitting the firing of 
projectiles through air at any pressure from 6 atmos- 
pheres to less than 1/100 atmosphere. The air tem- 
perature of the room in which the range is located is 
maintained at 74°F., and the temperature in the range 
is changed only a few degrees from this value by varia- 
tions of the range pressure. Photographs can be taken 
by the spark-shadowgraph method at any number of 
the 25 stations with which the range is equipped. The 
positions of these stations are known to an accuracy 
better than 0.002 ft., and readings of the shadowgraph 
plates can be made to the same accuracy. The relative 
times at which eight of these shadowgraphs are taken 
can be determined by the use of an RCA counter chrono- 
graph. The counters read to 10~ sec., and inaccura- 
Received March 4, 1953 
* Chief, Range Branch, Aerodynamics Division. 


' Project Engineer 


the time determinations are less than 0.2 


microsec. 


cies in 


EXPERIMENTAL DETAILS 


The spheres were fired from a smooth-bore gun having 
a caliber of 20 mm., with the sphere niounted in a light 
plastic sabot as shown in Fig. 2. After the sphere and 
sabot are ejected from the gun, the gas pressure that 
has built up between sphere and sabot causes them to 
separate, and the sabot is left behind because of its 
large deceleration. At very low air pressures, the drag 
on the sabot was inadequate to effect proper separation 
of sabot and projectile. This difficulty was obviated 
by using a sabot split lengthwise so that the halves of 
the sabot were blown apart rapidly. 

In order to obtain drag-coefficient values for a Reyn- 
olds Number range of 100 to 1 while the Mach Number 
varied only by a factor of 2.4, both sphere size and air 
The Reynolds Number, pvd/ xn, 


pressure were changed. 





Fic. 1. NOL Pressurized Ballistics Range permitting independ 
ent variation of Mach Number and Reynolds Numbet 





Fic. 2. Several sizes of spheres with their launching sabots. 
Left, half a sabot split to hasten separation from the sphere 
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TABLE | 


Drag Coefficients from the Firing of Spheres 


Pressure, Reynolds 
Diam- Mm. per Mach Number 
Round eter, In. Material Hg Number X 10~4 Cp 
362D 3/4 Steel $41 1.90 47.1 0.999 
362C 3/4 Steel 290 1.94 31.6 1.000 
361D 3/4 Steel 271 2.62 10.0 0.972 
364D 3/4 Steel 184 2.66 27.5 0.950 
363D 3/4 Steel 184 L.78 18.32 1.012 
358D 3/4 Steel 174 o:T2 27.0 0.947 
360D 3/4 Steel 127 2.62 18.61 0.947 
363C 3/4 Steel 124.5 1.60 tl... 1.019 
364C 3/4 Steel 109.2 2.71 16.53 0.941 
438C 3/4 Steel 62.3 2.85 9.79 0.9382 
441B 1/2 Steel 648.0 3.48 84.5 0.964 
440C 1/2 Steel 311.0 +. 20 18.9 0.934 
361C 1/2 Steel 185.0 1.44 9.96 1.012 
439B 1/2 Steel 124.1 1.25 5.80 0.955 
441A 1/2 Steel 120.5 4.69 21.0 0.906 
358C 1/2 Steel 105.0 2.15 8.40 0.942 
438D 1/2 Steel 101.8 2.96 11.28 0.955 
360C 1/2 Steel 76.5 2.62 7.45 0.937 
360A 1/2 Steel 75.5 1.66 1.67 0.942 
437D 1/2 Aluminum 76.6 1.52 4.32 0.923 
437A 1/2 Aluminum 34.4 2.57 3.31 0.923 
436A 1/2 Aluminum 13.0 3.01 1.518 0.938 
491C 1/2 Aluminum 6.85 2.94 0.810 0.978 
490B 1/2 Aluminum 6.84 2.45 0.671 0.985 
489B 1/2 Aluminum 6.65 3.26 0.870 0.957 
490C 1/2 Aluminum 6.55 2.70 0.712 0.986 
491B 1/2 Aluminum 6.32 2.19 0.559 0.995 
442B 1/2 Magnesium 24.4 1.28 1.161 0.850 
438B 1/2 Magnesium 21.32 1.46 1.155 0. 884 
440D 3/8 Steel 761.4 3.13 73.5 0.962 
435 3/8 Steel 81.6 3.01 6.88 0.934 
443C 5/16 Aluminum 88.0 4.58 9.38 0.906 
367B 9/32 Steel 143.3 1.66 5.04 0.948 
367A 9/32 Steel 103.9 1.79 3.89 0.9438 
362A 9/32 Steel 36.5 2.83 2.15 0.938 
366B 9/32 Steel 35.0 2.67 1.942 0.940 
266A 9/32 Steel 23.0 2.80 1.333 0.951 
365A 9/32 Steel 16.8 2.53 0.926 0.895 
361B 1/4 Steel 70.0 1.46 3.20 0.916 
439A 1/4 Bronze 181.8 3.71 9.23 0.940 
437B 1/4 Aluminum 30.2 3.71 1.514 0.955 
442A 1/4 Aluminum 9.01 1.66 0.449 0.986 
443A 1/4 Aluminum 7.60 1.17 0.166 0.956 
444A 1/4 Aluminum 7.60 1.07 0.142 0.789 
434C 1/4 Aluminum 7.40 1.80 0. 266 1.022 
432 1/4 Aluminum 7.40 1.72 0.262 0.980 
444B 1/4 Aluminum 7.00 0.82 0.114 0.632 
491A 1/4 Aluminum 6.75 2.27 0.308 1.033 
489C 1/4 Aluminum 6.70 3.00 0.405 1.003 
490A 1/4 2.38 0.315 1.022 


Aluminum 6.70 


is proportional to the air density p and the sphere 
diameter d. The values of the coefficient of viscosity u 
and of the speed of sound depend little on the air pres- 
sure, and their change with temperature is small in the 
firing range because of the approximate constancy ot 
temperature. Hence, the Mach Number, ./, is roughly 
proportional to the velocity, v, and the Reynolds 
Number is approximately proportional to the product 
M pd. 

The lowest values of pressure used (about 1/120 
atmosphere) were near the minimum easily attainable 
in the Pressurized Ballistics Range. Because of the 
small value of the drag force at such pressures, alumi- 
num spheres were used to increase the deceleration. 
On the other hand, elevated pressures were avoided 
since above 1 atmosphere, even with steel spheres, the 
change of speed becomes so great that a simple de- 
pendence of Cp on Mach Number cannot be assumed 
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Spark shadowgraph of a !/2-in. steel sphere traveling at 
M = 3.48, Re = 8.45 X 10°, in air at 0.85 atmosphere pressure 


Fic. 3. 


throughout the speed range. The sphere diameter was 
varied between !/ 4 in. and * 4 in. only. 

Fig, 3 is a print obtained from a typical spark shadow- 
graph (round 441B; see Table 1) of this program. This 
round consisted of a '/»-in. steel sphere fired at \/ = 
3.48 with an air pressure of 0.85 atmosphere. At 
lower pressures, the shock waves become less _pro- 
nounced in the photograph, and at the lowest pres- 
sures used nothing is seen in the shadowgraphs but the 
sphere itself. Measurements for drag determinations 
are made at the rear of the sphere to avoid the optical 
distortion in front of it. 

Spark shadowgraphs were taken only at the eight 
timing stations. Position-time data for determina- 
tion of the drag coefficient were obtained by measure- 
ments of the sphere position (relative to a fiducial 
mark) on the photographic plates and from readings of 
the timing counters. From these data, values of Cp 
were calculated by a least-square technique. 


PREVIOUS DRAG DETERMINATIONS 


On the Mach Number-Reynolds Number coordinates 
of Fig. 4, the shaded areas indicate those regions in 
which reasonably reliable drag data for spheres were 
already available. Values of Cp for low Mach Num- 
bers had been compiled by various authors':* for a 
wide range of Reynolds Numbers terminating just 
above the critical value. These are indicated by the 
shaded area immediately adjacent to the abscissa in Fig. 
4. 

The area that branches off from the one already 
described (near Re = 4 X 10*) was investigated by 
Charters and Thomas in the free-flight range at Aber- 
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Fic. 4. Contour plot of Cp of spheres 


Re 


The drag coefficient is constant along each line 


10° 10° 


Dashed lines are approximate only 


1o* 


Shaded areas are regions where reliable data were already available, and the dots show by their positions the Mach Numbers and 
Reynolds Numbers of the determinations made in the NOL Pressurized Ballistics Range 


deen, and their research is the subject of a compre- 
hensive report.* 

Other free-flight Mach 
Numbers and Reynolds Numbers were carried out at the 
New Mexico School of Mines‘ and in England.’ Zahm‘ 
investigated a small transonic region at high Reynolds 
Numbers by free-flight tests. The large supersonic 
region indicated at small Reynolds Numbers was in- 
vestigated in a wind tunnel at the University of Cali- 
fornia by Sherman and Kane’ and by Jensen.* 


determinations at higher 


EXPERIMENTAL RESULTS 


In the present program, 50 rounds were fired to 
provide Cp data. The Mach Numbers and Reynolds 
Numbers of these shots are indicated by the positions 
of the dots on Fig. 4+. The experimental data for the 
rounds are given in Table 1. 

The three-parameter data for these rounds (Cp, .\/, 
and Re) are represented in Fig. 4 by a contour plot 
using lines of constant Cp. These lines were deter- 
mined by inserting on the ./-Re coordinate axes the 
Cp values found in previous investigations and the 50 
values found in the present study and by estimating the 
best lines of constant Cp. In some areas, too little 
data were available to permit dependable estimates, 
and the contours in these areas have been dashed to 
indicate that they are qualitative only. It will be seen 
that the variation of Cp is only about 10 per cent 
throughout the region in which 1.6 < MW < 4.7 and 


4.0 X 10°< Re< 1.0 X 10°. 


There appears to be unmistakable evidence of a 
saddle point in the vicinity of 17 = 2.0, Re = 3 X 
10%. In passing along the line Re = 3 X 10', a maxi- 
mum is found near this point, while a minimum is 
encountered on the line 1/ = 2.0. 

No evidence was found of the existence of a critical 
Reynolds Number in that portion of the supersonic 
region investigated. Charters reported*® that he found 
no critical Reynolds Number in a portion of this region. 
The curves of Fig. 4 indicate that over the range of 
Mach Number and Reynolds Number in that figure, 
any critical Reynolds Numbers are restricted to \/ < 
0.7. 

The drag data of Fig. 4 are presented in more familiar 
Fig. 5 gives graphs of Cp against 
Re for various Mach Numbers. The lowest graph is the 
well-known plot! * for the incompressible case. Lack 
of data still permits only a guess as to the behavior 
At Reynolds Numbers 


form in Figs. 5 and 6. 


immediately above this region. 
between 10° and 10°, it is clear that Cp decreases as the 
Mach Number increases above about 1.6. From Fig. 4 
it appears probable that such a decrease of Cp with 
increasing Mach Number may be expected at higher 
Mach Numbers throughout the Reynolds Number 
range from 10 to 10°. 

Fig. 6 gives graphs of Cp against Mach Number for 
several values of Reynolds Number. The similar graph 
given by Charters and Thomas* was obtained by firing 
spheres at various speeds using fixed air pressure and a 


diameter range of 5 to 1. Consequently, their plot is 
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Fic. 6. Variation of Cp of spheres with Mach Number for 


several Reynolds Numbers. Dashed curves indicate that their 
positions are only approximate. The circles show values taken 
from a curve given by Charters and Thomas.* 


an average over Reynolds Numbers varying by 5 to | 
for a given Mach Number. Values taken from Char- 
ters’ curve are shown by the small circles superimposed 
on Fig. 6. 

The accuracy of the drag-coefficient values is difficult 
to determine. Even the known sources of error vary 
greatly throughout the pressure and velocity ranges 
investigated and are difficult to appraise. At very low 
air pressures (6 to 7 mm.) accuracy was probably 
poorest because of errors in pressure measurement. 
For example, from a study of the consistency of readings 
made during calibration of the gage used for round 491B 


(pressure = 6.32 mm.), it appears that the drag co- 


AERONAUTICAL 


The dashed curve is only approximate. 


SCIENCES—SEPTEMBER, 1953 
efficient for that round may be in error by as much as 
2 percent. For the */4-in. spheres, known inaccuracies 
would not cause errors larger than 0.6 to 0.9 per 
cent. 

The accuracy of the individual drag values is not 
easily estimated from a plot such as Fig. 4 since the 
placing of the contours depends on these values. Of the 
50 values of Cp found, 36 appear to be in error by | 
per cent or less. Of the others, seven seem to have 
errors between | and 2 per cent, and seven to have 
errors greater than 2 per cent. In only one or two 
cases does it appear sure that the error is larger than 3 
per cent. 

Twenty values in the rather ‘‘flat’’ region between 
0.03 and 0.95 contours (in general not shown) were 
The mean drag coefficient was 0.942, and 
It is thought that 


averaged. 
the standard deviation was 0.010. 
all contours shown by solid lines on Fig. 4 are accurate 
to 1 per cent and, in most regions, probably to !, 
per cent. 

An attempt to correlate the angle at which the flow 
separates from the sphere with drag coefficient, as the 
Reynolds Number changes, was unsuccessful. This is 
not surprising in view of the small magnitude of the 
drag coefficient variations. 
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Bending Tests of Thin-Walled Sandwich 


Cylinders’ 


GEORGE GERARD? 
New York University 


SUMMARY 


A series of tests was conducted on thin-walled circular sand 
wich cylinders under bending loads using aluminum-alloy faces 
and either cellular cellulose acetate or end grain balsa cores 
It was found that for cylinders with cellular cellulose acetate 
cores, the buckling strength averaged 32 per cent higher than the 
theoretical value. From previous tests under compressive loads, 
excellent agreement was obtained between the corresponding 


theoretical value and tests. Thus, it was concluded that, under 


bending loads, the buckling stress depends upon the average stress 
on the cross section rather than on the maximum stress. This 
behavior has been noted previously from tests on homogeneous 


cylinders 
INTRODUCTION 


AS PART OF A PROGRAM involving a theoretical and 
experimental study of the buckling behavior of 
sandwich cylinders, this paper is concerned with buck 
ling tests of thin-walled sandwich cylinders under bend- 
ing loads. Tests were conducted on cylinders con- 
structed with cores weak in shear (cellular cellulose 
acetate) and cores strong in shear (end grain balsa) to 
check the validity of a method for predicting overall 
buckling of sandwich cylinders under bending loads. 


TEST PROGRAM 


Test Specimens 

The large number of variables which influences the 
behavior of a composite element such as sandwich con- 
struction requires that in a test program attention be 
directed toward the significant parameters only. Thus, 
in order to study the overall buckling characteristics 
of sandwich cylinders, it is necessary that the sandwich 
thickness be as small as possible. The 1/s-in. thick 
sandwich used was considered to be the minimum thick- 
ness feasible. The nominal radius of 6 in. which was 
used for all cylinders was dictated by the capacity of 
the testing facilities. To ensure uniformity of intro- 
duction of bending load, the cylinders tested were 36 
in. long. 

The significant parameter in overall buckling of 
sandwich elements is the shear rigidity of the core ma- 


terial. Consequently, two different materials were 
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neering. 


selected for the core: cellular cellulose acetate (CCA) 
of approximately 6 Ibs. per cu.ft. density and end grain 
balsa (EGB) of approximately 5 to 5.5 Ibs. per cu.ft. 
density. For all cylinders, the face material used was 
0.010-in. 24S-T3 aluminum alloy. 

The cylinders used in this investigation were manu- 
factured by Skydyne, Inc., of Port Jervis, N.Y. The 
fabrication details were briefly as follows. The alumi 
num-alloy faces were coated with a metal primer prior 
to assembly to the core. The core was planed to the 
proper thickness, and a secondary adhesive was used 
to bond the faces to the core. Although the adhesive 
used will polymerize at room temperature, sufficient 
heat was applied to accelerate setting of the adhesive. 
Pressure necessary to bring the surfaces in intimate con- 
tact was applied, and the assembly was molded for a 
period of 2 to 3 hours. 

To facilitate installation of the test specimen in the 
loading jig, 2-in. wide hardwood blocks were bonded to 
both the inside and outside faces at each end of the 
cylinder (see Figs. | and 2). The end faces of these 
blocks were then turned down on a lathe to ensure 
planeness of the face, perpendicularity with the axis of 


the cylinder, and parallelism of the faces. 


Test Procedure 


The pertinent dimensions of each specimen were de- 
termined as follows. Four evenly spaced longitudinal 
stations were lightly scribed on the outside wall of the 


The outside radius was then measured to the 


cylinder. 
nearest 0.01 in. at each pair of stations at both ends. 
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TABLE | 
Average Dimensions and Failing Loads of Sandwich Cylinders 
Maximum Maximum 
Average Failing Failing Failing 
Core Mean Total Face Moment Load per In Stress 
Cylinder Material Length, In Radius, In Thickness, In Thickness, In In. Lbs Lbs. per Sq.In. Lbs. per Sq In 
| 8 R (h + t) t Wf, \ U/arR a 
I CCA 32.09 6.48 0.141 0.0105 67,900 514 14.300 
< CCA 32.13 6.50 0.125 0.0105 46,100 347 16,500 
3 CCA 32.13 6.49 0.125 0.0105 77.700 587 28.000 
4 CCA 32.16 6.53 0.141 0.0105 89 600 670 31.900 
5 EGB 31.29 6.47 0.141 0.0105 53,100 404 19,300 
6 EGB 31.09 6.55 0.156 0.0105 65,700 500 25.800 


The length of the cylinder between hardwood blocks 
was also determined at each station to the nearest 0.01 
in. 

The total thickness of the cylinder wall was measured 
at each station at both ends just inside the hardwood 
blocks and at the center by means of a special dial gage 
with a least division of 0.001 in. The thickness of the 
faces was measured with a ten-thousandths micrometer 
on samples of the face material. The average dimen- 
sions for each cylinder are given in Table 1. 

Installation of the specimens in the specially designed 
bending loading jig was accomplished by bolting the 
specimen to the steel end plates through the hardwood 
blocks at the ends of the cylinder. The loading jig 
is shown in Fig. | with a specimen installed and tested. 
Load was applied by means of a hydraulic tension jack 
and was measured by a calibrated strain-gage tension 
kink at the jack which was accurate to +2 per cent. 


Physical Property Tests 


The important physical properties in overall bending 
buckling of sandwich construction are the compressive 


material and 
The shear 


stress-strain characteristics of the face 


the shear modulus of the core material. 
moduli of the core materials were determined by testing 
as simple beams, strips of flat sandwich construction. 
The beams were fabricated by the manufacturer of the 
cylinders of the same materials and in the same manner 
as the cylinders themselves. The test consisted of ob- 
taining load-deflection data from which the shear modu- 
lus of the core material was determined from the dif- 


ference in deflection between the experimental value 


TABLE 2 
Data on Instability Patterns 
Longi Circum 


tudinal ferential 


Wave Wave 
Type of Length, Length 
Cylinder Instability Location In., Lr In., Le 
-1 '/sin. CCA Buckle Central 1 3 
2'/sin. CCA Buckle Central 1 3 
3 '/sin. CCA Buckle Central I'/4 3'/2 
4'/sin. CCA Antisymmetrical Both ends 7/16 rR 
wrinkle 
5 '/sin. EGB Buckle Central 1'/4 3! 
6 '/sin. EGB Antisymmetrical At end 7/14 rR 


wrinkle 


and that computed for the beam neglecting shear de- 


flections. From at least four tests of each core material, 


the following values of shear modulus were ob- 


tained: 
cellular cellulose acetate; G = 3,400 Ibs. per sq.in. + 
17 per cent 


st 


end grain balsa; G’ = 17,600 Ibs. per sq.in. + 


15 per cent 


The with-grain compressive stress-strain characteris 
tics of the 0.010-in. 24S-T3 aluminum alloy used for 
the faces of the sandwich cylinders was determined 
by use of a solid guide compression jig of the National 
Bureau of Standards type. The stress-strain data indi- 
cated that the proportional limit of this material was 
54,000 Ibs. per sq.in. Thus, the stress in the faces at 
failure of the sandwich cylinders (Table 1) was within 


the elastic range in all cases. 





Fic 2. Overall buckling of sandwich cylinder in bending, 








Test | 

Av 
ders | 
A ph 
the « 
tiona 
givel 
obtai 


Cor! 
TI 
unde 
whet 
anal’ 
a 
Tl 
there 
refer 
h, of 
aver 
cent 
grou 
acco 
load 
was 
the 
thar 
Ir 
whic 
theo 
ing | 
rigic 
7 pe 
lbs. 
was 
peri 
evi 
50 | 
app 
valu 
T 
thar 
sam 
gen 
that 
non 
stre 
stre 
stre 
' 
refe 
ZC Vi 


im 





0 
0 
0) 
() 


r de- 
‘rial, 
ob- 


TIS- 
for 
ned 
mal 
idi- 
was 
; at 
hin 

















TESTS OF 


Test Results 


Average measured dimensions of the sandwich cylin- 
ders together with the failing loads are given in Table 1. 
A photograph of a typical buckle pattern taken after 
the conclusion of the test is shown in Fig. 2. Addi- 
tional data on the buckling and wrinkling patterns are 
given in Table 2. The measured wave lengths were 
obtained in regions of relatively little distortion. 


CORRELATION OF THEORY WITH EXPERIMENTAL DATA 


Theoretical work on buckling of a sandwich cylinder 
under bending loads is given in reference | for the case 
where the core is weak in shear. The results of this 
analysis indicate that for this case, the buckling load, 
N.,, is equal to the shear rigidity of the core h G. 

The | s-in. CCA cylinders fall into this category, and, 
therefore, it is possible to check the theoretical work of 
reference | directly. The average centroidal height, 
h, of cylinders —1 to —4 is 0.118 + 6 per cent, and the 
average shear modulus is 3,400 Ibs. per sq.in. + 17 per 
cent. Therefore, the average shear rigidity of this 
group of cylinders is 400 Ibs. per in. + 23 per cent, which 
according to reference | is the theoretical buckling 
The average experimental failing loading 
Thus, 


loading. 
was found to be 530 Ibs. per in. + 33 per cent. 
the experimental failing loading is 32 per cent higher 
than the theoretical value, on the average. 

In reference 2, comparable tests were run on sand- 
which cylinders under axial compressive loads. The 
theory of reference 2 indicated that the buckling load- 
ing of cylinders weak in shear is also equal to the shear 
rigidity of the core. For an average h = 0.115 in. + 
7 per cent and an average shear modulus G = 2,000 
lbs. per sq.in. + 20 per cent, the average shear rigidity 
was 230 Ibs. per in. + 27 per cent. The average ex- 
perimental failing load of eleven !/s-in. CCA sandwich 
cylinders given in reference 2 was 235 Ibs. per in. + 
50 per cent. Thus, the experimental failing load was 
approximately 2 per cent higher than the theoretical 
value, on the average, indicating excellent agreement. 

The fact that the bending buckling load is higher 
than the compressive buckling load for cylinders of the 
same dimensions has been observed in tests on homo- 
geneous cylinders in reference 3. It has been suggested 
that this behavior may be due to the buckling phe- 
nomena responding to the average value of compressive 
stress on the cross section rather than the maximum 
stress. In the case of pure bending, the maximum 
stress would be 1.4 times as great as the average value. 

For sandwich cylinders in compression, the theory of 
reference 2 indicates that value of the parameter Et/GR 


governs the behavior of the cylinder. Thus, it has 
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been found that cylinders can be considered weak in 
shear when the value of ET/GR exceeds unity. To ef- 
fect a comparison between the cylinders loaded in com- 
pression and bending which were differently constructed 
as concerns G and R, it appears that, if the values of 
Et/GR were equal or at least greater than unity for 
each group, a direct comparison would be valid. 

The value of this parameter for the compression 
specimens averaged 5, and for the bending specimens 
kt/GR = 5 also. Thus, since the compressive test 
data check theory well, it can be concluded that in 
bending the 32 per cent higher experimental value is 
in agreement with the behavior observed for homogene- 
ous cylinders, thus accounting for the discrepancy. 


Although these conclusions are based on data with 
which large percentage variations are associated, it ap- 
pears that the trend is clearly indicated. The fact 
that the scatter of the experimental values is consider- 
ably higher than that for the core properties is based 
on the fact that the latter values were established from 
beam bending tests. In a beam test, the results are a 
statistical average of the core properties, whereas a 
structure subject to buckling responds to local weak 
spots caused by defects in the core properties. Thus, 
the scatter would tend to be greater for the buckling 
data. 

The '/s-in. end grain balsa cylinders, —5 and —6, 
can be considered to be in the weak-in-shear category, 
since, for this group, Et/GR = 14. These cylinders, 
therefore, failed at a small fraction of the theoretical 
load. This is evidenced by the fact that the average 
failing loading of this group was 452 Ibs. per in. as 
compared to 530 Ibs. per in. for the '/s-in. CCA cylin- 
ders, although the shear modulus of the EGB core was 
over five times as great as the CCA core. The fact 
that the EGB cylinders were not even as strong as the 
CCA cylinders is difficult to reconcile with the core 
properties. Unless some other mode of buckling is 
responsible for this anomalous behavior, it appears that 
poor bonding to the EGB core may be responsible for 
the low values of failing load. 
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where | 
A = 0/p¢ 


Fe THE PRESSURE DISTRIBUTION on a thin rectangular wing 
in supersonic flow at a constant angle of attack a (originally 
given by Busemann), it may be derived immediately that the lift 
dL of a spanwise strip (with span ) and width dx) at the chordwise 


station x is 
dL = 2BpU*a(b — Bx) dx (1) 


where 

B = (M? — 1)" 
In other words, the chordwise lift distribution decreases linearly 
from the Ackeret value, which is found at the leading edge, x = 0 
The formula is valid beyond the intersection of the two wing-tip 
Mach cones, but it is limited to chordwise stations ahead of the 
intersection of the wing-tip Mach cone with the opposite side- 
edge. Thus, the chord ¢ can be at most )/8. 
With the same restriction, the lift for a cambered rectangular 
At a chord 


wing can be given. Let a@ be now a function of x. 


wise station x;, the change in a over a length dx; isda = a’(x;)dx). 
The contribution of this change in @ to the lift at a station x > x, 
can be expressed by taking a rectangular wing with leading edge 
at x; and a constant angle of attack da = a’(x;)dx;: 


d(dL) = 2BpU2a'(x,)dx,[b — B(x — x,)]dx 


Superimposing these contributions, dZ at the station x becomes 
the Duhamel integral computed with this integrand, leading to 
the expression 


Px 
dL = 2Bpl?al0)\(b — Bu)dx + { Mal) 
x 


. 


~ 
= 28BpU*|ba(x) — B f, a(x,) dx,;| dx (2) 
The contribution of the same strip to the drag D is given by 
dD = ad(x)dL (3 
Put 
a(x) = d2y/dx (4) 


Let a = 0 
Introducing the nondimensional 


where 2 is the wing ordinate (positive downward. ) 
at the leading edge x = ‘0 
quantities 


the following expressions are obtained for the lift and drag coef- 


ficients: 
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is the ‘“‘reduced”’ aspect ratio > | 

The present note is concerned with the determination of a 
camber ¢(£), which renders the drag due to lift a minimum Ir 
linearized theory this drag is independent of thickness drag 
Because of the assumption of no twist, the minimum is only a 
restricted one. 

The question is: What function ¢(£) renders Cp a minimum 
for fixed C,? First, let ¢(1) be constant; a simple problem of the 
calculus of variation results which leads to the solution ¢” = 


constant, or a parabola, 
¢ = A TT Qé T Qe 


With this result in mind, put 


€ = awé — (Azo/c)4&(1 — E i 


where z)| is the maximum ordinate for ay = 0. Introduced 


in Eq. (5) and (6), this ieads to the result 


The effect of camber becomes evident in these formulas. A posi 
tive camber reduces the lift loss in the domain of the wing-tip 
Mach cones and thus increases Cy, but it also increases Cp on the 
whole wing. A gain in Cp at constant C, is possible. Evaluat- 
ing the optimal conditions by elementary means, it follows that 


Cpmin. Will be found for a camber 


Azo /C = ay/SA 10 


and that the gain over the flat plate (for the same lift) is 


Corp — Copmi 1 


Copp 1242? — 6A + 1 


rhe biggest gain is naturally expected at the lowest aspect ratios 
taking A = 1, which is the limit of the present theory, the gain is 
14.3 percent. For A = 2, the gain is only 2.7 per cent. 

The chordwise lift distribution for optimal conditions is sym 


metrical about the mid-chord, 


dC, 
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The results of this calculation are in excellent agreement 
with results of an approximate minimizing procedure given by 
Graham. ! 
standable general conclusion of Graham—namely, that camber 
yields only small reductions of the drag due to lift for higher as- 


The present example illustrates also the well-under- 


pect ratio supersonic wings 
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On Asymptotic Solutions for the Heat Transfer 
at Varying Wall Temperatures in a Laminar 
Boundary Layer with Hartree’s Velocity 
Profiles* 


Arthur N. Tiftord and Sheng To Chu 
The Ohio State University, Columbus Ohio 
May 22, 1953 


— IN A READERS’ FORUM NOTE of the same title! presents 
an independent derivation of the asymptotic solution de 
scribed in reference 2 and given in greater detail in reference 3. 
There are some ‘“‘fine’’ differences between the two enalyses, 
however, which lead to different evaluations of the ‘‘inherent”’ 
error. For that reason we reproduce some of our findings 
The differential equation under consideration is 


V5" + of Ya’ — of(2 — B)f/Va = 0 (1) 
with the boundary conditions, ¥7(0) = 1 and V%(-) = 0, 
where 
7’ = Hartree’s velocity function 
fii = exponent in the surface temperature representation 
Te — Ts = Shgx® 
7, = temperature acquired by the surface when insulated 
T,, = wall temperature 
¥; = thermal functions in the solution, 
T = T, + (:?/gce,J) ¥ + Dbz Vax” 
8 = pressure gradient parameter, 2m/(m + 1) 
i = independent variable, { {(m + 1)/2](u/xv)}" *y 
o = Prandtl Number 


We note that the second term, of l%’, is of the same order of 
magnitude as the portion of the third term, of’ }5. Therefore, 
when (2 — 38) >> 1—which is equivalent to 7 > 1 

term is negligible as compared with the third term. 
effect of the value of 7 on the thermal boundary-layer function, 
Vz, is exactly the same as that of the Prandtl Number, ¢. Just 
as for large Prandtl Numbers, the thermal boundary layer for 


* Research supported by the Aerodynamic Research Laboratory, USAF. 


the second 
Now the 
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large values of 7 is thin as compared with the velocity boundary 
layer. Consequently, the velocity profile, u/u:, may be repre- 
sented by the first few terms of the Taylor series, 


u/ = f’ = aon + (a;/2!)n? + 
where a, = f#"(0). (Note that the coefficients, do, a:, and ay, of 
the complete Taylor Series are zero for all laminar boundary 
layers on impervious surfaces. ) 

When 7% is extremely large, the first term, a.7, provides an 
adequate velocity profile representation for a nonseparating 


boundary layer. Under these circumstances, Eq. (1) reduces to 


Y®? —7yY =0 (2) 

where 7 = Ny, N = [o(2 — 8)asn|'/*, and the superscript (i) 
denotes the 7th order of differentiation with respect to 7. 

The solution satisfying the boundary conditions, }(0) = 1 and 


}( ©) = 0, is 
u\ /.. fs 4\ "4 ae ; 
Y =] 4r - y or (n/m) *{ Ki 3(2n /*/3) 
16// 3 3 


Y = 0.51769 '/?(K1/3(2n 2/3) | 


or 


where A;,/;3 is the modified Bessel function of the second kind of 
one-third order (see reference 4). 

At the time of our original analysis (1950) Dr. K. Pohlhausen, 
of the Aerodynamic Research Laboratory, USAF, suggested the 
investigation of a perturbative improvement upon the asymptotic 
solution. Such an improvement was developed. However, a 
more rigorous treatment was found to result from the use of an 
asymptotic series expansion for the function, }%. The latter 
treatment is given below. 

Let the asymptotic series expansion for }j% be 


l 
Fg f + — AS — BS | 


Substituting this series and the Taylor series expressions for f and 


7’ into Eq. (1), we obtain 


de a3 ‘ . I 
o (,*, ls 31Nt 7 t+ \Q » + no » 4 ) _ 
l a3 1 
ay («3 tow + Ne reg? ) a 


the independent variable having been changed throughout to 7. 
Collecting terms of the same power of V, the following system 
of equations is obtained: 
(s—1) (s—3) 
ye? —ay = ZT Cry -— = 
j=0 q=0 


d,n*y, t (4) 


where all the subscripts are positive integers, and 


TABLE | 
Values of (D> Y%)o for a Flat Plate 
o = 0.72 
Asymptotic | ” ” = “Exact’’ Values 
ad é (D35y3) Di ys) (Daye) (D5) (Daye): ( yo) 
Result | (Dz¥)o + 7Y3)0 (D5Y¥)o + (Anyao nyeo (D5¥)o + nya)o 4 = 6)o rt 4 t) poe ee 
DynY¥)o N3 N3 Ns N3 N6 N? 
10 —1.378 —1.404 —1.405 —1.405 —1.384 —1.423 
5 —1.094 —1.135 —1.137 —1.136 —1.106 -1.130 
4 —1.016 -1.063 —1.065 —1.065 —1.031 —1.053 
3 —0.923 —0.980 —0.984 —0.983 —0.942 —0.967 
2 —0.806 —0.881 —0.889 —0. 886 —0.832 —0.845 
1 —0.640 —0.759 0.785 —0.762 —0.682 —0.691 
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i=(s —j +1), k=(s —q-1) 
i+ 1 
C. a9 & d, ne ye = y 
a2(i!) (k!) 
The first equation in Eqs. (4)—i.e. s = O—is identically the 


asymptotic equation, Eq. (2), obtained previously. The remain- 


ing equations have the common formulation 


ye? — 1¥e = Fi[0; Yo-1, Vo-2 Yo; Ve-a, Ve-a™, yo? ] 


They can be solved analytically in succession. 

The solution for the case of the flat plate has been worked out 
in detail for o = 0.72. 
cally because in addition to ao, ai, and ay, the coefficients dz, de, 


This case is particularly simple numeri- 
47, Ag, Ayo, etc., are also zero. The resulting values of (Dy \7)o are 
given in Table 1. The corresponding ‘“‘exact’’ values obtained 
on the REAC analog computers at Wright Field and obtained by 
Chapman and Rubesin in reference 5 are listed for comparison. 
The original asymptotic analysis contains two different sources 
of error for applications at finite values of 7—-namely, (1) the 
assumption that the second term of the differential equation has 
a relatively small influence on the solution, and (2) the additional 
assumption that the velocity profile may be adequately repre- 
sented by the first term of a Taylor series. It is clear that the 
error due to (1) is independent of the Prandtl Number and de 
pends only on the value of (2 — 8) as compared with unity 
Errors due to (2), on the other hand, depend on both the value of 
7 and the Prandtl Number. Although the improvement ob 
tained by means of the asymptotic series approach tends to de 
crease both types of error, it inherently cannot fully eliminate 
the error arising from the first assumption. In other words, at 
large vaiues of n(2 — 8) the complete treatment is essentially 
exact; at smaller values of this parameter—e.g., values smaller 
than 10 
for the inherent error in the original asymptotic approach 


even the improved analysis cannot fully compensate 


ADDENDUM 


The May, 1953, issue of the JOURNAL OF THE AERONAUTICAL 
SCIENCES containing the criticism by Tribus® of Schuh’s Readers’ 
Forum note has just reached us. We offer the following com 
ment: 

It is true that the results of Lighthill’ are somewhat more 
general than those given by Schuh. Nevertheless, the implica- 
tion that the asymptotic analysis is therefore hardly worthy of 
mention is not justified. Despite appearances, the asymptotic 
analysis does not cover the same ground as did Lighthill. 

In addition to offering an independent evaluation of the surface 
rate of heat transfer, the asymptotic analysis contributes the 
following: 

(a) A derivation showing explicitly that the surface rate of 
(Although 


“* at large values of 7. 


heat transfer varies as 7 
Lighthill’s analysis may be developed further so as to affirm this, 
the published analysis does not contain such an analytic finding. ) 
(b) A novel application of ‘‘order-of-magnitude”’ analysis. 
(c) An analytic specification of the temperature profile within 
the boundary layer at large values of 7% (see references 2 and 3). 
Furthermore, the asymptotic analysis is more readily subject 
to analytic improvement through the use of a more correct specifi- 
cation of the boundary-layer velocity profile (see asymptotic 
series solution given above). Such an improvement should 
prove particularly significant when strong surface pressure gradi 
ents —especially adverse ones—are present. 
As a final remark, we bring to the reader’s attention the im- 


provement upon Lighthill’s analysis given in reference 8. 
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On the Behavior of a Sharp Leading Edge 


Y. C. Fung 
California Institute of Technology, Pasadena, California 
May 21, 1953 


B’ SHOWS THAT in a supersonic flow an airfoil having a 
sharp leading edge is unstable chordwise, in the sense that 
the leading edge will curl up. The supersonic part of the argu 
ment is unessential, the sharp leading edge will also curl up ina 
It appears that the idealization of a sharp leading 
As the 
concept of sharp leading edge is prevailing in aerodynamics, it 


subsonic flow. 
edge is unsatisfactory as a model for aeroelastic analyses. 


may be worth while to point out the elastic behavior of such an 
edge; it is important in comparing aerodynainic theory and ex 
periments that require an aeroelastic model. 

Consider a two-dimensional wedge of perfectly elastic material, 
as shown in Fig. 1. If one side of the wedge (6 = 0) is loaded bya 


normal pressure distribution 
p(r) = Pym ] 


while the other side (@ = a@) is stress-free, the Airy stress function 
o(r, 8) can be expressed in the following form: 


(1) Him = 0, 
@ = r™*? la cos (m + 2)0 + b sin (m + 2)0 + 
ccos m6 + d sin m6} (2 
where 
a = —Km [m+ 2 — cos2(m+ l)a - 
(m + 1) cos 2a] /2A 
b = Km {sin 2(m + 1l)a + (m + 1) sin 2a] /2A 


c = K —a, d = —(m + 2)c/m 
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A = —m(m + 2) — cos 2(m + 1l)a + (m + 1)? cos 2a 
K = P/(m + 1)(m + 2) 
(2) If m = 0, 
¢ = Pr? [—tan @ cos? 8 + (1/2) sin 20 + a — 6]/2(a — tan a) 
(3) 
where r, 9 are polar coordinates and a is the wedge angle. m = 0 
corresponds to a uniform pressure on the wedge; m = l, a 


linear pressure distribution that decreases to zero at the leading 
edge. Eq. (2) is a simple generalization of a result given by 
Levy;? it is valid for real or complex values of m for which A # 
0 

The elastic behavior of the wedge at the leading edge, given 
by the limiting case r — 0, is as follows, where a6, oy, tye denote 
the normal and shearing stresses and u, v the components of dis 
placement in the direction of increasing r and 6, respectively: 

(a) oy, o@, Tre tend to zero if m > 0; tend toa finite value if 
m= (0; and tend to infinity if m < 0. 
(b) The slope ov/Or is finite if m 

m <0. 
(c) The curvature 0°7/0r? is finite if mm > 1, but tends to in 


> 0; it tends to infinity if 


finity if m <1. 

Note that the conclusions (b) and (c) agree with those of 
simple beam formulas, by considering the wedge as a beam of 
variable cross section. Hence, the above more elaborate analysis 
shows that, for small a, the beam formulas are justified either 
when the tip of the wedge is sharp or is cut off. 

For example, in a subsonic flow the leading edge is a stagnation 
point, but small changes in local slope do not strongly change the 
< 0, and it is 


pressure distribution. Hence, one may take m 


seen that an accidental disturbance will cause a sharp leading 


edge to curl up. 
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The Effect of Dissociation in the Stagnation 
Region of a Blunt-Nosed Body 


Ivan E. Beckwith 
Langley Aeronautical Laboratory, NACA, Langley Field, Va 


May 21, 1953 


I* A RECENT PAPER BY Moore,! the laminar boundary-layer 
characteristics with equilibrium dissociation were computed 
for a flat plate. An interesting conclusion obtained was that the 
heat-transfer rate per unit area is essentially unaffected by dis- 
sociation when the plate temperature is below 2,500°R. The 
question naturally arises as to whether the same conclusion can 
be made concerning the heat transfer in the stagnation region of 
a blunt-nosed body. 

The purpose of this note is to discuss an approximate method 
for the calculation of the heat transfer on an arbitrary two- 
dimensional body in the presence of equilibrium dissociation. As 
a numerical example, the heat transfer is computed in the stag 
nation region of a blunt-nosed body. 

The method used is similar to that of Kalikhman,? except for 
some details that have been modified so as to include dissociation 
The boundary-layer equations are written as follows: 


Momentum, 


FORUM 645 
Ou Ou Op ra) Ou 
pu + pv =-— i Me 
Ox oy Ox oy oy (1) 
O = Op/oy 
Continuity, 
(0/Ox)(pu) + (0/dy)(pv) = 0 (2) 


Energy, 
oH 


_ °( oH 4 ) 6 
; oy oy M oy oy , Oy\a 


OH 
“ 
Ox 


pu + p 


- i) | (3) 


where H = h + (u?/2) is the stagnation enthalpy, A is the static 
enthalpy, and a = (u/k)(dh/dT) is the Prandtl Number. 

Integration of Eqs. (1) and (3) across the boundary layer and 
the use of Eq. (2) result in 


d am i u? pe tes u? 
f, ( - Jan + 2 | - ) an + 
dxJ° ly uu" 0 uy u,? 
a 1 1 di UwPw| Of1 
f (2 - “an ¢ 1 a (" (4 
0 p My u, dx yp, **|_ On \u, " 
d e a 1 du, 1 dh, Pm og 
f, 1 — dyn + _ J 
dxJ0 it) ag u,dx t,* dx 0 My 
Sea wiu oO /t* 
(: - Jan = =~ | ( ) | (5) 
t,* Opi 2, | On \t,* e 


where 7 is the Dorodnitsyn variable 


y p ' 
dy (6 
0 pi* 


Hl — h, is the difference be 


p* is the stagnation density, t* = 
tween the local stagnation enthalpy and wall enthalpy, and sub 
scripts 1 and w refer to local free-stream values (behind the bow 
shock) and wall values, respectively. 

The usual fourth-degree polynomial is used for the velocity 


u 12 +A " rr 3A — 12 6—A F. 2 
- E yy -—-9™ + : 9" “aah, ie (7) 
My 6 = 6 6 


n/A (Ais the nominal velocity boundary-layer thick- 


Pi Pi . 2 du; - 
= A? 
HMw \Pu dx 


Eq. (8) satisfies Eq. (1) evaluated at the wall with the restriction 
2,500°R., allowing the additional assumptions that, 


where n* = 


ness) and 


(8) 


that 7, < 
in the immediate vicinity of the wall, 


wp = (pe/Te)T 


Gul A1/T) 


(9) 


p= (10) 


i.e., the wall temperature is always less than the dissociation 
temperature so that the conventional perfect gas assumptions 


can be used here. 
The stagnation enthalpy profile is also taken as a fourth-degree 


polynomial in 7* 


a 
- «Be? + Bae F(t — 3B, — 2Bz)n** 4+ 
ty 
(2B, + Bz, — 3)n** (11) 
where 
(12) 


uy? (= + *) 

B. = (1 — ow) - 
2t,* 6 

which satisfies Eq. (3) evaluated at the wall 
used in Eqs. (11) and (12) is that the thermal or enthalpy bound- 
ary-layer thickness is equal to the velocity boundary layer thick 
ness A. Calculations made for a perfect gas indicate that this 
assumption gives results in close agreement with the exact com 
putations of Brown and Donoughe,’ as long as the ratio of the 


The assumption 
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surface temperature to stream temperature is less than about 0.7. 
Thus, instead of using the ratio of thermal to velocity boundary- 
layer thickness and \ as the unknowns in the solution of Eqs. 
(4) and (5) as in reference 2, B; and ) are used as the unknowns. 
Before the solution can be carried out, some approximation 
must be made for the density profile p;/p appearing in the third 
integral of Eq. (4). The dissociation data given by Moore can 
can be closely approximated by a fifth-degree polynomial in 4 


ee eee b( =) 4 1 a) 4 » (2 )'+6(7 )’ 
p he he he hep hp, 

(13) 
where hg = 96.8 B.t.u. per lb. and, for hy = Cp,Tw, 


h Cool a u,? fu \? 4 (7: he + u,? \t* 
hg hg Qhe\m he hg 2he/t,* 
in terms of the velocity and stagnation enthalpy profiles. 


Application of the above results to Eqs. (4) and (5) for the 
limiting case of u, — O (stagnation flow) results in the following 


22(2 é r n 
pi \10 120 
ppl12+2X pp 2 


_ rd 5a? 
-(37 --- (14) 
315 3 144] °* 


two simultaneous equations for \ and B;: 


k + kB, + k,B,? + k;B,8 + k,B,* + kB = 


Pe Gr Pl 
owrdla + (A/6)as3] : 
B, = : (15) 
(pi/Pw) — owdlar + (A/6)as) 
k, ki, ke, , are constants depending on the values of },, be, : 
and the wall temperature and stream enthalpy; a, a, , are 


also constants. 
The convective heat g transferred from a unit area of surface 


per unit time is 
dw = —k,(O7T/0y) 


or using the definitions of ¢*, ¢,*, A, Bi, and Eq. (6) 


Re p, L2du, | 
——§,° 5, 
Cs. ey h OS 


Defining a wall. Nusselt 


qwulb = — (16) 


where L is a characteristic dimension. 
Number as 


Nuw = quwlhe/RwlTw — (hi/Cp,)] 


I 

there results 
Nix 
3, Pi/ Pu 


Via 


This is plotted in Fig. 1 versus h;/hg for 
Ty» = 2,000°R. with and without dissociation. 


; (17) 

V (pwl.?/tw)( du; /dx) 

for stagnation flow. 

The correspond- 
ing density variation is also shown. 

The that 

equilibrium dissociation has only a small effect on the actual 


conclusion obtained here for stagnation flow is 
heat-transfer rate per unit area, provided the wall temperature is 
less than the dissociation temperature (3,000—4,000°R. ). 

The validity of this conclusion can be indicated by using the 
above methods to compute the heat transfer on a flat plate with 
and without dissociation and comparing the results with those of 
Moore.!' The expression obtained by the present method for 


the heat transfer on a flat plate is 


2 
a (i. + 
kv 2 


Gu = 2” Pu PTS DES 1/ 
“te 37 315 pop x \ '/2 
eee 2 — 
630¢,, 37 i) 


Thus, since the only values appearing in Eq. (18) depend on 
the wall temperature, stream velocity, and stream enthalpy, 
there can be no effect of dissociation so long as both the stream 


- Cs) + 2a(l — oy)? 
(18) 
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at the forward stagnation point of a blunt-nosed body with and 
without dissociation, The lower curves show the corresponding 
density changes. 

and wall temperatures are below dissociation values. This re. 
sult is not surprising, since by the present method the velocity 
and enthalpy profiles depend only on wall and stream conditions, 
It is to be noted, however, that the temperature and density pro- 
files would be different for dissociation even with the same wall 
and stream conditions, since the enthalpy profile is invariant 
Apparently, the temperature and density profiles either have little 
effect or compensating effects on the heat transfer at the wall of 
both a flat plate and blunt body. 
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Heat Transfer in a Transitional and Turbulent 
Boundary Layer at Supersonic Speeds 


William B. Fallis* 
Institute of Aerophysics, University of Toronto, Toronto, Canada 
May 11, 1953 


7 VARIATION IN TEMPERATURE RECOVERY and heat trans- 
mission with Reynolds Number in the laminar, transitional, 
and turbulent boundary layers on cone and cone-cvlinder models 
has been recently reported by Eber.' His temperature measure- 
ments were made under transient conditions in the NOL 40- by 
40-cm. Aeroballistics Tunnel. des Clers and Sternberg? have de- 
termined local recovery factors in the transitional and turbulent 
boundary-layer regions under temperature equilibrium condi- 
tions in the BRL tunnel also using cone and cone-cylinder 
The Addi- 


tional recovery factor and heat-transfer data® are available from 


models results do not support Eber’s findings. 


transient temperature measurements made at a Mach Number 
of 2.5 in the transitional and turbulent boundary layers of a flat- 
plate model mounted in the UTIA 16- by 16-in. tunnel. Good 
agreement with Eber’s results has been obtained 

Two geometrically similar models were used in the experi- 


ments: one a thermally insulated flat plate and the other a con- 


Aircraft 


* Now, Vultee 


Corporation, Fort Worth, Tex 


Senior Aerodynamics Engineer, Consolidated 





ductin 
10 in 
along 
yals b 
compl 
The 
made 
coppe! 
surfac 
tially 
were | 
tempe 
speed 
then f 
statio 
were ¢ 
Loc 
ature 
as a | 
consic 
curve 
in rec 
recov 
lamin 
transi 
repre: 
f= 
the fi 
ture 
phote 
layer 
creas 
the v 
out t! 
Th 
cond 
to its 
The | 
surro 
coup! 
surfa 
heate 
temp 
flow 
tunn 
moco 
it pre 
Th 
of th 
abati 
Num 
respc 
selt 
for a 
turbt 
prese 
flat-y 
appr 
in th 
is in 
valu: 


E 
Missi 
Wind 
Janu: 


08 





06 


PS Pw 


04 


02 





Y ratio 
ith and 
vonding 


‘his re 
elocity 
litions, 
ty pro- 
1€ wall 
ariant 
€ little 
vall of 


ns fora 


Journal 
NACA 


minar- 
ies Are 


ent 











Ti- 
on- 


aft 











READERS 


The model span was 8'/2 in. and the chord length 


ducting plate 
10 in. Eight plugs were inserted, flush with the surface, 


along the centerline of each plate and were spaced at 1-in. inter- 
yals beginning 1'/> in. aft of the leading edge. The plugs were 
completely interchangeable at all stations. 

The thermocouple plug used with the nonconducting plate was 
made of bakelite with a thin copper disc set into its surface. A 
copper-eureka thermocouple junction was soldered to the under 
surface of the copper insert. This thermocouple plug was ini 
tially installed in the farthest aft station, and nonconducting plugs 
were inserted in the remaining stations. The insulated surface 
temperature was recorded for a series of tunnel runs by a high 
speed photoelectric potentiometer. The thermocouple plug was 
then removed and interchanged with the plug in the next forward 
station. The process was repeated until temperature-time traces 
were obtained at each station. 

Local recovery factors determined from adiabatic wall temper 
ature measurements on the insulated plate are plotted in Fig. | 
as a function of Reynolds Number. The recovery factors are 
considered accurate to within +1 per cent. The experimental 
curve may be divided into two regions: one where the variaiion 
inrecovery factor is large and the other where it isconstant. The 
recovery factors in the first region are higher than the theoretical 
laminar recovery curve and indicate that the boundary layer is 
transitional. In the second region the experimental data are 
represented to within 1.4 per cent by the theoretical relation 
r = (Pr). * derived for turbulent boundary-layer flow, where 
the fluid properties are evaluated at the adiabatic wall tempera- 
ture. A liquid film method of flow visualization and schlieren 
photographs were also used to identify the turbulent boundary- 
layer region. The recovery factor in the transition region in- 
creased to a peak value approximately 4.5 per cent greater than 
the value in the turbulent boundary layer. These results bear 
out the work of Eber! and Stalder. 

The heat-transfer flat-plate model differed from the non 
conducting model in that a sheet of '/s-in. thick copper was fitted 
to its surface. The centerline plugs were also made of copper. 
The thermocouple plug consisted of a '/s-in. thick copper core 
surrounded by an ebonite insulating jacket, with the thermo 
couple junction located midway between the upper and lower 
surfaces of the core. The conducting surface was uniformly 
heated to a temperature of about 40°C. above the adiabatic wall 
temperature. Heat-transfer measurements were made with heat 
flow occurring over the entire length of the plate. A series of 
tunnel runs and temperature records were made with the ther 
mocouple plug installed in the trailing-edge station and then with 
it progressively moved forward toward the leading edge. 

The local heat-transfer coefficient, calculated from a knowledge 
of the rate of cooling of the conducting model and the local adi 
abatic wall temperature, is plotted in Fig. 2 as the local Nusselt 
Number versus Reynolds Number, using the fluid properties cor 
responding to the adiabatic wall temperature. The local Nus 
selt Number predicted by the methods of Eckert and Drewitz 
for a laminar boundary layer and von Karman and Colburn for a 
turbulent boundary layer are shown for comparison. The 
present experimental data indicate that the boundary layer on the 
flat-plate model was transitional for Reynolds Numbers less than 
approximately 4 & 10°. The Nusselt Number for heat transfer 
in the fully developed turbulent boundary layer on the flat plate 
is in accord with Eber’s cone-cylinder turbulent boundary-layer 
values and is well represented by the relation due to Colburn. 
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On the Operational Form of the Linearized 
Equation of Supersonic Flow 


L. E. Fraenkel 

Department of Aeronautics, the University, Glasgow, 
Scotland 

May 25, 1953 


AS IS WELL KNOWN, the linearized equation of supersonic flow, 


, Oo 0° 0*o (1 
: _ 7 = ) 


"Ox? oy" Os? 
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may be written in operational form 
B*p* — (0*%6/dy?) — (0°/0s?) = 0 (2) 


where $(p, y, z) is the Laplace transform of the perturbation po- 
tential ¢(x, y, z). However, the rigorous deduction of Eq. (2) 
from Eq. (1) is not so simple as it might at first appear, and the 
only justification given by many writers (see, for example, refer- 
ences i, 2, and 3)—namely that ‘‘¢ and 0¢@/0Ox vanish at x = 0” 
is both insufficient and incorrect in the present case and may well 
mislead the uninitiated (as it did the present writer). 

The first point to be noted in considering the transformation of 
Eq. (1) to Eq. (2) is that the relevant boundary values are those 
on the downstream side of the foremost Mach surface, which 
bounds the field of disturbed flow; and on this side @ vanishes, 
but its derivatives in general do not. 
fields there are a number of surfaces across which ¢ is continuous 
while its derivatives are discontinuous, and all these discontinui- 
ties must be taken into account when the operational form is 


Furthermore, in most flow 


derived. 
We take the free stream to be in the direction of increasing x. 


The surfaces of discontinuity are denoted by x = fi(y, 2), (7 = 
0, 1 The Heaviside operator is defined by 


p'= f ; dx’ 
foly, 2) 
z) 


and the corresponding Laplace transform of any function F(x, y, 2 


«#8. 


1S 


F(p, y,s) = p {, e-™ F(x, 9,3) dx 
of Joly, 2) 
We also adopt the notation 


F y, zs] = AF(f;, y, 2) 


fly, 2) +, 9,2] — Fifty, 2) —. 3 





We now multiply Eq. (1) by pe?” and integrate from x = 

foly, s) to ©, assuming that the resulting integrals exist for all 

the lines y = constant, : = constant on which Eq. (1) is satisfied 

forallx. [Wecannot expect to justify Eq. (2) for values y = y,, 

z = 2, such that there are points x = x1, y = Wi, 5 = 21, where the 

solution will prove singular and Eq. (1) itself is not satisfied. | 
Integrating twice by parts, we obtain 


° —1 "fii 9 
Oh ° O° 
p e pr dx = > p 4 ” dx + 
JS fo Ox? i=0 fit x" 
F 7 O° “ ? Op 
p eo? — dy = = pe Pl. Fis. He 2) + 
fick Ox? 2d, Ox 


p*o(p, y, 5) (3) 


Differentiating ¢ twice under the integral sign, we also find 


7 


‘ 1 f , 
ot. TS ci 
2. p e "pdx + p eo dx 
oy? Oy? | ;=0 fi aes 
oo 0% n ; do of, 
=p ” eal - dx — 2 pe Pri.a (Sis ¥, 3)° (4) 
s oy? i=0 oy oy 


the expression for 0°¢/0=s* being similar. Hence Eq. (1) becomes, 


by Eqs. (3) and (4), 


n 


O*p " Om . 
> pe] Bta— (fi, ¥, 5) + 


Ov? 


B*p*s — _— a at 
Oy" 0=* i=0 Ox 
Op Of; Od Of; 
A (Jie ¥ @)° + A (f;, y, 5)° = (5 
oy” oy Os j oz | 


We now complete the justification of Eq. (2) by proving that, 


for each /, 


,, Ob Op Of; OP. Of; 
B7A G39 2) + A (fi, ¥ Z)° + A (fis ¥, 2) = () 
Ox oy oy Os Oz 
(6) 
We introduce variables (£, 7, ¢) such that — = O on the Mach 
surface x = fi(y, 2), 
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té=x-—fily,2), n=y FF =32 
and we also write $(x, y, z) = ¢(£,7, ¢). Then, 
dp/Ox = Iy/dE 
Oo O¢g Of; o¢ 
= — oe ery 
oy Of Oy On 7 
Oo » _ Of; m oy 
Os OF OS oe 


and, since the equations of continuity and momentum ensure 
that @ and the tangential velocities are continuous across the 
surface x = f;(y, s), we have 


o¢ d¢ 
A (O,”,¢) = A 
On of 


(0,7, ¢) = 0 g 


Hence, the left side of Eq. (6) becomes, by Eqs. (7) and (8), 


o¢g of;\? of;* 
4 (0, 7, ¢) | B? — - - 9 
Of oy Oz 


Now the equation of any Mach cone may be written 


2@antey (y — 91)? + (s — 2)? 
so that 
Ox ste Bly — yi) = Ox = 
vy VWiy—n)?+(s—a)? % t 
B(2 — 3) 
Vivg—-AP Tes 

whence 

(Ox/Oy)? + (Ox/Os)? = B? 10 
But the surface x = fi(y, s) is the envelope of a family of such 


Mach cones, so that Eq. (10), which is a relation between local 
surface slopes valid at every point of these constituent Mach 
cones, must also be valid on the Mach surface. Hence, the ex- 
pression (9) is zero, Eq. (6) is proved, and the operational form, 
Eq. (2), is justified. 

Finally, it may be worth mentioning that, when the generalized 
Fourier transform is applied to the differential equation (1), the 
above discussion still applies with only very slight changes. 
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On the Velocity Distribution in the Laminar 
Sublayer and Transition Region in Turbulent 
Shear Flows 


Francis R. Hama 

lowa Institute of Hydraulic Research, State University of lowa, lowa 
City, lowa 

May 4, 1953 


I HE WELL-KNOWN RELATION on the mixing length 
l = ky i 
is chiefly based upon Nikuradse’s experiments with flow in pipe 
His result, however, is somewhat ambiguous in the region close 


to the wall. First, the measurement is done outside the pipe’s 
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exit but not in the pipe. In fact, his result shows a finite velocity 
Second, the viscous term is neglected in the calcula- 


Furthermore, the logarithmic formula tends to minus 


aty = 0 
tion of / 
infinity near the wall 

The writer! has calculated the mixing length from the velocity 
distribution in the turbulent boundary layer along a flat plate by 


using 
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l 
u re) ‘ ul 
( l : / ( 1——)dy (2) 
UJ ox Jo l s 


r = pl*(du/dy)? + p(du/dy) (3) 


and obtained a universal relation 


i* = (ay); 7” ie, ¥° yu*/v, v*? to/p (4) 


in the viscous region (laminar sublayer and transition region in 
which the viscous term is not negligible), a being a numerical 


constant, as shown in Fig. 1. 
Once the universal relation, Eq. (4), is established, the velocity 


distribution in the viscous region can easily be obtained. Putting 
Eq. (4) into Eq. (3) and assuming +r = to = constant, we have 
the differential equation 

s'(dg/dz)? + 4(d¢/ds) {=0 (5) 


for which the solution is 


21i1-V 1 +325 1 
c= + 7 l (6) 
: | 33 : (. 2’ ) | 


in which ¢ V/ 2au*, s V 2ay*, and F is the elliptic integral 


of the first kind defined by 


~ dy f dz (7 
=2 
o0 V1 — (1/2) sin? ¥ o Vi¢+s 


cos ¥ (1 — s?)/(1 + 2? 


The numerical value of ¢(z) is tabulated in Table 1 and plotted 
in Fig. 2. This solution includes the laminar sublayer profile 
” a 7s u* u/v* 8 


i“ 


for very small y* and can be connected smoothly to the log 
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TABLE 1 

4 ¢(z) 
0 0 

0.2 0.200 
0.4 0.397 
0.6 0.595 
0.8 0.785 
1.0 0.960 
1.4 1.250 
ee 1.464 
2.2 1.621 
2.6 1.738 
3.0 1.829 
$0) 1.981 
5.0 2 O77 
6.0 2.141 
7.0 2.187 
8.0 2.223 
9.0 2.250 
10.0 2.212 
20 2.472 


arithmic distribution. The experimental results for the flat 
plate boundary layer coincide well with Eq. (6) if the value a = 
0.102 is chosen. These facts and calculations have already been 


given by the writer.! 


Now, if the relation (4) is also valid for channel or pipe flow, 
the velocity profile, Eq. (6), must be a good representation for 
that case, since the other assumption rt = 7) = constant is allow- 
able in the viscous region for channel or pipe flow. Comparison 
with some of the experimental results? is made in Figs. 3 and 4, 
in which @ is chosen so that Eq. (6) connects smoothly to the 
logarithmic distribution for the channel flow or given arbitrarily 
for the pipe flow, since no logarithmic formula is given in the lat- 


ter case. 


In conclusion, the universal relation (4) is now well confirmed, 
and Eq. (6) gives a good representation for the velocity distribu- 
tion in the viscous region either for boundary-layer or for channel 
or pipe flow. By use of the velocity distribution given here, a 
refinement in the calculation of friction or heat transfer can be 
expected. The numerical constant a@ cannot be determined 
definitely at the present stage, because of the lack of a consistent 


logarithmic law. 
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The Medium in Fluid Mechanics 


L. S. Dzung 
Brown, Boveri & Cie., Baden, Switzerland 
May 29, 1953 


I IS USUAL TO ASSUME certain idealized properties of the me- 
dium to simplify the mathematical theory of fluid mechanics. 
The most common assumption, in case the fluid is gaseous, is 
that of ideal gas with constant specific heat—i.e., polytropic 
gas.! 
Traupel? shows that the equations of a polytropic gas 
pu = RT (1) 


K Cp/Cy = const (2) 
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can in most cases be relaxed and substituted by the less stringent 
ones, 
Z pv/RT = xs) (3 
k —(0 log p/O log v), = const (4 
While Eq. (3) requires only that the ‘‘compressibility factor” zis 
function of the entropy s alone, Eq. (1) requires that ¢ is a con. 
stant. The isentropic exponent k becomes the ratio of specific 
heats « for an ideal gas. 
It can be shown, however, that Eq. (3) alone is sufficient to 
define a gaseous medium similar to the ideal gas defined by Eq 
(1). The additional condition, Eq. (4), corresponds to the addi- 


tional Eq. (2) 


Define 
a (O log v/O log 7), 
8 (O log p/O log 7) 
¥ —(0O log p/O log v)7 


From thermodynamic identities, it is known that 


z (28) = kK ; . 
olan, «-1' ” 


(**) , (°') , ; 
= ( — wy, = (p — ) (6 
Op) 7 ov / 7 P 


where /# and u denote enthalpy and energy, respectively. Now 
differentiate pu = sRT, 
dp dv dz dT 
a aa 
p v : 7 


1\d 1: iT 
(: - )? -' = (s = const 


k-1 (2) ee (=) 
k par). | s\ar/, 


Substituting Eqs. (5) into the left-hand side of the last equation 


ay=f8, k= KY, 


or 


and denoting the right-hand side by Z results in 


k-1 K k—1 : - 
° B= ‘a = Z (7) 
k k— | k— 1 


For an isothermal, Eqs. (6) give 
vdp = [1/(1 — a) dh 
pdv = [1/(8B — 1)]d(h — pr) 


or, after adding, 


l ] I 
i(pv) = dh — 1( pv) 
ap) (, — a ' 8 - ) B— !1 - 


and by substituting a and 8 from Eq. (7), 


cm ] 
fi - 254] 
k—1 


-d(pv) (T = const.) (8) 
b-« aed ; 


dh = 
Now for an isentropic, dh = vdp. This can be written, to- 
gether with the definition of k, as 


dh = -d(pv) (s = const.) (9) 
k l 
Eqs. (8) and (9) become identical if Z = 1-—i.e., if (07/07), = 0 


This is Eq. (3). In this case 


dh = [k/(k — 1)]-d(pv) (10) 
Without restriction as to the constancy of a parameter 
Define 
T pv/R et, AX k/(k — 1) (11) 
and consider 7 and s as the independent variables of state. Eq. 
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(10) can be written more precisely as 


(Oh/Or), = AR, (Oh/ds), = 0 


Differentiating the last equation with respect to 7 and using the 
first one results in 
07h/O07rds = (OA/Os),R = O 


The last equations show that / and X (or k) do not depend on s 
and are, therefore, functions of 7 alone, 


h = h(r), k = k(r) = R(h) (12) 


These properties are exactly those of an ideal gas. 
Define further a quantity o corresponding to the entropy by 


do = ds/2(s). If now Eqs. (10) and (11) are substituted into 


Tds = dh — vdp, one obtains 
do = \R(dr/r) — R(dp/p) 


Because of Eq. (12), o is ‘‘separable’’ as in case of the entropy 


of an ideal gas. 

In fact, a medium defined by Eq. (3) has all the mathematical 
properties of an ideal gas as discussed in reference 1, Section 4, 
if T is substituted by 7, s by o, c,p/R by A, c./R by X — 1, and x 
by k&. It seems, therefore, appropriate to call it an ideal vapor 
or a perfect vapor, following Leib.‘ If, in addition, Eq. (4) is true, 
then the medium can in most cases be treated as a polytropic 
gas and may be called a polytropic vapor. In this case, Eq. (10) 
can be integrated explicitly to give h = Apv, corresponding toh = 
cpl for a polytropic gas. The integration constant has been set 
equal to zero. 

Traupel? mentions that the equivalence breaks down if heat 
conduction is present. This situation can be saved sometimes, 
at least formally, if the gradients of JT and of h can be considered 
as parallel. It is then possible to define an ‘‘enthalpy diffusivity” 
which relates the heat flux to grad. h as the thermal conductivity 
is related to grad. J. This diffusivity, unlike the conductivity, 
is, however, not a unique function of state unless the state path 


is prescribed 
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Beitradge zur Thermodynamik der realen Gase, Schweize, 
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Comments on Location of Detached Shock 
Waves Ahead of Plane Bodies 


Robert Wesley Truitt 

Professor of Aeronautical Engineering, Virginia Polytechnic 
Institute, Blacksburg, Va. 

May 25, 1953 


I’ HIS COMMENT on my Readers’ Forum paper,' Mr. Moeckel 
(in reference 2) correctly points out that both papers of refer- 
ences | and 3 are based on the same assumption concerning the 
body sonic point and, further, that experiment and theory shown 
in reference 1 confirms his geometrical method. My conclusions 
made in reference 1 were not intended to disprove Mr. Moeckel’s 
Original theory but rather to indicate the direct effect of nose 
angle on shock detachment distance. As stated in reference 1, 
the purpose of that paper was ‘‘to derive an expression for the 
shock detachment distance which is an explicit function of the 


opening angle of the wedge,’ which was done. 
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In reference 4, Laitone erroneously states that the results of 
references | and 3 are only valid for relatively blunt bodies at the 
higher transonic Mach Numbers. Indeed, it is quite well known 
that for bodies with well-defined shoulders, such as slender simple 
wedges, the body sonic point is located at the shoulder, and, con- 
sequently, the characteristic distance is obviously from the fore- 
most point of the shock to the shoulder of the wedge. This fact 
was the basis of the methods of references 1 and 3 and was shown 
in Fig. 2 of reference 1 to be in excellent agreement with experi 
ment for very slender bodies. Laitone further claims that the 
only results valid for very slender wedges (@ < 5°) at Mach 
Numbers near unity are his, given in Eq. (3) of reference 4 
Laitone’s results are plotted in Fig. 1 along with theoretical and 
experimental data of reference | for a 4'/2° wedge. It is clear 
that Laitone’s results are in very poor agreement with experi 


ment. On the other hand, theories of references 1 and 3 show 


excellent agreement. 
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On the Transonic Drag of a Cone Cylinder 


John W. Miles 
Department of Engineering, University of California, Los Angeles 
Consultant, Douglas Aircraft Company, Inc., Santa Monica 


June 3, 1953 


_ DIFFERENTIAL EQUATIONS governing the axisymmetric 
flow past a body of revolution in an otherwise undisturbed 


stream of near sonic velocity l’ are 
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upy—v, = 0 (1) r 
r—-\rv), = [((M? — 1) + (7 + I)uju, (2) \ 
where the total axial (x) and radial (7) velocity components are \rzf ( x) 


U (1 + #) and Uv and where y is the usual specific heat ratio. 
We consider an approximate solution of these equations for a 
slender body (6 << 1) consisting of a conical head of semicone 
angle 6 and unit altitude and a cylindrical afterbody of constant 
radius 6, as shown in Fig. 1. 


If it is assumed that the free-stream flow is only slightly super- 
sonic or, more precisely, that the parameter 


k = (M2 — 1)/(y + 1)8 (3) 


is a small positive number, there will be a detached shock far 
upstream of the body, and the flow will be subsonic between this 
shock and a sonic line, r = f(x), which terminates at the shoulder 
(x = 1,7 = 5) of the cone-cylinder. The boundary condition at 
this sonic line, obtained by requiring the /ocal Mach Number to 
be unity there, may be written 


1 
u+-—v2 = —ké?, r = f(x) (4) 


Moreover, it is evident that the flow at the shoulder is locally 
Prandtl-Meyer, whence the sonic line must be perpendicular to 
the conical surface at the shoulder, viz. 


f(1) = 6, fl) = -é"! (5) 

In the neighborhood of the body the right side of Eq. (2) may 
be neglected as 0(6'); moreover, it is never large compared with 
the left in the region 6 << x < 1,0 <r <f (x), except possibly: 
(a) near x = O, where there is a stagnation point at the nose of 
the cone; and (b) in the neighborhood of the shock. 
ingly, we pose the slender body approximation 


aes ee r | _ ax*l 
“u= és ) k + In Li | - oy? ( (6a ) 


v = 6*x/r (6b) 


Accord- 


which satisfies Eq. (1), Eq. (2) with right side neglected, Eq. (4), 
and the boundary condition v = 6 on the cone (r = 6x). 
The pressure on the cone is given by 


pb = — pU2[u + (1/2)0?}, <4, (7a) 


a | l (= 
= pl%?<k +1 = 9a> =b) 
. ) . bx + 2 f2 ; { (7b 


The corresponding drag coefficient, based on the free-stream 
dynamic pressure (1/2)pl and the base area 76%, is (after inte 
gration by parts) 


* 
. f"(x) 
Cp = 26 [ ~~ alt dx (8) 
0 f(x) 


The term 6*x?//? in Eq. (7) makes only a higher order (in 5) con 
tribution to Eq. (8). 


The integral in Eq. (8) can be evaluated approximately by 
noting that f(x) is 0(6~') except near x = 1, where it is 0(6), so 
that the contribution of the integrand is concentrated in the 
neighborhood of x = 1. It follows that only a small error will 


be committed by replacing the slowly varying numerator x?f'(x) 


by its value at x = 1, viz. —6~!, and replacing f(x) by the first 
two terms in its Taylor expansion about x = 1, whence 
1 , 
CS) , dx 
x? dx =~ —§™! ~In 6? (9) 


0 f(x) 9 €— 6-Kx — 1) 


\ 
~~ * 
a, 


U Utieu) \ ' 
or am eee, 
i 














(SHOCK) /” 
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Fic. 1. Cone-cylinder in slightly supersonic flow 


|We emphasize that the approximation (9) does not assume f(x 
to be everywhere a straight line.] Substituting Eq. (9) in Eq 


(8) yields 
Cp = 6°(2k — 4 1n 5) (10 


The result, Eq. (10), may be subjected to two general checks 
Firstly, provided that y is held constant,* it satisfies the similar 
ity law proposed by Oswatitsch and Berndt,!' viz., 6~? Cp + 41né 
depends only on the similarity parameter k. Secondly, the slope 


of the Cp versus 1/ curve at J = 1 is given by 
OCp ou we = 4 ly + 1) (11 


as should be expected from the general considerations advanced 


by Liepmann and Bryson.2. More directly, 


0u/OM| wer = —([2/(y + 1)] 

The present problem has been solved numerically by Yoshi- 
hara’ for AJ = 1. 
rule to the integration of the pressure distribution over a 5° cone 
(Fig. 7 of reference 3) is Cp = 0.0758, while Eq. (10) yields 0.0747 
Since the difference (1'/. per cent) between these two numbers 
is of the same order as 6‘, the agreement is within the expected 


The result obtained by applying Simpson's 


accuracy of the small disturbance theory, and it appears that 
Eq. (10) is a sufficiently 
small values of & and 6. 
that this is in 


satisfactory approximation for 
Nevertheless, it must be emphasized 
fortuitous, and Eq. (10 
0( 6"), though 


moreover, it is not 


some measure 


probably still neglects some term of even 
the term of 0(6? In 6) is correctly predicted; 
reliable for values of & sufficiently large to bring the intersection 
of the sonic line and bow shock close to the body. Since shock 


attachment occurs for values of & near unity (1.4 for a 5° cone), 


this last restriction implies k< 1. 
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WADC TR 52-295, Wright Air Development Center, Dayton, Ohio, No 
vember, 1952 

* It is not clear whether the similarity law of reference 1 is valid with 
respect to changes in y in the present problem due to the local Prandtl 
Meyer expansion at the shoulder. On the other hand, the requirement 
that the sonic line be normal to the cone surface is not entirely consistent 
with the slender body approximation 
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Similar Solutions of Compressible Boundary- 
Layer Equations 


Ting-Yi Li and Henry T. Nagamatsu 
California Institute of Technology, Pasadena, Calif. 
lune 1, 1953 


rather 


They are, in terms of the conventional sym- 


HE COMPRESSIBLE BOUNDARY-LAYER EQUATIONS are 
[ well known 


bols, as follows: 


Ou Ou Op re) Ou 
pu + pu = . Hu (1) 
Ox oy Ox oy oy 
Op/oy = VU (2) 
(0/Ox)( pu) + (0/Oy)( pu) = O (3) 
or oT Op . ou \* O or 
IC, | pu + pv = u + u + k (4) 
: Ox oy Ox oy/ oy oy / 


Assume Pr (Prandtl Number) = 1 and JC, = constant, and in 


troduce 
ICgT + ('/2)u* 


= soe — 
TCpT1 + (1/2)u;? 


l (5) 


where suffix 1 refers to the main stream conditions. Then, the 


energy equation, Eq. (4), can be put into the following conven 


Os Os Oo oS 
pli + pv = M (6) 
OX Oy oy oy 


Perform, next, Stewartson’s transformation! from (x,y) plane 


ient form: 


to(Z,)’) plane where 


If we take p T, then Eqs. (1) and (6) become, respectively, 


ory Oy doy O7y oy dv, 
vs — 4 ="— — (7 + 3) (9) 
oY OY OYOZ 0}? OZ dZ 
wWosS aSody oS i 
—_ = (10) 
OY 02 OY OZ oFY* 
where y is the stream function defined as follows: 
po oy a, OY 
“= Vv = (11) 
p oy ay O) 
Po oy 
Vo 
p Ox 
p (* 3y—1 1) oy _ ov o} - 
au ; (12) 
» © EN oZ OY ox 
Vi = (ao/a;)u, = anM, (13) 


Therefore, in the (Z, Y) plane the boundary-layer equations be 


come Eqs. (9) and (10). The usual boundary conditions 


y = 0, u=v=0, T=T, (14) 
yo &, u—> MK, ae 
can be represented as 
ov 
Taw: 4 =y 0 S Su 
oO} ‘ 
(15) 
ov 
ae E age B, so | 
oO} 


where the suffix w refers to the surface conditions. The condi 
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tion of heat insulation at the surface is simply S = S, = 0; in 
that case, Eq. (10) is an identity, while the problem reduces to 
the solution of the incompressible boundary-layer equation with 
a main stream velocity V;.! For the general cases with heat flow 
across the surface, we must, however, consider Eqs. (9) and 
(10), together with the boundary conditions, Eqs. (15). 

If, in Eq. (15), S~ = constant, we can construct similar solu 
tions for the system, Eqs. (9), (10), and (15). To find such simi 
lar solutions of the compressible boundary-layer equations, we 


introduce f, g, 7 as follows: 


y = Z°V;,"f(n) (16) 
S = Swg(n) (17) 
Y = 2°Vi% (18) 


where a, b, p, and g are constant parameters satisfying the fol 


lowing conditions simultaneously : 


dV, a 
: Kee V, “2 (19) 
dZ 
dV, ee ec, , 
+ « Z-%y2—-2b _- Z-1V, (p #0) (20) 
dZ p p 
dV, Kw = a ih ial b ‘ 7. 
| = Z-%4y,2-2? — ZV, (q #0) (21) 
dZ qd q 


The partial differential equations, Eqs. (9) and (10), then become 


ordinary equations 


d*f _ fdf\? _ df ‘ : 
- — Ky + Kof - = —K3(1 + Aug) (22) 
dn? dn dn? 
(d?g dn?) + Kof(dg dn) = 0 (23) 
The boundary conditions, Eqs. (15), become 
n = 0, f = df/dn = 0, g=] (24 
n> ©, df/dn = Kg, g=0) = 
In Eqs. (19)-(24), An» are constants (mn = 1-10). They are re 


lated simply as follows: 


K; Ky _K a i 
K,=—, Ke==> Ko=-=, Kw = Ke — Ky 

Ks Ks Ks (25) 
K, = Sy, Ks; = 0 \ 


By examining the integrals of the differential equations [Eqs 
(19)-(21)], we conclude that the following classes of similar solu 


tions are possible: 


(ljp=1+q, a =hb=0, Vi; = [(1 + 2q)(K3Z 4 C)}} 1+2¢ 
(26)* 
Eqs. (22) and (23) become 
d'f _ af : df \? ; 
+ (1 + g)Kif — 5 = Ai — (1 + Kag) | (22a) 
dn dy? dn 
(d2g/dn?) + (1 + q)Kif(dg/dn) = 0 (23a 
The boundary conditions are 
n = 0 f = df/dn = 0, g= 1 (25 
pore ys | 
n> &, df/dy = 1, g=0 
(2) p=1/2, gq= —-(1/2), a=b=0 Vi = ceX4 (27)* 
Eqs. (22) and (23) become 
d*f l d*f . df\? 2 
T Kif = K;, -(1 + Kyxg) (22?b) 
dn} 2 dn? dn 
(d2g/dn?) + (1/2)Kif(dg/dn) = 0 (23b) 
The boundary conditions are the same as given in Eq. (25a) 


* C is an arbitrary constant 
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. Gg = b l —_ 2h l = 2a . Ky S l l l 1 - 
(3) = = =m, V,= Z™ (28) (4) p=-, q=--, a=-, b=-, Y =CZK 
i-~-q=— p 2g+1 2p-—! m 2 2 2 2 
(31 
(m being a constant ). : is eae 
Eqs. (22) and (23) become Eqs. (22) and (23) become 
; ; : d*f l 5 eer . df \? : 
d5f ey : df \? 2 ; - + -(1 + KA3) f — = Ks; — (1 + Kyg) (22d 
+ Kof = Kh, — Kx(1 + Kg) (22c) dy® 2 dy? dn 
dn? dn? dn 
(d*g/dn*) + (1/2)(1 + K3)f(dg/dn) = 0 (23d 


(d*g/dn?) + K2f(dg/dn) = ( 


bo 
ww 
x 


The boundary conditions are the same as given in Eq. (25a) 
Phe boundary conditions are In cases (3) and (4) above, the expressions for V; are similar, 


usually it is more convenient to deal with case (3), with, indeed, 


9 = 0, f = df/dy = 0, ges (99) little sacrifice in generality. Hence, similar solutions are possible 
i> 2, df/dn = Kg, g=0 if Sw = constant and if a, b, p, q, and V, are defined as in Eqs 


(26), (27), (28), and (31). 
K,, Ke, K;, and Kg are mutually related as follows: For a limited range of values of S,, and Ay, numerical integra 


tion of the ordinary differential equations, Eqs. (22d) and (234d), 
5 1 1 


o— - 1+ [m/(a—b)} A 7 i - _ 
K; = mK with the proper boundary conditions, Eq. (25a), has been carried 


>. Pe .-14m(2p—1)/(e—d)] out or the Reeves Electronic Analogue Computor.t To con 
K, = Ke , ue ; P 
2p —1 clude the present note, we have shown in Fig. 1 a diagram of 
— (30)  2Cy,/C;, versus K; obtained from analyzing these numerical 
K, = op —(b-4 mo | K,i t+ Imp —1) /(a—6)] results. C,, and C;, are the local coefficients of heat transfer 
ae = 1 
Ky\@~-? +1 * C is an arbitrary constant 
K, = t The authors are indebted to Arthur Presson, of the Jet Propulsion Labo 
m ] 


ratory, for his close cooperation in obtaining these numerical solutions 
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READERS 


ind skin friction, respectively, defined as 
ky (OT /OY) w 


Cx ee oe T T 
' Cppit( Tw — Twins.) 


Bw OU/OY) w 
(1/2) pr? (32) 
where Twins. is the temperature on an insulated surface, 
Twins. = T1}1 + [(y — 1)/2] 403} 
It can easily be shown that 


2C,, dg\ / {df 
=— r (33) 
C, dn / dn? } |n=0 : 


The parameter A; is related to the pressure gradient in the free 


stream as follows: 


(37 —-1)/( 1 
a Y Y Y- 
) 1 4. 


ay 


K;= - 


1 dp * fa, 


piu,” dx ay 
AO 


Evidently, A; = 0 corresponds to the case of a flat surface with 


zero pressure gradient, for such case, 2C,,/C,;, = 1 as expected 


For a flat surface with favorable pressure gradient, dp/dx < 0, 
K; > 0, it is found 2C,,/C,;, < 
steeper as S,, (which for the range of S,,’s computed are positive ) 
This seems to indicate that the effect of a favorable 


1; the curves in Fig. 1 become 


increases. 
pressure gradient is to increase the skin-friction coefficient such 
This effect 


becomes more acute when the surface is much heated above the 


that the ratio 2C),,/C,, becomes smaller than unity. 


For a flat surface with adverse pres- 


0, it is found 2C,,/C;, > 1; the 


free-stream temperature. 
sure gradient, dp/dx > 0, K 
negative A; branch of the curve given in Fig. | is computed for 
S, = 1. It is well known that a laminar boundary layer will 
separate under the influence of an adverse pressure gradient of 
For a separating flow, C;, ~ 0, this seems 
curve for A 


certain magnitude 
the reason for the rapid rise of the 2C),/Cy, Qin 
Fig. 1. 
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SUMMARY 


For many reasons, it is desirable to obtain reliable and complete informa 
tion about the turbulence structure of the atmosphere 


correct interpretation of wind-tunnel results on laminar profiles, maximum 


For instance, the 
lift devices, etc., depends upon the knowledge of the structure of the turbu 
lence encountered in flight Also, gust loads and vibrations can be better 
predicted if the turbulence structure of the atmosphere is thoroughly known 

Many flight tests have dealt with the measurement of the turbulence en 
countered in flight. However, it seems that the most versatile method specifi 
cally suited to study of the problem at great flight speed has not been used 
or not sufficiently so, to furnish all the information required. The following 
explains this comparatively simple method and presents results obtained by 


its use 


THE MEASURING METHOD 


pressure, as 


T. RBULENCE CAUSES OSCILLATIONS of the ram 
If (Aé) is the aver 


well as oscillations of the angle of attack. 
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age speed vector of the turbulence, (v) is the flight speed, and 
(p) is the air density, the average amplitude of the ram pressure 
oscillation is 


Ap = pxvxAv (1) 


and the average amplitude of the angle of attack oscillations is 


(2) 


A& = Ai/% 


The equations indicate that, for equal turbulence, the ram-pres 
sure oscillations increase with increasing flight speed, while the 
angle of attack oscillations decrease with increasing flight speed 
Thus a study of the turbulence carried out at high speed must 
concentrate on the measurement of the ram-pressure oscillations 

Today’s flight speeds are close to sonic speed or approximately 
1,000 ft. per sec. 
sity, the unit turbulence speed vector can be expressed as a ram 


When assuming this speed and sea-level den 


pressure noise level: 
Aé = 1 in. per sec. = 114 db. acoustic pressure (3) 
(acoustic reference level = 0.0002 dyn. per cm.?) 

One inch per second is certainly a small speed vector of the 
atmospheric turbulence, and yet this small speed vector generates 
a high noise level when acting on a microphone that is exposed 
The 


pressure of 114 db. mentioned is, for instance, comparable to the 


to the ram-pressure oscillations at high speed. acoustic 
peak pressure generated by a conventional single propeller air 
plane at a 100-ft. distance. Thus a microphone fed by the ram 
pressure is a comparatively simple means to measure the turbu 
lence encountered in flight if the noise generated by the carrier 
airplane does not interfere with the measurement. The elimina 
tion of the carrier noise interference is the greatest problem of the 
method. There exist two basic possibilities to eliminate this 
interference. 

(1) Ifthe carrier speed is slightly below sonic speed, mounting 
the microphone on the tip of a boom protruding in front of the 
carrier is one possibility. The geometric distance between the 
microphone and the noise source will thereby be acoustically in 
creased with the factor |1/(1 — W))] when (1/9) is the speed of 
the carrier expressed in Mach Number. 

(2) If the carrier speed is approximately half the sonic speed, 
the microphone must be mounted inside a towed body that trails 
the carrier at an appropriate distance. Naturally, the weight 
of the towed body and the drag of the towing cable must be so 
that the towed body flies below the turbulent wake of the carrier 


THE RESULTS OF FLIGHT MEASUREMENTS 


A comparatively slow single propeller airplane was used for the 
flight measurements. Its maximum diving speed was limited to 
Mach 0.6. Consequently, the microphone was mounted inside 
a towed body which, in turn, was released and retrieved in flight. 
The towing cable had a length of 400 ft. and an outside diameter 
of '/i¢ in. The overall length of the towed body was 3 ft. and 
its total weight was 12 lbs. 
cycles per sec. and had a band width of approximately 10 per 


The microphone was tuned to 175 


The microphone output was fed 


through the towing The 
voltmeter was calibrated in decibel center frequency pressure at 


cent of its center frequency. 
cable into an electronic voltmeter 
sea level. 

All flight tests were flown under fair weather conditions. The 
first flights were to prove whether a carrier noise interference 
Maximum power flights and minimum power dives 
Since the 


existed. 
were repeated at various flight speeds and altitudes 
power setting did not affect the microphone output, it was felt 
that the carrier noise did not interfere Nevertheless, most of 
the results were obtained at a low power setting. 

The first actual measurements were to determine the relation 
between the the flight 
flights at 10,000 ft. altitude indicated that, within the range from 


150 to 450 m.p.h. true speed, the noise pressure is proportional to 


ram-pressure noise and speed. Two 


the flight speed. An essential conclusion can be drawn from this 
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result. Since the microphone was tuned to a fixed frequency, a 
variation of the flight speed means that the noise is generated by 
eddies of various sizes. For instance, at 150 m.p.h. or 220 ft. per 
sec. flight speed, eddies of 1.26 ft. in diameter generate the 175 
cycles per sec. noise; while at 450 m.p.h. or 660 ft. per sec. flight 
speed, eddies of 3.78 ft. in diameter generate the 175 cycles per 
sec. noise. Thus, according to Eq. (1), proportionality between 
the flight speed and the noise pressure means the average fair 
weather turbulence speed vector (Ab) is independent of the eddy size. 
Naturally, this conclusion is limited to the comparatively narrow 
range of the measurement. Greater speed ratios, as well as the 
use of several center frequencies, could broaden this range and 
might still disclose a variation of the fair-weather turbulence vec- 


tor for greatly different eddy sizes. 
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The next flights were to determine the relation between the 
Fig. 1 compiles the 
results of seven flights at 400 m.p.h. true speed and various alti- 


Four flights covered the range from 1,000 to 10,000 ft 
The line (1) in Fig. 1 indicates the averaged result of 


ram pressure noise and the flight altitude. 


tudes. 
altitude. 
these four flights. 
up to 30,000 ft. altitude. The averaged result of these two flights 
is indicated by the curve (2) in Fig. 1. Both flights were fiown 
Finally, one flight was flown 
Curve (3) in Fig. 1 indicates the 


Two more flights extended the measurement 


on clear windy days in summer 
on a calm hazy day in the fall 
result of this last flight. 
Fig. 1 shows the well-known fact that wind and sunshine cause 
This ground effect dis- 
altitude 


increased turbulence near the ground. 
10,000. ft 


the measurement on a calm hazy day seems to indicate that, out- 


appears above approximately Especially 
side of the ground effect, the ram noise is approximately propor- 
tional to the ambient air density. This, according to Eq. (1), 
would mean the average fair weather turbulence speed vector (di) 
is constant above 10,000 ft. altitude. This is an important con- 


clusion if it can be verified by more extensive flight 


tests. 

Surprising is the actual value of the average turbulence speed 
vector. 
comparing the result with that of Eq. 3, we find the average fair- 
10,000. ft 


When extrapolating curve (3) in Fig. 1 to sea level and 


weather turbulence vector above altitude to be as 


1,000 l : 
1.0 = (0.6 in. per sec (4) 
580 2.83 


CONCLUSIONS 


small as 


(1) A narrow-band microphone fed by the ram pressure was 
used to measure the ambient turbulence in flight. The method is 
comparatively simple and informative, particularly in connection 
with high-speed carrier airplanes 

(2) The results obtained seem to indicate that, under fair- 
weather conditions, the average turbulence speed vector is inde- 
pendent of the eddy size, as well as of the altitude, when measured 
above 10,000 ft The fair-weather turbulence speed 
vector was 0.6 in. per sec. above 10,000 ft Near the 


altitude 
altitude. 


ground, the vector may be up to ten times stronger. 


(3) When using a wide-band microphone and a tape recorder 
aboard a high-speed carrier, a future investigation would cover a 
large frequency range, as well as a wide altitude range, in a single 
flight. Such contribute to the 


knowleage of the turbulence structure of the atmosphere 


investigations would greatly 
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